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NATIONAL  DEVELOPMENTS 


COMPUTER  RETRIEVAL  SYSTEM  PROPOSED  FOR  STANDARDS  DATABASE 

Beijing  ZHONGGUO  BIAOZHUNHUA  [CHINA  STANDARDIZATION]  in  Chinese  No  3,  Mar  86 

pp  19-20 

[Article  by  Yu  Shunzhong  [0151  7311  0022],  Qinghai  Provincial  Bureau  of 
Standards:  "Regarding  Certain  Problems  in  Overall  Design  of  Computer  Retrieval 
Systems  for  Standards"] 

[Text]  There  have  been  precedents  both  within  this  country  emd  outside  it  for 
using  computers  to  retrieve  standards.  Initial  explorations  have  also  brought 
us  certain  results.  To  keep  those  in  this  profession  from  making  the  same 
mistakes  and  from  duplicating  our  efforts,  we  are  sharing  what  little 
experience  we  have  accumulated  and  our  views  on  the  subject,  all  for  the 
reference  of  others  and  to  hear  the  opinions  of  others. 

1 .  The  problem  of  classifying  stzuadards  documents 

From  one  point  of  view,  standards  documents  can  be  classified  into  two 
categories:  one  being  the  standards  documents  themselves  (hereafter.  Just 
"standards")  and  the  other  reference  books  and  periodicals  for  formulating  and 
revising  standards,  as  well  as  writing  materials  (hereafter.  Just  "standards 
materials" ) . 

Because  in  comparison  standards  and  standards  materials  are  characterized  by 
few  categories  and  clear  subclassifications,  they  should  be  handled  specially 
in  classifying  Euid  retrieving.  We  feel  that  standards  should  be  stored 
according  to  issuing  sequence  number  and  that  retrieval  by  category  number  is 
most  appropriate;  for  standards  materials  it  is  the  same  as  with  ordinary 
books,  where  it  is  most  appropriate  to  store  emd  retrieve  by  the  book 
classification  number.  This  is  to  say  that  computer  retrieval  of  standards  is 
retrieval  on  a  small  scale  in  which  it  is  best  not  to  Include  the  retrieval  of 
standards  materials,  especially  when  using  microcomputers. 

2.  The  problem  of  writing  records  and  their  formats 

Since  a  database  is  built  in  a  computer  for  retrieval,  it  may  be  said  that 
more  should  be  available  from  what  is  retrieved  than  when  retrieval  is  manual, 
and  we  have  listed  22  items  for  the  record  table  and  used  the  contents  of 
these  record  tables  for  retrieval.  There  should  be  no  problems  on  non- 
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microcomputers,  and  on  micros  this  can  be  managed  even  if  using  only  floppy 
disks.  However,  because  disks  must  be  exchanged  we  may  not  retrieve  according 
to  free  association  of  logical  expressions.  Therefore,  I  feel  that  it  would 
be  better  when  using  microcomputers  if  databases  are  not  built  directly 
according  to  the  standard  issuing  sequence  number  in  accordance  with  the 
record  table,  but  rather  each  item  is  only  noted  by  the  most  essential 
subordinate  material.  We  have  experimented  with  this  system,  where  for  each 
item  is  only  noted  the  most  essential  information,  such  as  standard  sequential 
number,  the  standard  name,  book  format  and  page  number,  the  location  of  the 
documents,  and  the  date  of  implementation,  as  well  as  special  category 
numbers.  For  example,  GB2867.7-81  cold  stamp  thrusting  movable  punch  shank 
[lengchongmo  mobing  tuirushi  huodong  mobing]  A4-6  (HAN)  81-01  DAJGFBHE,  where 
the  fewest  number  of  items  satisfy  the  search  for  the  most  essential 
information  desired  for  a  standard. 


In  summary,  whether  or  not  the  record  table  is  to  be  used  is  a  problem  worthy 
of  diligent  deliberation,  especially  when  using  microcomputers. 

3.  The  problem  of  formats  for  standard  classification  and  retrieval 

As  we  have  stated  above,  the  classification  framework  for  standards  should 
only  be  undertaken  according  to  the  issuing  unit  and  the  issuing  sequence 
number.  Because  these  two  items  both  appear  in  the  most  conspicuous  positions 
on  the  covers  of  standards  documents,  they  are  most  easily  recognized.  Also, 
within  a  particular  location  the  classification  framework  can  only  be 
undertaken  according  to  a  single  sequence,  for  otherwise  the  arrangement  of 
standards  will  be  disorderly  and  unsystematic.  Many  different  ways  of 
retrieval  may  be  used,  as  for  example  in  the  "microcomputer  standards 
retrieval  system"  that  we  developed,  where  seven  retrieval  methods  are  used: 
stemdard  sequence  number,  special  classification  number,  UDC  number,  freely 
associated  logical  expressions,  key  words  in  items,  first  words,  and 
identifiers.  Actually,  the  number  of  retrieval  methods  is  determined  by  the 
entry  contents  within  a  record,  where  the  actual  sequential  format  of  the 
standard  is  not  Important.  Even  though,  as  we  said  above,  we  do  not  use 
record  tables,  the  three  retrieval  methods  by  the  standard  sequential  number, 
special  classification  number,  and  freely  associated  logical  expressions  may 
also  be  used.  As  to  a  detailed  explanation  of  retrieval  methods,  one  may 
refer  to  other  specialty  papers,  because  here  we  wish  only  to  explain  that  the 
standards  classification  method  and  retrieval  methods  are  mutually  related  and 
should  be  separately  studied  and  used. 


4,  The  problem  of  subclassifying  standards 


As  for  other  methods  of  retrieval,  the  problem  of  subclassifying  standards 
would  not  appear  to  be  too  important,  but  when  searching  according  to  special 
classification  numbers  there  is  simply  no  way  to  proceed  unless  the  standards 
are  subclassified.  Because  the  person  who  is  doing  the  search  can  only  search 
according  to  special  classification  numbers,  that  cannot  be  as  quick  as 
searching  directly  by  the  standard  sequence  number.  The  most  important 
characteristic  when  searching  for  special  classification  numbers  (and  the 
greatest  advantage)  is  that  the  search  can  be  made  by  hierarchical  reduction 
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of  scope.  On  the  computer,  the  screen  hierarchical  prompting  mode  may  be 
used,  only  by  which  may  the  goals  of  speed,  accuracy,  and  completeness  be  met. 

What  is  the  most  suitable  way  to  subclassify?  After  subclassification,  how 
should  numbering  be  done?  These  are  all  problems  that  should  be  explored 
specially.  Here,  I  will  only  mention  the  essentials  of  my  own  view:  there 
should  be  large  scale  distinctions  made  between  the  international  and 
national,  and  between  the  national  and  the  various  ministry  commissions, 
provinces,  municipalities,  and  prefectures,  according  to  issuing  units.  Also, 
progressively  finer  distinctions  should  be  made  according  to  specialized 
classifications  so  that  a  standard  article  at  the  lowest  level  could  be  fully 
displayed  on  only  one  or  two  screens.  Among  classification  numbering  methods 
there  are  currently  the  alphabetic  numbering  system,  the  digital  numbering 
system,  and  the  alphanumeric  numbering  system,  and  the  "international  decimal 
classification  system"  in  broad  use  is  a  digital  plus  carry  classification 
numbering  system.  Seeing  that  there  are  advantages  and  disadvantages  to  the 
various  methods  currently  in  use  in  systems  for  the  computer  retrieval  of 
standards,  where  it  would  be  difficult  to  achieve  the  most  satisfactory 
results  from  the  use  of  just  one,  I  advocate  the  adoption  of  a  new  alphabetic 
with  carry  classification  numbering  method,  which  is  a  "classification 
numbering  method  having  26  carry  positions"  composed  of  the  26  letters  of  the 
English  alphabet  that  I  have  attached  for  the  reader’s  perusal.  By  using  this 
classification  method  various  results  such  as  a  clear  hierarchical  concept, 
sufficient  capacity  at  each  level,  and  easily  recognized  symbols  can  be 
simultaneously  attained. 

5.  The  problem  of  what  type  of  computer  to  use 

Computer  usage  is  primarily  determined  by  time  (speed)  and  space  (internal  and 
external  storage),  as  well  as  by  cost.  Regarding  time  and  space,  non¬ 
microcomputers  will  naturally  be  more  advantageous  than  microcomputers,  but 
regarding  price,  the  microcomputer  will  be  much  less  expensive.  As  for  the 
retrieval  of  books  or  of  general  scientific  and  technical  materials,  due  to 
great  demands  on  capacity,  non-microcomputers  should  generally  be  chosen.  As 
for  the  locations  of  standards  materials  of  a  national  nature,  and  since  the 
total  number  of  world  standards  exceeds  1.2  documents,  if  they  are  to  be 
stored  as  completely  as  possible,  then  the  large  scale  mainframes  must  be 
used.  But  for  a  small  to  medium  size  standards  information  office  at  the 
province  or  city  level,  where  databases  are  built  in  accordance  with  the 
format  of  our  record  tables,  and  where  a  location  stores  from  5  to  10 
standards  documents,  then  use  of  a  microcomputer  is  entirely  possible.  If  an 
IBM-PC/XT  Apple  [unclear  how  this  is  intended]  is  chosen,  with  512K  RAM  and 
fitted  with  a  10  (or  20)  megabyte  hard  disk,  the  total  price  will  be  about 
10,000  yuan;  in  the  case  of  the  domestically  produced  Great  Wall  0520C 
compatible,  where  the  machine  already  has  a  20  megabyte  hard  disk,  another  20 
megabyte  hard  disk  can  be  added  for  a  total  current  price  of  about  40,000 
yuan.  If  our  record  table  is  not  used  and  the  standard’s  sequential  numbering 
is  used  directly,  some  information  will  have  to  be  added  to  the  database,  but 
the  number  of  standards  that  can  be  stored  will  be  greater  by  perhaps  10 
times.  And  at  the  same  time,  when  space  allows,  use  of  a  microc<xnputer  will 
have  a  series  of  advantages  such  as  less  investment,  quick  start  up,  low 
requirements  for  the  operating  environment,  easy  security,  and  ease  in 
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learning.  Therefore,  we  believe  that  one  should  proceed  from  the  actual 
conditions  of  the  total  number  of  standards  needed  to  be  stored  at  a 
particular  location.  Where  a  microcomputer  can  be  used,  it  should  be.  We 
need  not  think  it  necessary  to  have  the  biggest  and  most  perfect. 

6.  Regarding  the  necessity  of  setting  up  auxiliary  functions  other  than 
inquiry 

In  a  standards  retrieval  system,  there  should  be  other  auxiliary  functions 
than  just  the  two  important  ones  of  establishing  databases  and  retrieval, 
functions  such  as  revising,  security,  and  printing. 

Regarding  the  revising  function,  whether  the  error  is  discovered  when  setting 
up  the  database  or  during  its  use,  it  must  be  correctable.  At  this  time  there 
must  first  be  deletion,  then  addition;  there  must  also  be  random  entry  after 
the  Issuance  of  new  standards,  and  this  naturally  requires  extra  space  that 
can  be  used.  If  these  functions  do  not  exist  there  will  be  all  sorts  of 
disorder  in  altering  and  there  will  be  many  errors. 

As  for  the  security  function,  this  is  done  by  setting  up  encrypted  passwords, 
where  the  encryption  can  be  randomly  altered  as  needed  to  prevent  unauthorized 
personnel  from  operating  the  system.  Otherwise,  it  is  possible  that  when  not 
careful  relevant  data  in  the  database  could  be  erased,  or  there  could  be  a 
malfunction  in  the  computer  that  would  render  the  system  unusable. 

As  for  printing,  that  is  for  printing  out  the  information  that  has  been 
retrieved,  because  the  ability  to  print  is  not  only  convenient  for  recording 
that  which  has  been  retrieved,  but  various  management  functions  like  library 
contents,  lists  by  special  topics,  contents  by  date,  and  cards  for  manual 
retrieval  can  be  printed  out  according  to  the  results  of  freely  associated 
local  expressions. 

In  summary,  the  auxiliary  functions  of  editing,  security,  and  printing  should 
also  be  available  in  a  standards  retrieval  system,  for  without  these  functions 
the  system  will  hardly  be  usable  in  a  normal  way. 

12586 
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PHYSICAL  SCIENCES 


MICROCOMPUTER  METHOD  FOR  SIMULATING  TRANSIENT  RESPONSE  OF  LINEAR 
DIFFERENTIAL  SYSTEMS 

Chongqing  CHONGQING  DAXUE  XUEBAO  [JOURNAL  OF  CHONGQING  UNIVERSITY]  in  Chinese 
Vol  8  No  3,  Sep  85  pp  55-61 

[Article  by  Huang  Guangping  [7806  03^2  1627],  Zheng  Tikuan  [6774  7555  1401], 
and  Wang  Airen  [3769  5676  0086],  Department  of  Thermodynamic  Engineering] 

[Text]  The  augmented  matrix  method  is  an  effective 
technique  for  solving  high-order  differential  systems  with 
stiffness,  but  the  existing  literature  on  the  subject 
describes  only  the  single-input  form.  The  present  paper 
describes  a  multiple-input  form.  The  calculation  method 
used  combines  time-step  progression  and  iteration,  achieving 
greater  accuracy  and  greatly  decreasing  the  amount  of 
computation  needed,  so  that  it  is  highly  suitable  for 
numerical  simulations  on  microcomputers. 

I.  Introduction 

System  simulation  generally  requires  the  solution  of  sets  of  differential  and 
algebraic  equations  of  the  form 

X  =  +  5U 

Y  =  CX  +  DU 


where  A,  B,  C  and  D  are  constant  coefficient  matrices.  When  the 
dimensionalities  of  the  equations  are  rather  large  (e.g.  10  dimensions  or 
more)  and  the  stiffness  is  rather  great  (a  factor  of  several  hundred  or  more), 
such  problems  are  rather  difficult  to  solve. 

Analog  computers  are  effective  in  solving  such  systems  of  equations,  but 
because  systems  of  high  dimensionality  require  large  numbers  of  amplifiers, 
large  analog  devices  become  necessary.  In  addition,  the  accuracy  of  analog 
machines  is  low  and  it  is  generally  hard  to  assure  high  accuracy.  As  digital 
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computers  have  developed  and  numerous  high-level  languages  have  appeared  and 
come  into  widespread  use,  digital  simulation  has  increasingly  replaced  analog 
simulation.  For  ease  in  making  the  transition  from  analog  to  digital 
simulation,  numerous  digital  computer  simulation  languages  have  been 
developed,  such  as  DAS,  MIDAS,  MIMIC,  DSL,  PACTOLUS,  CSMP  and  CSSL  [1-4]. 

These  languages  can  be  utilized  directly  by  the  user  without  spending  large 
amounts  of  time  on  programming,  and  they  are  convenient  to  use.  But  the 
simulations  can  be  run  only  on  large  computers  and  minicomputers,  while 
microcomputers  have  been  useless  for  the  purpose. 

The  most  widely  used  technique  is  direct  numerical  integration,  e.g.  the 
Runge-Kutta,  Merson  and  Ji'er  [phonetic]  methods  [5,  7].  But  most  integration 
methods  are  limited  by  the  order  of  magnitude  of  the  smallest  time  constant  in 
the  equations.  A  stable  system  whose  stiffness  is  not  very  great  is  readily 
solved  by  the  above  method.  But  if  the  system  has  high  stiffness  and  high 
dimensionality  or  if  it  is  unstable,  the  usual  integration  methods  such  as  the 
Runge-Kutta  method  may  fail  or  have  a  high  degree  of  error,  and  in  addition 
will  be  very  slow.  Such  problems  are  even  more  difficult  to  solve  on 
microcomputers  because  of  the  short  word  length.  Consequently  microcomputer 
simulation  requires  increased  simulation  accuracy  and  decreased  amounts  of 
computation  in  order  to  increase  simulation  speed. 

Taylor's  augmented  matrix  method  [8]  is  an  effective  method  of  solving  stiff 
linear  differential  systems  on  microcomputers  under  certain  specific  input 
conditions.  But  his  technique  is  only  a  single-input  form,  and  it  too 
requires  large  amounts  of  computation  and  a  rather  large  amount  of  internal 
memory  in  the  case  of  high-dimensionality  systems.  Below  we  expand  the 
augmentation  method  to  multiple  inputs,  subject  to  the  requirement  that  every 
component  of  the  input  vector  is  the  usual  type  of  step,  slope,  acceleration 
or  exponent  input.  In  addition,  some  augmentations  are  possible  when  the 
multiple  inputs  are  not  simultaneously  perturbed.  In  the  case  of  high- 
dimenslonallty,  large  time  scale  problems,  in  order  to  increase  numerical 
efficiency,  the  time-step  progression  method  can  be  used  at  the  beinning  of 
computation,  with  later  replacement  by  the  iteration  method,  and  if  the 
optimum  iteration  time  is  determined  the  amount  of  computation  needed  can  be 
minimized . 

II.  The  Augmentation  Method  for  Multiple-Input  Equations 

Suppose  that  the  system  to  be  simulated  is 

r  X(  /  )  =  ^X(  O  +  5U(  M 

■  Y<  #  )  =  CX(  n  +£)U(  i  )  (  I  ) 

X(o)  =  X. 

where  X  is  n-dimensional ,  Y  is  m-dimensional ,  and  U  is  r-dimensional.  In 
other  words,  the  system  is  a  multiple-input,  multiple-output  system. 
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We  know  from  control  theory  [9]  that  the  solution  of  a  differential  equation 
consists  of  two  parts:  the  free  terms  and  the  forced  terms.  In  numerical 
integration  of  the  forced  terms,  assuring  accuracy  requires  a  rather  large 
amount  of  computation  (see  Ref.  10).  If  homogeneous  equations  are  involved, 
then  the  solution  has  only  free  terms  and  accuracy  can  be  satisfactorily 
guaranteed  by  choosing  a  suitable  number  of  terms  in  the  expansion  of  e  . 
The  principle  of  the  augmented  matrix  method  is  that  the  forced  input  terms 
are  converted  to  a  state  vector  and  nonhomogeneous  equations  are  converted  to 
homogeneous  equations.  As  in  the  single-input  augmentation  matrix  of  Ref.  8, 
when  every  component  of  U  is  a  step,  slope,  acceleration  or  exponent,  the 
augmentation  is  easily  realized. 

Taking  as  an  example  the  case  in  which  the  input  vectors  are  all  constants,  we 
write 

U(  n  =  '<  i  t>0 


where 

K  =  CKi,Ki,’“,Kry 


We  define  the  n+1-th,  n+2-th, . . . ,n+r-th  state  vectors  as 

Xr.,S  t  )  =Mz(0  =«r(0  =A', 


Then, 


O  =  o  J 


x,+i(  0  >=. 

«■+*<  o )  =Ktt 


•  •  ••• 


•••  ••• 


i  )=  o> 


x,,,(  o  )  *'iC, 


Then,  following  augmentation,  the  state  equation,  output  equation  and  initial 
conditions  are 


'kit)  ; 

'  1 

' X  it)  ' 

■  -  ' 

to 

x»*iit)  1 

: 

^  ii  +  i  (  t  ) 

1  t)  \  ' 

i 

+  t  ) 

!  .  ! 

• 

i  *  ! 

o 

o 

• 

• 

• 

;  ^ii  +  r(0  / 

;  . 
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Y(  #>  =  CCiI» 


X(  /  ) 


*.t»<  # ) 


^«+rC  t  ) 


X  <  0  )  ■ 

r  :  .  l 

*.♦1(0) 

"■  V’-' -  , ..  '  • 

K, 

0  ) 

Kt 

• 

•  • 

1  •  .  .  ■ ' 

t 

[:■  'Kr-: 

<S) 


Let  X(,f  =  CJt,4i»  •*»♦»»  •••»  .  Then  the  above  equations  can  be 

simplified  to 


where 


fx  (#)  =3?  X(#) 

Y  (#)  =  c  x(n 

Wo)  »  ^0 


A 


{JALA 

i  0  i  o  , 


is  the  augmented  transfer  matrix, 


c  (;c|b> 


is  the  augmented  observation  matrix 


X 


is  the  augmented  state  vector. 


Many  of  the  systems  that  are  simulated  are  multiple-input,  multiple-output 
models.  Some  of  them  have  numerous  influencing  factors,  so  that  the  input  may 
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have  tens  or  hundreds  of  dimensions.  If  we  use  the  augmentation  techniques 
described  above,  then  the  order  of  the  matrix  will  be  much  greater  after 
augmentation,  which  creates  difficulty  in  calculating  the  matrix  coefficients 
and  takes  up  large  amounts  of  internal  memory.  Actually,  when  a  model  is 
being  investigated  or  a  certain  term  is  being  simulated,  far  from  all  of  the 
input  quantities  will  be  used;  in  addition,  some  input  vectors  will  be  zero. 
In  order  to  decrease  unnecessary  calculations  and  increase  numerical 
efficiency,  we  can  confine  ourselves  to  the  relevant  inputs  in  the 
augmentation  process.  If  the  input  vector  is 

V 


where  Ktt  •••>  is  an  -dimensional  constant  vector 

and  UssOsCOt  0,  •••,  0  D*"  is  an  r  -  r^  =  r^-dimensional  zero  vector. 


and  if  for  the  corresponding  B-matrix, 

B  *  C  Bi  i  Bt  3 


where  is  an  n  X  r^  matrix  formed  by  the  column  vectors  representing  the 
nonzero  inputs  of  U  and  Bg  is  the  n  X  r2  matrix  formed  by  the  remaining 
columns  of  U  after  elimination  of  B^,  then  the  equation  has  the  form 


X  =  AX  +  CBt  I 


..... 


or 


X  =  ^X  +  Bi  (16  ) 


We  now  again  perform  augmentation  by  the  above  method;  here,  A  will  be  an 
augmented  transfer  matrix  of  dimensionality  (n  +  r  )  X  (n  +  r  ),  the 
number  of  elements  in  which  is  smaller  than  that  in  the  augmented  full-input 
matrix  by  (n  +  r)*  -  (n  +  r.)^  =  2(n  +  r.)‘r2  +  r®,  which  greatly 
decreases  the  number  of  unnecessary  calculations. 

III.  Combination  of  the  Time-Step  Progression  and  Iteration  Methods 

If  the  coefficients  of  the  augmented  transfer  matrix  are  exp(Ah),  where 

exPUh)  =  J  +'Ah  +  ^  +  •.*  +  +  ...^ 
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with  h  representing  the  time  step  (which  is  generally  chosen  rather  small  in 
order  to  assure  convergence  of  the  coefficients),  then  the  augmented 
homogeneous  equations  can  be  solved  by  iteration,  obtaining 

X(h)  =  expCAh)X<o) 

X(2h  =  expiAh)Xih) 

•••  ••• 

X(nft)  «  expCAh)  XCin  ~  l)h:i 


and  by  using  the  "hopping"  [kuitiao  6^23  6426]  technique  we  can  obtain  the 
time  step  progression  solution 

)«/»)=  exp(Jh(X(o) 

T(2h)  =  exp(Ah)»exp{Ah)  *X{o) 

=  expi2~Ah)X(.o) 

Xi4h)  =  exp  (.2Z4h)-exP(2Ah)  •  X(o) 

•••  ••• 

Xi2nh)=  exPin'Ah)  •  exp  (nAh)  •  X(  o  ) 


The  first  of  these  systems  is  the  solution  of  the  iteration  problem,  while 
the  second  gives  the  increase  in  the  time  step  in  the  matrix  coefficients 
obtained  as  a  result  of  squaring.  In  the  former  solution  the  iteration  erors 
will  accumulate;  and  although  the  latter  has  no  iteration  errors,  as  the  time 
lengthens,  the  steps  will  gradually  increase,  and  eventually  the  steps  will  be 
too  large  and  the  resulting  plots  will  be  inaccurate.  In  terms  of  the  amount 
of  computation  required,  although  the  time-step  progression  method  involves  a 
small  number  of  computation  steps,  each  step  entails  multiplication  of  one 
matrix  by  another  followed  by  multiplication  by  a  vector,  so  that  in  the  case 
of  a  high-order  matrix  the  total  number  of  calculations  will  be  rather  large. 
Although  in  the  iteration  process  each  step  involves  simply  the  multiplication 
of  a  matrix  by  a  vector,  for  a  large-time-scale  problem  the  amount  of 
computation  required  for  the  multiple  iterations  will  also  be  rather  large. 

An  effective  way  of  overcoming  these  deficiencies  is  to  combine  the  two 
techniques,  using  the  time-step  progression  method  at  the  outset  for  small 
time  scales,  then  changing  to  the  iteration  method  after  the  time  step  reaches 
a  certain  length.  If  the  initial  step  is  h,  the  iteration  interval  is  T,  and 
the  solution  time  scale  is  [0,  T^],  then  the  time  progression  is 
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A,  2A,  4A,  Bh,  i-r,  r,  27’,  37’,  Tx-Ty  Ti 

K 

where  T  =  2  h;  K  is  the  number  of  cumulations  in  the  first  stage.  In 
order  to  minimize  the  amount  of  calculation,  an  optimum  iteration  time  T  can 
be  chosen.  Let  there  be  a  set  of  homogeneous  equations  (obtained  by 
iteration  as  described  above  if  the  problem  initially  involved  non-homogeneous 
equations) : 


i  =  AY. 
Y  =  CX 


where  X  and  I  are  an  n-dimensional  and  an  m-dimensional  vector 
respectively. 

The  number  of  steps  in  the  time-step  progression  is  K  =  loggCT/h).  The  step 
size  in  the  iteration  calculation  is  T./T  -  1.  The  number  of  multiplications 
performed  in  the  time-step  progression ^calculation  is 

i 

Pi  =  (n*  +  n*  +  n  X  m  )•  (052  ^ ) 


The  number  of  multiplications  needed  in  the  iteration  calculation  is 

Pi  =  («*  +  »  X  m  )  (  -p-  -  1  ) 


The  total  number  of  multiplications  performed  in  the  entire  calculation 
(excluding  the  determination  of  exp(Ah))  is 

P  -  Pi  +P*  =  (n*  -hn  +.m)*n.  /o52  ( ^  +  <n +  m)  •«.  ^  -pr--  1  ) 


We  find  the  poles  of  the  above  equation;  then  P  is  at  a  minimum  when 

T  -  ^  )/n2 

“  V  +  r»  +  m 
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Because  T  must  meet  the  condition  T  =  2*^h,  if  we  first  specify  h,  then  the 
time  obtained  from  equation  (7)  may  fall  between  two  values  of  T.  Eithter  the 
number  of  entries  for  the  larger  of  the  T's  can  decreased,  making  the  grid 
somewhat  sparser,  or  the  grid  for  for  the  slightly  smaller  T  can  be  made  more 
dense,  but  in  the  initial  stage  both  cases  are  the  same;  the  only  difference 
is  that  one  involves  more  steps  than  the  other.  As  regards  asymptotic 
stability,  the  former  approach  is  more  suitable;  but  if  the  duration  of  the 
process  is  too  long  and  if  peaks,  oscillations  and  the  like  may  occur  in  the 
latter  part  of  the  curve,  the  smaller  value  of  T  should  be  chosen. 

Conversely,  if  the  optimum  value  of  T  is  chosen  from  the  formula  and  is  used 
in  turn  to  choose  h,  the  best  method  is  to  work  backward  in  a  geometric 
progression,  obtaining  a  value  of  h  that  is  related  to  the  convergence  of 
exp (Ah) . 

Ref.  11  is  an  example  of  a  solution  using  the  method  described  here.  It 
involves  a  model  of  a  200,000  kW  fossil-fired  power  generating  unit  operating 
in  the  vicinity  of  the  rated  characteristics,  with  a  36-dimensional  state 
vector,  a  9-dimensional  input  vector,  a  non-state  vector  requiring  7 
observables,  and  a  dynamic  process  duration  of  3000  seconds.  When  the  input 
vector  is  all  step  values,  the  augmented  matrix  method  is  used  to  convert  the 
non-homogeneous  equations  of  the  numerical  model  into  homogeneous  equations. 
Thus  n  =  36  +  9  =  ^*5,  m  =  7,  and  T.  =  3000  seconds.  We  find  from  equation 
(7)  the  optimum  iteration  step  T  =  52.06  seconds.  In  order  to  make  the 
observations  occur  at  suitable  times,  we  choose  T  =  50  seconds.  From  the 
characteristics  of  the  augmented  transfer  matrix,  in  order  to  make  exp(Ah) 
converge  rapidly,  h  0.1  seconds.  For  T  =  50  seconds,  we  work  backward  in  a 
geometric  progression,  choosing  the  suitable  value  h  =  0.097^^  seconds.  Thus, 
we  expand  from  0.09766  seconds  to  50  seconds  in  a  9-step  progression  and  then 
proceed  to  3000  seconds  in  59  iterations  with  a  50-second  time  step.  In  an 
iteration  process  with  a  time  step  of  0.0766  seconds,  a  total  of  30,719 
iterations  would  be  required  to  reach  T..  If  the  number  of  observations  is 
the  same  in  both  methods,  then  for  the  differential  equations  needed  to  solve 
this  problem  the  number  of  multiplications  required  (not  including 
determination  of  exp(Ah))  would  be  9.5985‘10^  for  the  mixed  method  and 
6.2206*10'  for  the  iteration  method. 

The  latter  method  requires  65  times  as  many  computations.  Another  extremely 
important  point  is  that  because  the  mixed  method  requires  only  59  iterations, 
while  the  iteration  method  requires  30,791  iterations,  the  roundoff  error  of 
the  iteration  method  will  be  greater  than  that  of  the  mixed  method. 

IV.  Conclusions 

The  use  of  the  augmented  matrix  method  to  solve  systems  of  linear  constant 
differential  equations  is  free  of  limitations  on  equation  dimensionality 
and  number  of  inputs  and  can  be  used  to  solve  stable  or  unstable, 
stiff  or  nonstiff  systems.  In  particular,  its  advantages  over  ordinary 
numerical  methods  show  up  most  clearly  when  solving  stiff  systems.  When  the 
combination  of  the  time-step  progression  and  iteration  methods  described 
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above  is  used,  the  computation  is  rapid  and  accurate,  internal  memory 
requirements  are  low,  the  method  is  easy  to  program,  so  that  it  is 
particularly  well  suited  for  rapid  numerical  computations  on  microcomputers. 
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PHYSICAL  SCIENCES 


ADHESIVE  FOR  URANIUM  IN  NATURAL  SEAWATER 

Beijing  HAIYANG  XUEBAO  [ACTA  OCEANOLOGICA  SINICA]  in  Chinese  Vol  7  No  3 
15  May  85  pp  313-316 

[Article  by  Li  Xuanru  [2621  6513  1171],  Chen  Lianzhi  [7115  6647  1807],  Liu 
Shlyu  [0491  4258  3768],  and  Hong  Yuying  [3163  3768  5391];  "Use  of  Alkaline 
Zinc  Carbonate  Adsorbent  for  Adsorption  of  Uranium  in  Natural  Seawater"; 
paper  received  15  November  1983] 

[Text]  The  experimental  results  in  this  paper  show; 

!•  The  following  relationship  exists  between  the  concentration  of  uranium  in 
the  adsorbent  and  the  uranium  contained  in  the  natural  seawater  used; 

C  =  8.51  X  10“l(u*)0-^9. 

2.  When  the  amount  of  adsorbent  basically  does  not  change,  as  the  amount  of 
seawater  passed  through  the  adsorbent  column  Increases,  the  amount  of  uranium 
adsorbed  by  the  adsorbent,  the  uranium  concentration  in  the  adsorbent  and  the 
concentration  of  uranium  which  remains  in  the  natural  seawater  after  uranium 
adsorption  constantly  increases;  the  recovery  rate  of  uranium  constantly 
decreases;  after  the  distribution  coefficient  decreases  to  a  certain  value 
as  the  volume  of  seawater  passed  through  the  adsorbent  column  increases  it 
rises  again  slightly. 

There  are  currently  very  few  papers  abroad  dealing  with  the  adsorption  of 
uranium  in  seawater  using  alkaline  zinc  carbonate.  This  paper  used  an 
experiment  passing  water  through  a  cylindrical  tube  to  bring  the  adsorbent 
directly  into  contact  with  a  large  quantity  of  natural  seawater  and  at  set 
intervals  sampled,  analyzed  and  measured  the  amount  of  uranium  adsorbed  by 
the  adsorbent,  and  from  the  results  of  the  experiment  explain  the  laws  of 
adsorption  of  uranium  in  natural  seawater  using  alkaline  zinc  carbonate. 

Let  u  be  the  amount  of  uranium  adsorbed  by  the  adsorbent  (in  grams) ;  u*  is 
the  amount  of  uranium  contained  in  the  natural  seawater  used  (in  grams); 

G  is  the  weight  of  the  natural  seawater  used  (in  grams)  ;  T  is  the  amount  of 
adsorbent  used  (in  grams) ;  Cq  is  the  concentration  of  uranium  in  the  natural 
seawater  (gram/gram) ;  the  average  concentration  of  uranium  in  seawater  from 
the  Xiamen  harbor  is  2.8  x  10“^  grams  of  uranium/grams  (of  seawater);  C  is 
the  concentration  of  uranium  in  the  adsorbent  (gram/gram) ; 
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amount  of  saturated  uranium  adsorbed  by  the  adsorbent  (gram/gram);  Q  is  the 
amount  of  seawater  used/the  amount  of  adsorbent  used  (gram/gram) ;  R  is  the 
uranium  recovery  rate;  kd  is  the  distribution  coefficient;  k  and  n  are 
constants. 

On  the  basis  of  the  above  assumptions,  we  can  obtain  the  following  series  of 
equations : 


u=CT  : 
u*=CM  , 

Q=G/T.. 

u*-u  CoG-CT 


Co-C/0=C,(l-i?)  ; 


G  G 

C=k,C*=kAC.-C/Q)  ; 


P-JL-SIL _ 

u*  ~  C.G  ~  C,Q  ' 


rC _ ^  _1_ 

'  C*  1-J?^  ~  CjC-f/G  ’ 
Cm.x  =  lim  C=lim  ki{C o  —  C/Q')=kiC 0 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 


I.  Experiment  Part 

(1)  Adsorbent  manufacturing  conditions 

A  specific  quantity  of  0.6M  ZnS04  solution  was  added  to  a  small  quantity  of 
such  organic  additives  as  5-nitro-2~amino  phenol,  and  while  boiling  and 
constantly  stirring  it,  a  l.OM  Na2C03  concentration  of  the  same  mass  as  the 
ZnSO^  concentrations  was  added,  and  the  boiling  continued  for  20  minutes, 
then  it  was  cooled  and  aged  for  22  hours.  The  adsorbent  was  thoroughly 
washed  until  there  were  no  sulphuric  acid  roots,  then  filtered,  dried,  and 
ground  through  a  screen  to  40-80  mu  [4158]  adsorbent. 

(2)  A  water  pump  was  used  to  draw  natural  seawater  into  a  holding  pool,  and 
after  clarification,  it  was  drawn  into  a  sand  filtration  pool  and  filtered, 
then  the  filtered  seawater  was  drawn  into  plastic  foam  filter  and  further 
refined.  The  refined  seawater  was  passed  through  the  adsorption  column  with 
its  flow  velocity  controlled  to  a  set  speed  by  a  rotameter,  brought  in 
direct  contact  with  a  set  amount  of  adsorbent,  and  after  exchange  adsorption, 
the  seawater  was  expelled  from  the  upper  end  of  the  adsorption  column. 

The  amount  of  adsorbent  used  in  this  experiment  was  419  grams,  the  seawater 
flow  velocity  was  8.0  kg/min.  After  analysis  and  measurement,  it  was  deter¬ 
mined  that  the  average  amount  of  uranium  contained  in  Xiamen  harbor  seawater 
was  2.8  yg  uranium/kg  (seawater),  then  each  day  about  32.26  mg  of  uranium 
was  passed  through  the  adsorption  column,  and  each  day  samples  were  taken 
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and  the  amount  of  uranium  adsorbed  by  the  adsorbent  was  analyzed  and  measured. 
The  results  are  given  in  Table  1. 


Table  1.  Relationship  of  the  Number  of  Days  Water  Was  Passed  Through  the 
Adsorbent  (419  g)  and  the  Amount  of  Uranium  It  Adsorbed 

Number  of  days  water  passed  Amount  of  uranium  adsorbed 

_ (days) _  (yig  uranium/grams) _ 


4 

272 

5 

302 

6 

336 

7 

340 

8 

354 

9 

382 

12 

429 

13 

437 

14 

469 

15 

514 

16 

518 

From  the  experimental  results  it  can  be  proved  that  the  alkaline  zinc 
carbonate  adsorption  of  uranium  in  natural  seawater  is  little  influenced 
by  the  seawater  temperature.  Thus,  on  the  basis  of  the  results  of  this 
experiment,  the  relationship  between  the  amount  of  uranium  adsorbed  by  the 
adsorbent  and  the  amount  of  uranium  in  the  seawater  which  was  passed  through 
the  adsorption  column  are  in  Figure  1. 


c 


Figure  1.  Relationship  of  C  and  u* 

From  this  it  is  clear  that  the  alkaline  zinc  carbonate  adsorption  of  uranium 
in  natural  seawater  basically  conforms  to  the  temperature  equation  of 
Freundlich  et  al.,  and  its  adsorption  equation  can  be  expressed  as  follows: 

C  =  k(u*)’^, 

in  which  the  two  constants  k  and  n  can  be  solved  on  the  basis  of  the  data  in 
Table  2  serving  as  the  diagram  of  the  relationship  of  logC  and  logu* 

(Figure  2) . 
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Table  2.  Relationship  of  Number  of  Days  Water  Passed  Over  Adsorbent  and 
G»  Uj  u*  j  Qj  C)  Cqj  Rj  and 


1 

I 

HI 

m 

|H 

R 

(xio-*) 

hi 

(xio®) 

5 

67,60 

419 

126538 

.161280 

137.47 

302 

0.603 

78.46 

500.70 

6 

69.12 

419 

140780 

193536 

165.00 

336 

0.763 

72.76 

440.50 

7 

80,64 

419 

142460 

225792 

192.46 

340 

J,033 

63.09 

329.02 

9 

103.68 

419 

160058 

290304 

247.45 

382 

1.266 

55.13 

304.08 

12 

138.24 

419 

179751 

387072 

329.93 

429 

1.500 

46.44 

286.10 

13 

146,76 

419 

183103 

419328 

367.42 

437 

1.677 

43.67 

277.05 

14 

161.28 

419 

196511 

451584 

384,92 

469 

i.582 

43,51 

269.54 

16 

184.32 

419 

217042 

516096 

439.90 

518 

1.622 

42.05 

319.26 

Note:  T  is  actually  not  constant  because  when  the  amount  of  uranium  con¬ 
tained  in  the  adsorbent  was  analyzed,  each  time  we  needed  to  take  a  sample 
of  about  1  gram  for  analysis,  but  its  Impact  on  the  calculations  in  the 
table  is  not  substantial. 

Key: 

1.  Number  days  water  passed  over  adsorbent 

2.  Grams 

3 .  Grams/ grams 


Figure  2.  Relationship  Diagram  of  log  C  -  log  u* 
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Linear  slope: 


2, 8600 --2, 3700 
6.0000-5.0000 


0.49, 


n 


intercept:  logk  =  -0.07,  k  =  8.51  x  10^1. 

Therefore,  the  equation  is:  C  =  8.51  x  10”! (u*) 0-^9^ 

Carrying  out  test  calculations  on  C  on  the  basis  of  the  experimental  data, 
the  results  are  as  follows: 

C  (original)  272  302  336  340  514  518 

C  (calculated)  272  303  331  357  517  536 

C  (original)  is  experimental  test  data;  C  (calculated)  is  data  calculated 
on  the  basis  of  the  above  equations. 

II.  Results  and  Discussion 

1.  When  the  amount  of  adsorbent  (T)  used  does  not  change,  as  the  amount  of 
water  (G)  which  is  passed  through  the  adsorbent  column  increases,  the  amount 
of  uranium  (u)  which  is  adsorbed  by  the  adsorbent,  the  concentration  of 
uranium  (C)  in  the  adsorbent  and  the  concentration  of  uranium  (C*)  which 
remains  in  the  seawater  after  the  uranium  is  adsorbed  increases;  thus  the 
uranium  recovery  rate  constantly  decreases.  However,  after  the  distribu¬ 
tion  coefficient  (kd)  increases  and  decreases  to  a  certain  point,  it  rises 
again  slightly. 

2.  The  adsorption  of  uranium  in  natural  seawater  by  the  alkaline  zinc 
carbonate  adsorbent  under  the  conditions  in  this  paper  basically  conforms 
to  the  Freundlich  equation: 

C  =  8.51  X 
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SPECULAR  DIRECTION  OF  SOUND  REFLECTION  ON  ROUGH  SURFACE 

Beijing  HAIYANG  YU  HUZAO  [OCEANOLOGIA  ET  LIMNOLOGIA  SINICA]  in  Chinese  Vol  17, 
No  2,  Mar  86  pp  111-116 

[Article  by  Qi  Yirang  [2058  6095  6245],  Lu  Deming  [6424  1795  2494],  Xu 
Liezhong  [1776  3525  1813],  Wei  Xiaoyuan  [5898  1321  3123],  Li  Jiangji  [2621 
3068  Q679] ,  and  Wang  Meifen  [3769  2734  5358],  Shandong  College  of 
Oceanology,  Qingdao:  "The  Sound  Reflection  in  Specular  Direction  on  a 
Rough  Surface";  draft  receipt  date:  28  April  1984] 

[Text]  English  abstract:  The  reflected  sound  in  specular  direction  on  a 
rough  surface  was  Investigated  in  laboratory  wind-driven  tank.  The  prob¬ 
ability  density  function  (PDF)  of  the  rough  surface  deviated  from  the 
Gaussian  function  to  different  extent.  The  measured  results  were  compared 
with  the  results  calculated  hy  the  aid  of  the  Gaussian  function  and  Gram- 

Charlier  series.  There  is  a  linear  relationship  between  In  ^  and  2(,ya)^  in 

the  range  of  small  value  of  yo(,ya  <  1)  where  p  and  pQ  are  the  average  ampli¬ 
tude  of  sound  pressure  reflected  on  rough  surface  and  the  amplitude  of 
sound  pressure  reflected  on  smooth  surface,  respectively  and 

y  =  sin0  where  6  is  the  glancing  angle  and  X  is  the  sound  wave  length. 

The  slope  of  the  linear  relationship  is  smaller  than  that  predicted  by  the 
Eckart  theory.  The  standard  deviation  Ow  measured  by  wave-height  meter  is 
close  to  the  value  obtained  from  our  experimental  relationship  but  larger 
than  that  calculated  by  the  aid  of  the  Gaussian  type  PDF. 

The  reflection  characteristics  of  rough  water  surface  on  sound  waves  is 
determined  by  the  probability  distribution  of  the  rise  and  fall  of  the 
surface.  Eckart^ used  the  probability  density  function  (PDF)  of  rough 
surface  to  calculate  the  average  reflected  sound  pressure  and  hypothesized 
that  the  PDF  was  Gaussian.  But  the  Gaussian  PDF  with  regard  to  wind-driven 
water  surfaces  is  an  approximate  description^.  Glay,  Medwin  and  Wright^ 
carried  out  analog  sea  surface  tests  in  a  water  tank.  They  used  plastic 
spread  on  the  water  surface  to  create  an  air  passage  about  15.2  cm  high  and 

Note:  This  work  was  carried  out  with  the  support  of  Professor  Yang  Youfan 
[2799  2589  2868].  While  measuring  waves.  Comrade  Gao  Shenjnie  [7559  1957 
2588]  of  this  institute's  Department  of  Physical  Oceanography  provided 
enthusiastic  assistance  and  we  express  our  deepest  gratitude. 
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produced  waves  using  a  fan.  Their  work  showed  that  the  average  reflected 
sound  pressure  is  very  sensitive  to  PDF  and  the  very  slight  deviation  of 
the  form  of  the  PDF  with  regard  to  the  Gaussian  function  will  clearly 
influence  the  reflected  sound  pressure. 

We  researched  the  problem  of  reflection  of  wind-driven  water  surface  on 
sound  waves  and  on  the  basis  of  the  average  reflected  sound  waves  measured 
during  the  experiments  we  predict  the  standard  deviation  of  the  undulation 
of  the  water  surface. 

I .  Experiment  Method 

The  experiments  were  carried  out  in  the  wind-driven  water  tank  of  our  insti¬ 
tute's  Dynamics  Laboratory.  The  water  tank  is  120  m  long  and  2  m  wide.  By 
changing  the  rotation  speed  of  the  wind  machine  different  wind  velocities  can 
be  produced  and  thus  obtaining  different  water  surface  situations.  To  elimi¬ 
nate  reflection  waves  from  the  two  ends  of  the  water  tank  we  placed  wave 
elimination  material  at  the  two  ends.  In  the  region  where  the  sound  measure¬ 
ments  were  conducted,  the  two  sides  and  the  bottom  of  the  tank  were  covered 
with  rubber  wedge  acoustical  material.  The  overall  length  of  the  measurement 
area  was  5  m,  and  a  0.5  m-high  damper  was  used  to  divide  the  area  laterally 
into  emission  and  reception  areas.  The  damper  was  also  covered  with  rubber 
wedges  thus  effectively  eliminating  interference  from  direct  sound  waves . 
During  the  measurements,  the  water  depth  was  maintained  at  about  0.75  m  and 
the  top  of  the  damper  was  about  0.25  m  above  the  water  surface.  The  water 
surface  sound  photographing  area  was  in  the  middle  of  the  square  measurement 
area  on  the  damper.  Figure  1  is  a  photograph  of  the  measurement  area. 


Figure  1.  Test  Area 

The  acoustical  emitter  was  made  up  of  an  emission  transducer  and  a  reflector 
cover.  The  transducer  component  was  a  ring  made  of  PZT-5  material  with  an 
outer  diameter  of  20  mm  and  an  inner  diameter  of  10  mm.  The  ring  was  placed 
in  a  shell  of  organic  glass  filled  with  organic  silicone  oil.  The  double¬ 
layered  right-angle  cone  type  reflector  cover  was  made  of  2  mm  thick  aluminum 
plate.  Between  the  two  layers  of  aluminum  plate  was  air.  The  diameter  of 
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the  opening  was  250  mm,  and  the  transducer  was  placed  on  the  symmetrical  axis 
of  the  right-angle  cone.  The  directionally  characteristic  acuteness  of  the 
angle  of  the  acoustical  emitter  was  about  12°.  The  electrical  excitation 
signal  was  a  single  frequency  packed  pulse  signal,  with  a  pulse  width  of 
0.5  ms,  a  repeated  frequency  of  20  Hz,  and  a  packed  signal  frequency 
change  range  of  30-170  kHz. 

The  acoustical  receiver  was  a  hollow  spherical  hydrophone  with  an  outer 
diameter  of  10  mm.  During  the  measurements,  the  hydrophone's  position  had 
to  be  adjusted  so  that  the  smooth  surface  reflected  signals  it  could  receive 
were  maximized.  When  the  glancing  angle  was  30°  and  45°,  the  distance 
between  the  hydrophone  and  the  emitter  was  2.4  m  and  1.5  m  respectively. 

The  received  signal  was  observed  directly  with  an  oscilloscope  and  recorded 
photographically.  Examples  of  smooth  surface  and  rough  surface  reflected 
signal  are  illustrated  in  Figure  2. 


Figure  2.  Photograph  of  Rough  and  Smooth  Water  Surface  Sound  Reflection 
Signal 

Glancing  angle  30°,  packed  signal  frequency  100  kHz.  On  top  is 
the  smooth  water  surface  sound  reflection  signal;  below  is  the 
rough  water  surface  sound  reflection  signal  (water  surface 
undulation  standard  deviation  0=5  mm) 

Measurement  precision  using  a  tantalum  capacitor  wave-height  device  to 
measure  water  surface  undulations  was  ±0.3  mm,  when  the  measurement  probe 
was  within  20  cm  of  the  sound  photographing  area.  The  measurement  system 
was  corrected  by  using  the  method  of  gradually  submerging  the  probe  by 
changing  the  height  of  the  water  surface.  The  acoustical  signal  and  water 
wave  signal  were  recorded  at  the  same  time,  taking  about  1  minute  of  record 
ing  time  for  each  record.  The  water  wave  record  curve  was  sampled  at 
intervals  of  0.025  seconds,  the  wave  surface  coordinate  read  and  then  they 
were  processed  using  a  computer.  Figure  3  is  a  block  diagram  of  the 
measurement  system. 

II.  Measurement  Results  and  Discussion 
1.  Wave  Measurement 

On  the  basis  of  the  wave  surface  measurement  data  we  computed  the  standard 
deviation  a,  bias  coefficient  Cg  and  the  peak  coefficient 
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Figure  3.  Block  Diagram  of  Sound  Signal  and  Wave  Measurement  System 
l--Sound  emitter;  2 — ^Water  listening  device;  3— Wave  height 
meter  probe;  4 — ^Emitter;  5 — Oscilloscope;  6 — ^Wave  height 
meter;  7 — Photographic  system.  (The  circle  in  the  diagram  is  2) 
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in  which  N  is  the  data  where  N  >  1600;  is  the  wave  surface  coordinate 
reading}  where  tends  downward  it  is  positive  (see  Figure  3).  H  ±s  the 
average  water  surface  position.  Figure  4a  is  an  example  of  the  wave  record 
curve.  Figure  4b  is  the  corresponding  wave  surface  coordinate  probability 
distribution,  the  experiment  points  are  the  results  of  wave  height  measure¬ 
ment,  and  the  Gaussian  probability  distribution  curve  is  also  drawn.  From 
the  figures  it  can  be  seen  that  the  actual  wave  surface  coordinate  prob¬ 
ability  distribution  is  clearly  different  from  the  Gaussian  curve.  The  peak 
of  the  former  is  near  Z/a  =  0.5. 

2.  Acoustical  Signal  Measurement 

We  calculated  the  ratio  of  rough  water  surface  average  reflected  sound 
pressure  "p  and  the  smooth  water  surface  reflected  sound  amplitude  po  on  the 
basis  of  oscilloscope  photographs  of  the  rough  and  smooth  water  surface 
reflected  signals  under  the  same  conditions.  The  relationship  of  InCp/po) 
and  2(yo)2  is  as  Illustrated  In  Figures  5  and  6,  in  which  y  =  k  sin  0; 

k  =  0  is  the  glancing  angle;  X  is  the  sound  wave  wavelength  and  the 

velocity  of  sound  is  1500  m/s.  When  measured  at  different  frequencies. 
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since  it  is  very  difficult  to  make  the  water  surface  conditions  completely 
identical,  for  each  frequency  we  recorded  the  corresponding  wave  curve  and 
found  the  corresponding  value  of  a.  Within  the  range  of  yo  <  1,  the  experi¬ 
mental  results  can  be  expressed  approximately  as  a  straight  line: 


=“  —  2A(Y<jy 


(4) 


in  which  A  =  0.34. 


W(Z) 


Figure  4.  Wave  Record  Curve  Example  and  Its  Corresponding  Wave  Surface 
Coordinate  Frequency  Distribution 

a — Wave  record  curve;  b — Wave  surface  coordinate  frequency 
distribution 

•••  is  the  wave  height  meter  measurement  results,  measurement  con¬ 
ditions:  a  =  0.44  cm,  Cs  =  0.24,  Cg  =  2.6;  -  Gaussian  frequency 

distribution  curve,  for  Gaussian  distribution:  Cg  =  0,  Cg  =  3 

On  the  basis  of  Eckart's  theory^ 

p  =  Fo  j"^exp(-2/rZ)Pf'(ZVZ  (5) 


in  which  W(Z)  is  the  probability  density  function  of  a  random  surface; 

Ott 

Y  =  k  sin  0,  k  =  — ,  0  is  the  glancing  angle,  and  X  is  the  sound  wave  wave- 
X 

length.  For  Gaussian  surfaces 


^/ 7-jra  2(7  ) 


2}c(t 

Substituting  Eq.  (6)  in  Eq.  (5),  we  get 

ln(p/Po)  “  —  2(ray 

and  in  Figures  5  and  6,  Eq.  (7)  is  expressed  as  a  dotted  line. 


(6) 
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Huang  and  Long  point  out^  that  the  first  four  terms  of  the  Gram-Charlier 
series  can  be  used  to  express  the  probability  density  function  of  random 
surfaces: 


in  which  m3  =  Cs/2;  1114  **  (Ce-3)/2.  And  Hp  is  the  Hermitian  multi-term 
expression: 

3/;  + 

H»  —  /*  —  15/'  +  45/*  —  15. 


In  which  t  »  Z/o. 

Substituting  Eq.  (8)  in  Eq.  (4),  we  get 

i.  «  (l+±  ^  mfrV  -  m,rv)  (9) 

P,  V  3  9  3  / 

whose  amplitude  ratio  is 

Jli  -  I  +  ±  ^ 

Po  '•  3  9 

_  ill  0  0  ) 

27  81  / 

The  results  of  the  calculations  according  to  Eq.  (10)  are  plotted  in 
Figures  5  and  6»  with  Cs  and  Ce  taken  from  Table  1. 


Figure  5.  Relationship  of  ln(p/po)  and  2(yo)2  When  Glancing  Angle  is  30® 

•••  is  test  result,  each  test  point  is  the  average  of  100  pieces 

of  data;  -  represents  Eq.  (7) ;  000  is  the  results  of  calcula- 

2ir 

tlons  using  Eq.  (10);  y  “  ksin0  (in  which  k  =  — *  6  the  glancing 
angle;  X  is  the  sound  wave  long-wave) 
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Figure  6.  Relationship  of  ln(p/po)  and  2(yo)^  When  Glancing  Angle  is  45“ 
(Examples  in  the  illustration  are  the  same  as  in  Figure  5) 

From  Figures  5  and  6  it  is  clear  that  the  actual  measurement  of  the  theoreti¬ 
cal  ratio  of  rough  water  surface  reflected  sound  pressure  is  small.  This  is 
primarily  due  to  the  fact  that  the  wave  surface  coordinate  probability 
distribution  tends  to  be  a  Gaussian  distribution.  Below  it  can  also  be  seen 
that  calculated  according  to  the  Gaussian  probability  density  function 

is  smaller  than  the  Owave  measured  by  the  wave  height  meter. 

3.  Prediction  of  a 

On  the  basis  of  the  relationship  of  InCp/p^)  and  2(Ya)  within  the  range 
YO  <  1»  acoustical  methods  can  be  used  to  measure  the  standard  deviation  of 
the  wave  surface  undulations.  Table  1  gives  the  reflected  sound  pressure 
measured  in  the  experiment  and  the  a  and  crjiigh  ‘^^Iculated  using  Eqs.  (4)  and 
(7)  and  compares  them  with  the  (fwave  measurea  using  the  wave  height  meter. 

Table  1.  cr^ave  Measured  by  the  Wave  Height  Meter  and  the  p/po  Based  on  the 


Measurement  and  a  and  Ohigh  Calculated  Using  Eqs.  (4)  and  (7) 


Key; 

1.  Glancing  angle  3.  <Jwave 

2.  Sound  wave  frequency  4.  Ohigh 
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From  Table  1  it  is  clear  that  the  Ohigh  calculated  using  the  hypothesis  of 
the  Gaussian  probability  density  function  is  generally  larger  than  the 
offset  of  cTwave  the  a  measured  using  acoustical  methods  is  rather  close 
to 

III.  Conclusion 

1.  The  wave  surface  coordinate  probability  distribution  of  rough  water 
surface  clearly  influences  the  reflected  sound  pressure.  In  our  experi¬ 
ments,  the  glancing  angle  was  larger  and  the  shielding  £ffect  should  not 
play  an  important  role.  The  difference  between  the  ln(p/po) ~  2(yo)^ 
relationship  obtained  using  the  Gaussian  probability  density  function  and 
the  experimental  results  is  primarily  due  to  the  fact  that  the  actual  wave 
surface  coordinate  probability  distribution  tends  to  be  a  Gaussian  distri¬ 
bution. 

2.  The  standard  deviation  of  the  measured  reflected  sound  pressure  and 
the  wave  surface  undulation  calculated  using  the  Gaussian  probability 
density  function  is  clearly  smaller  than  the  value  measured  using  the  wave 
height  meter  but  the  value  of  a  obtained  by  calculation  according  to  Eq.  (4) 
is  basically  the  same  as  the  latter.  When  the  acoustical  methods  are  used 
to  measure  standard  deviation,  the  corresponding  sound  reflecting  charac¬ 
teristics  of  different  water  surface  situations  should  first  be  found. 
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CHARACTERISTICS  OF  1.3  MICRON  InGaAsP/lnP  DC-PBH  LASER 

Beijing  BANDAOTI  XUEBAO  [CHINESE  JOURNAL  OF  SEMICONDUCTORS]  in  Chinese  Vol  7, 
No  3,  May  86  pp  324-326 

[Article  by  Zhao  Songshan  [6392  1529  1472],  Wang  Dechao  [3769  1795  6389], 

Wu  Youyu  [0702  0645  1342],  and  Wang  Yuzhang  [3769  3786  4545]  of  Wuhan 
Research  Institute  of  Post  and  Teleconmunication:  "A  1.3  Micron  InGaAsP/InP 
Double-Channel  Buried  Heterostructure  Laser";  paper  received  22  April  1985] 

[Text]  Abstract;  A  1.3  micron  InGaAsP/lnP  double  channel  planar  buried 
heterojunction  laser  device  (DC-PBHLD)  was  fabricated  by  a  super-cooling 
growth  technique  involving  two  liquid  phase  epitaxial  (LPE)  steps.  The 
minimum  threshold  current  at  room  temperature  is  15  mA  and  the  typical  value 
is  20  mA.  The  maximum  continuous  working  (CW)  temperature  is  80“C.  The 
output  power  is  2  mW.  At  4  times  the  threshold  current,  a  stable  longitudi¬ 
nal  mode  output  could  still  be  maintained. 

The  double  channel  planar  buried  hetero junction  laser  device  (DC-PBH  LD)^ 
has  an  excellent  current  limiting  P-N-P-N  structure  which  is  suited  for 
continuous  operation  at  low  threshold  current,  high  quantum  efficiency  and 
high  temperature. 

The  process  flow  diagram  of  this  DC-PBH  laser  is  shown  in  Figure  1  and  its 
structure  is  shown  in  Figure  1(e). 
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Figure  1.  Process  Flow  of  DC-PBH  Laser 

We  employed  a  super-cooling  growth  technique.  The  first  LPE  step  is  to  grow 
three  layers  on  a  Sn-InP  (100)  substrate:  N-InP  buffer  layer  (l-2xl0l8  cm“3) , 
undoped  InGaAsP  source  layer  and  P-InP  (2-5x10^7  cm“3)  overlay.  The  wide 
contact  threshold  current  density  is  approximately  0.7-2  KA/cm^  and  the 
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lasing  efficiency  is  over  50  percent.  Si02  is  grown  on  top  of  the  first  LPE 
surface.  The  surface  is  photollthographlcally  etched  along  the  <011'> 
direction.  A  2  ym  wide  2.5  ym  deep  platform  was  created  by  etching  with  HCl 
and  bromomethyl  alcohol2»3.  On  either  side  of  the  platform,  there  is  a  7  ym 
wide  smooth  bottom  channel.  The  second  LPE  step  grows  four  layers:  P-InP 
(1-3x1017  cm-3)  and  N“-InP  (l-5xl0l7  cm-3)  blocking  layers,  P-InP  overlay 
and  P+-InGaAsP  (more  than  lO^S  cm“3)  contact  layer.  A  conventional  ohmic 
contact  technique  was  used.  The  P  side  is  facing  down  and  is  soldered  on 
a  gold-plated  copper  heat  sink. 

The  characteristics  of  this  DC-PBH  laser  device  are  as  follows: 

1.  Photoelectric  Characteristics 

The  minimum  threshold  current  at  20®C  is  15  mA,  the  typical  value  is  20  mA 
and  the  linear  laser  power  output  is  2  mW;  as  shown  in  Figure  2. 


Figure  2.  Light  Intensity  vs.  Current 
2.  V-A  Characteristics 

The  positive  voltage  drop  at  1  mA  is  0.8  volt.  The  minimum  serial  resistance 
is  Rs  =  2.7fi  and  it  is  typically  3-40.  The  reverse  voltage  drop  at  1  yA  is 
over  2  volts. 
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3.  Spectral  Characteristics 


Figure  3  shows  the  laser  spectra  of  the  DC-PBH  laser  fabricated  at  different 
power^output  levels.  When  the  injection  current  is  between  1.3  ^  and 

4  ^  (i.e.,  power  output  between  1  and  12  mW) ,  it  maintains  a  purg  longi¬ 

tudinal  lasing  mode.  The  half  width  of  the  spectrum  is  less  than  2  A. 


Figure  3.  Spectra  of  DC-PBH  Laser  at  Different  Output  Power 
Levels.  Power  Output  1-12  mW,  Pure  Longitudinal 
Mode. 

4.  Far  Field  Distribution 

When  the  width  of  the  DC-PBH  laser  is  1.5-2  pm,  the  far  field  angle  0^  and  0^^ 
are  almost  equal.  The  light  spot  is  circular.  Figure  4  shows  the  far  field 
distribution  parallel  to  the  junction  surface  at  1-10  mW.  0^,  13.5“ .  When 

coupling  to  an  0.2  aperture  50  pm  diameter  multi-mode  optical  fiber,  the 
coupling  efficiency  reaches  90  percent. 

If  the  width  is  less  than  1  pm,  the  behavior  generally  becomes  poorer.  The 
parallel  and  perpendicular  far  field  angles  0//  and  0.^  increase  and  become 
unsmooth.  In  a  few  cases,  double  peak  spectra  were  observed. 

The  authors  wish  to  express  their  gratitude  to  Comrades  Xu  Guizhen  [6079  2710 
3791],  Song  Xlanzeng  [1345  2009  2582],  Lin  Liren  [2651  4409  0088],  and  Wang 
Suqin  [3769  4790  3830]  for  their  assistance  and  support  in  the  investigation. 

The  authors  also  wish  to  thank  Comrades  Wang  Qimin  [3769  0796  3046]  and  Peng 
Huaide  [1756  2037  1795]  for  reviewing  this  manuscript. 
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Figure  4.  Distribution  of  DC-PBH  Laser  Parallel  to  the  Junction  (T  =  20®C) 
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CHARACTERISTICS  OF  1.3  MICRON  TRANSVERSE  InGaAsP/lnP  BH  LASER 

Beijing  BANDAOTI  XUEBAO  [CHINESE  JOURNAL  OF  SEMICONDUCTORS]  in  Chinese  Vol  7, 
No  3,  May  86  pp  337-339 

[Article  by  Wang  Wei,  Zhang  Jingyuan,  Tian  Huiliang,  Wang  Xiaojie,  and  the 
Device  Technology  Group  at  Institute  of  Semiconductors,  Chinese  Academy  of 
Sciences:  "A  1.3  Micron  Low  Threshold  High  Power  Fundamental  Transverse  Mode 
BH  InGaAsP/lnP  Laser";  paper  received  18  December  1985] 

[Text]  Abstract:  By  using  a  two-step  liquid  phase  epitaxial  (LPE)  process, 
we  successfully  developed  a  1.3  micron  low  threshold  current,  high  power  out¬ 
put,  fundamental  transverse  mode  InGaAsP/lnP  BH  laser.  Its  CW  threshold 
current  at  room  temperature  is  10  mA  and  its  single  side  differential  quantum 
yield  is  31  percent.  The  maximum  linear  power  output  is  20  mW*/facet.  It 
can  operate  in  a  stable  transverse  mode  at  2.5  times  threshold  current. 

1.  Device  Structure 

The  structure  of  the  BH  laser  device  is  shown  in  Figure  1.  The  source  is 
buried  in  InP,  which  has  a  lower  index  of  refraction  among  InGaAsP  materials, 
to  control  the  transverse  model,  xn  the  buried  area,  an  InP  junction  in  the 
reverse  direction  was  prepared.  A  high  impedance  P  layer  is  in  contact  with 
the  source  in  order  to  have  a  good  control  of  the  current.  The  cross- 
sectional  micrograph  after  the  second  epitaxy  is  shown  in  Figure  2.  The 
source  area  of  the  laser  is  0.3  ym  thick  and  2  ym  wide.  The  cavity  is  200-250 
micron  long.  The  preparation  of  the  buried  zone  was  reported  in  reference  2. 


Figure  1.  Schematic  Diagram  of  the  InGaAsP/lnP  Laser  Structure 


*The  power  output  was  determined  by  the  Chinese  Academy  of  Measurements. 
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Figure  2.  Cross-Sectional  Micrograph  of  the  BH  InGaAsP/lnP  Laser 
2 .  Device  Characteristics 

The  minimum  CW  threshold  current  of  the  1.3  micron  BH  InGaAsP/lnP  laser  under 
development  is  15  mA.  Typically,  it  is  20-30  mA.  The  maximum  linear  power 
output  is  over  20  mW/facet.  The  maximum  single  face  differential  quantum 
yield  is  31  percent  and  it  is  typically  20  percent.  The  laser  wavelength  is 
1.29-1.31  micron.  The  intensity  vs.  current  curves  at  various  ambient 
temperatures  are  shown  in  Figure  3.  At  100®C,  it  still  has  a  power  output 
near  3  mW.  Figure  4  shows  the  near  range  distribution  of  the  device  at 
different  operating  currents.  It  shows  that  the  device  can  still  operate  in 
a  stable  transverse  mode  at  2.5  times  the  threshold  current. 


Figure  3.  Light  Intensity  vs.  Current  at  Different  (Ambient)  Temperatures 
3 .  Conclusions 

Recent  experimental  results  show  that  our  technique  is  repeatable.  The  yield 
of  this  low  threshold  device  is  high.  Figure  5  shows  the  statistic  CW 
threshold  current  distribution  at  20"C  of  devices  (256  in  total  which  is 


55  percent  of  the  446  core  tubes  made  out  of  an  epitaxial  wafer  (BH-135)). 
Our  devices  were  aged  at  constant  power  at  50®C.  Some  have  already  passed 
700  hours  of  test*  equivalent  to  10^  hours  of  life  at  room  temperature^*^. 


Figure  4.  Near  Range  Distribution 
at  Different  Operating 
Currents 


Figure  5.  Statistical  Distribution  of 
Threshold  Current  at  20“ C  of 
Devices  Hade  From  BH-135 
Epitaxial  Wafer 
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hermitian  operator  in  tensor  product  space  discussed 

Shanghai  ZIRAN  ZAZHI  [NATURE  JOURNAL]  in  Chinese  Vol  9.  No  5.  May  86  p  394 
[Article  by  Hu  Shuan  [5170  6615  13441  PLA  p  n 

3  May  1985]  Engineering  College;  Paper  received 

[Text]  V  is  an  n-dlmenslonal  unitary  space.  And  is  thp  vn, 

duct  space  with  the  induced  inner  prLu^t  defl^^d  o!  Lco.ptabL  tenso-^’^^" 

JC®-;fj0...0jf*,y®-y,(g)...0;,^by  the  formula  n 

K.  Fan  and  MarcusU,  2]  ^ave  defined  the  numerical  range  of  c<  the  Une^r- 
operator  on®*^,  as  follows:  ®  linear 


WHSf)  -{(SJ'x®,  |;f„  o.  n.}, 

and  studied  a  number  of  its  basic  properties. 

Recently.  Wang  Boylng  (3769  0130  5391]  131  has  provad  that  If  A,  ^.6 

MK),  I— I,  — ,  k,  k<n,  thgp  ^  Jg  hermitian  If  when  k  =  n  =  7 

eti-t- rr 

We  proved  Wang's  conjecture  by  using  the  following  lemmas: 

^  non-zero,  non-scalar  complex  matrices 

and  arri!oi!rerr“°’'’'  ““  coefficients  of  0*AU.  i*W. 


Lemma  2.  LetS'-a4,0^,0^3_g<^»0^.0^^.^  Then 

^i€L(V),  dirnFats. 

It  is  not  difficult  to  prove  the  following  by 
Theorem  1.  Let  0a4,-^,*0...0 

lfsr-0- 


-0,  if  S^t-o  in  which 
Induction. 

which  i4,6l(K),  «>3  .  Then 
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Theorem  2.  Let  In  which  -^16  £(F),  «>3.  Then  ffn  (^£R  If  £^; 

Is  hermitian. 
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COMPUTER  SIMULATION  OF  LiF-CsCI  SYSTEM  BY  MONTE  CARLO  METHOD 

Shanghai  ZIRAN  ZAZHI  [NATURE  JOURNAL]  in  Chinese  Vol  9,  No  5,  May  86  p  395 

[Article  by  Chen  Nianyi  [7115  1819  6318]  and  Xu  Chi  [1776  7459],  Shanghai 
Metallurgy  Institute,  Chinese  Academy  of  Sciences,  and  Jiang  Nalxiong  [3068 
0035  7160],  Shanghai  Computational  Technology  Institute;  Paper  received 
7  October  1985] 

[Text]  The  Monte  Carlo  method  is  a  very  effective  approach  to  study  the  struc¬ 
ture  of  liquids  and  non-crystalline  solids.  Because  ion  pair  potentials  can  be 
easily  derived  accurately,  the  computer  simulation  of  ionic  system  is  particu¬ 
larly  successful.  Though  the  simulation  of  two-dimensional  system  differs  from 
the  real  system,  it  is  valuable  because  its  calculations  are  manageable  and  its 
results  can  be  displayed  on  a  plane  so  that  changes  can  be  observed  directly. 

We  have  used  the  Monte  Carlo  method  to  simulate  64  two-dimensional  ionic  LiF- 
CsCl  systems  and  observed  their  liquid  structures  and  the  process  of  the  forma¬ 
tion  of  "non-crystalline  solid"  by  quick  cooling.  The  Cs"*",  Li+,  Cl~,  and  F“ 
ions  (treated  as  charged  spheres  with  Pauling  radii),  16  each,  form  an  "ion 
cluster"  with  free  boundary.  Using  the  metropolis  calculation  method, the 
system  at  1,300K  reaches  equilibrium  after  105  operations.  There  is  an  obvious 
order  in  the  distribution  of  ions  at  equilibrium:  it  can  be  roughly  described 

as  a  mixture  of  the  Li^F“  and  CsnCV  clusters.  Most  of  the  LinFm  clusters  are 

tightly  packed  with  little  space  between  Li"*"  and  F".  The  Cs^Clm  clusters  are 

relatively  loosely  packed  with  many  holes  inside.  Individual  ion  clusters  are 
not  paired  and  there  are  many  holes  among  them.  Therefore,  their  free  volume 
distribution  is  very  uneven.  At  1,300K,  the  total  potential  energy  of  the 
system  is  -719.08  kJ/mol. 

Based  on  these,  the  process  of  the  formation  of  non-crystalline  solids  by  the 
quick  cooling  of  liquids  was  simulated  by  computer.  When  the  "liquid"  reached 
equilibrium  at  1,300K,  the  temperature  term  (TK)  in  the  Boltzmann  factor  was 
substituted  with  ZOOK  and  the  calculations  were  continued  to  simulate  the  quick 
cooling.  The  calculations  show  that  the  arrangement  of  ion  changes  rapidly 
when  the  temperature  drops.  After  only  170  operations,  the  total  energy  drops 
from  -719.08  to  -750.73  kJ/mol  and  to  -755.24  after  10,030  operations.  After 
that,  the  change  slows  down  rapidly.  After  72,420  operations,  it  drops  to 
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-759.34  kJ/mol.  At  the  same  time,  the  "area”  of  the  ion  clusters  shrinks  by 
about  10  percent  and  it  becomes  a  "non-crystalline  solid"  with  extremely  slow 
relaxation. 

It  is  generally  believed  that  the  structure  of  non-crystalline  solids  is  the 
structure  of  liquid  frozen  by  quick  cooling.  But  actually  only  the  "activation 
process"  (e.g.,  the  disruption  of  the  existing  clusters  in  liquids)  that  in¬ 
creases  the  total  energy  is  frozen  during  quick  cooling  because  the  Boltzmann 
factor  decreases.  On  the  other  hand,  processes  that  cause  a  drop  in  total 
0^0i-gy  ,  opposite  ions  move  closer  to  each  other  within  a  cluster  to 

eliminate  spaces  between  them)  to  speed  up.  The  result  is  that  the  unpaired 
shrinking  ion  clusters  are  packed  into  a  non— crystalline  solid  that  lacks  a 
long-range  order. 


FOOTNOTES 

I,  L.V.  Woodcock  and  K.  Singer,  TRANS.  FARAD.  SOC.,  Vol  67,  1971  p  12. 
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INJURY  THRESHOLD  OF  SKIN  IRRADIATED  WITH  308-NM  EXCIMER  LASER  LIGHT 

Shanghai  ZHONGGUO  JIGUANG  [CHINESE  JOURNAL  OF  LASERS]  in  Chinese  Vol  12  No  12 
Dec  83  p  4 

[Article  by  Li  Zhaozhang  [2621  0340  3864],  Wu  Jianu  [0702  1367  1166]  and  Gai 
tookang  [5556  1405  1660],  Laboratory  of  Laser  Medicine,  Shanghai  No  2  Medical 

College,  and  Chou  Zhengzhuo  [0719  2398  0587],  Shanghai  Institute  of  Laser 
Technology] 

[Text]  ABSTRACT.  Preliminary  results  of  the  use  of  an  XeCl 
excimer  laser  to  study  ultraviolet  erythema  of  the  skin  are 
reported.  Using  white  pigs  as  test  animals,  it  was  found 
that  the  minimum  reddening  dose  (MRD_q)  was  45.97-63.05 
mJ/cm®  (95%  confidence  range).  ^ 

Biological  research  on  the  effect  of  ultraviolet  irradiation  of  the  skin  began 
in  the  early  20th  century.  The  Danish  physicist  Niels  Finsen  was  the  first  to 
use  an  artifical  ultraviolet  light  source  to  treat  patients  suffering  from  the 
skin  disease  lupus.  The  West  Germans  Hausser  and  Vahle  used  a  monochromatic 
light  source  consisting  of  a  low-pressure  mercury  lamp  and  two  quartz  prisms 
for  a  detailed  investigation  of  the  effect  of  ultraviolet  light  of  various 
wavelengths  on  human  skin;  they  noted  that  erythema  of  the  skin  was  produced 
most  readily  in  the  vicinity  of  300  nm.  In  the  1930's  Coblentz  et  al. 
measured  the  ultraviolet  erythema  producing  spectrum  and  found  that 
ultraviolet  light  between  200  and  280  nm  produced  distinct  instantaneous 
erythema,  while  light  at  254  and  297  nm  most  readily  produced  delayed 
erythema.  In  1965,  Everett  et  al.  used  a  monochromatic  source  consisting  of  a 
xenon  arc  lamp  and  an  optical  grating  to  determine  the  erythema  effect  curve 
for  humans  in  the  range  between  250  and  310  nm.  They  discovered  that  as  the 
wavelength  of  the  ultraviolet  light  increased,  the  minimum  amount  of  energy 
needed  to  produce  erythema  also  rose.  When  the  wavelength  was  above  300  nm, 
the  erythema-producing  effect  clearly  decreased.  In  the  1970 's,  Parrish  used 
a  molecular  nitrogen  laser  to  irradiate  the  skin  of  white  subjects  and 
determined  that  the  minimum  exposure  from  a  pulsed  Ng  laser  that  produced 
erythema  was  21.6±4.8  J/cm®.  Aufmuth  (1979)  reported  the  effect  of  a  dye 
laser  (tunable  between  260  and  345  nm)  on  human  skin;  he  determined  that  the 
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Ught  wa=  5 

pulse  uUravlouriSer  Ueit  „f h  '•*=“'•'='>  O"  fe  effect  of  short- 

induced  chemical  reactlons*Dlav  anTmLrt*'’t  ’?*>',  discovered  that  light- 
reaction.  '•actions  play  an  Important  role  In  the  ultraviolet  erythLa 

JaCthre's^oM^rt'h;  SoL're''ors"irto°™5®"'  ‘he 

Besearoh  in  this  area  is  much  n»Lf5  ^  308-nm  exclmer  laser  light, 
standards  for  ChC!  User  safety 


1.  Materials  and  Method. 

?L';uu“  ia:V:s'  aT?h“°"  “U 

mJ;  the  dispersion  angirwai  Ig^g  ^  T"  =‘"?l'-?“Ue  output  was  100 
mrad  In  the  Lrt-axL  dlr^tlon  Ung-axls  direction  and  7.1 

2peSSs,°ihrJrradlUon^IeveJ“arcaSuted  ^JortJ”  '' 

and  the  beam-splitting  ratio  (K  =  2  ^±00Qf)iZZ  reading 

had  been  calibrated  by  the  St2e  liirol^glcli  InstltSS®'  "" 

79-3±8.4lt?"™*"^^^  tonditions  were:  room  temperature  25S1.1®  c,  humidity 

As  experimental  animals  we  used  thr*pp  v.*4. 

S,  I^d^ra^Sl^-  " 

sre“:?frtS2ero™  a*xil?- 

sides  of  the  back  2  cm  from  the  sninp  were  drawn  on  both 

center  of  each  square  through  a  diaphragm  5  L  irdlaSter 

!LS2aUon“S2S  S'nofLscJlSd'j;  "f  2" 

study  in  this  experiment  Is  Indlcate^bv  a  oo^lf”’®  threshold  under 

we  sought  the  irradiation  level  at  wh^nh  *'ee'^lt  for  erythema  (+); 

results  for  erythema  2i|  hours  after  exDosurP^*  would  be  50  percent  positive 
examined  under  a  an  ontlcai  m^fv».v^  posure.  Seven  exposure  groups  were 

exposure;  the  groups  that  were  close^to\h^^rf^*'^"^  Intervals  following 
examined  under  the  electron  microscope.  ®  threshold  were  also 

22l12lonf°tKh:r-|,rr2str:''s2‘^’J  -^^tologlcal 

group.  The  weighted  tegrIssu2m:tL"(:n  fpppf, 

Obtain  the  minimum  reddening  dose  for  whii-p  n-j  i  ^  computer)  was  used  to 
User  and  the  95*  ^n^ldLTr^^!  *  ^ 
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Fig.  1 .  Optical  circuit  of  experimental  apparatus 


Table  1 .  Instantaneous  erythema  following  exposure  of  pig  skin 
to  308-nm  excimer  laser  light 


Exposure 

group 

Avg.  exposure 

±std. deviation 1 

mJ/cm^ 

Number  of 
points 
exposed 

Instantaneous 

erythema 

Outcome 

Time  to  appear, 
seconds  (”) 
after  exposure 

Persistence, 
min  ('),  sec  (•') 

1 

22.28d:2.26 

32 

4^ 

— 

— 

— 

2 

29.44±1.80 

82 

— 

•  — 

— 

6 

89.76±S.85 

80 

83.88 

48"~2a7" 

4 

ss.sedbs.io 

81 

41.04 

7"^12" 

1'28"~2'26" 

.i. 

5 

69.14db4.61 

27 

70.87 

1'54"^5' 

+  - 

6 

91.27±1.06 

83.88 

45"-8'42'' 

+ 

7 

106.17±3.M 

^5 

100.00 

instantaneous 
to  15” 

80"--l(y 

+  +  - 

Table  2.  Delayed  erythema  following  exposure  of  pig  skin 
to  308-nm  excimer  laser  light 


Exposure 

group 

Avg.  exposure 
istd.  deviation, 
mJ/cm* 

No.  of  points 
exposed 

No.  of  erythema 
points  2^  hours' 
after  exposure 

Delayed  erythema 
rate,  24  hrs,  % 

Outcome 

1 

22.88db2.26 

82 

5 

15.68 

2 

29.44dbl.80 

82 

9 

2S.18 

± 

8 

89.76±d.86 

80 

11 

86.67 

4 

98.S6db9.19 

81 

14 

45.16 

5 

69.14±4.61 

27 

15 

.  56. W 

+  4--, 

6 

91.27±1.96 

80 

22 

78.88 

+  -f- 

7 

10G.17±8.18 

80 

25 

83.88 

+  +  + 
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2.  Results  and  Analysis 

Tables  1  and  2  present  data  on  instantaneous  and  delayed  reddening  of  skin 
after  irradiation  with  308-nm  laser  light.  The  delayed  erythema  generally 
appeared  8  to  9  hours  after  exposure,  but  in  groups  6  and  7  it  appeared  at 
some  irradiated  locations  5  after  exposure.  The  number  of  erythema  spots 
reached  a  peak  in  all  erythema  groups  24  hours  after  irradiation,  with  a 
rather  deep  color;  the  spots  generally  persisted  for  2  to  3  days.  The  spots 
in  the  sixth  group  had  already  become  pigmented  48  hours  after  exposure,  with 
the  pigmentation  persisting  for  6  to  14  days.  Scabs  formed  on  about  63 
percent  of  the  delayed  erythema  spots  in  group  7. 

The  results  of  standard  histological  examinations  of  the  skin  of  the  exposure 
groups  were  as  follows:  sections  of  the  delayed  erythema  points  24  hours 
after  irradiation  showed  no  evident  signs  of  reaction  in  the  first  group;  in 
groups  2-4  the  capillaries  in  the  dermis  were  dilated  and  congested  to  various 
degrees  and  small  numbers  of  erythrocytes  had  been  exuded.  In  the  fifth  group, 
in  addition  to  the  above  phenomena,  lymph  cells  appeared  in  the  vicinity  of 
the  blood  vessels.  In  group  6  the  dermis  contained  exterior  acervate 
Inflammation  cells;  in  group  7  the  basal  cells  showed  intracellular  edema 
accompanied  by  small  numbers  of  vacuoles,  and  the  cells  were  rather  loosely 
arranged.  More  detailed  information  is  presented  in  the  photographs  (a-d). 

In  groups  3»  4  and  5  the  instantaneous  erythema  was  studied  under  the  optical 
and  electron  microscopes.  The  basic  structure  was  normal.  Photographs  e  and 
f  are  from  the  electron  microscope  examination. 

The  MRD(._  values  and  95  percent  confidence  limits  for  irradiation  with  308- 
nm  laseP^light  based  on  visual  determination  and  weighted  regression  analysis 
are  presented  in  Table  3  and  Fig.  2. 

The  X  axis  is  the  logarithm  of  the  exposure,  while  the  y  axis  gives  the 
delayed  erythema  occurrence  probability.  A  chi-square  test  of  the 
regression  formula  shown  in  the  table  gave  satisfactory  results. 

3.  Discussion 

a.  Two-Phase  Erythema  Reaction  Produced  by  Short-Pulse  Ultraviolet  Laser 
Light 

Ultraviolet  Irradiation  causes  a  two-phase  (instantaneous  and  delayed)  skin 
erythema  reaction,  which  was  clearly  evident  in  our  experiment.  The 
instantaneous  erythema  produced  by  308-nm  laser  light  appeared  early  and 
cleared  up  rapidly,  persisting  no  longer  than  15  minutes.  This  phenomenon  may 
be  produced  by  the  pressure  effect  of  short-pulse  ultraviolet  light:  only 
when  the  pressure  reaches  some  threshold  value  is  blood  vessel  permeability 
altered,  producing  the  instantaneous  reaction.  The  delayed  erythema  appears 
only  after  a  certain  time  interval:  this  is  a  charatceristic  of  optically 
induced  chemical  reactions.  Some  foreign  investigators  have  reported  that 
histamine,  serotonin,  kinins  and  prostaglandins  (E,  F  and  G)  may  be  the 
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Photographs,  (a)  Histologic  examination  of  tissue  of  group  4,  2^ 
hours  after  exposure  (HE  stain).  Capillaries  of  dermis  dilated  and 
congested,  (b)  Histologic  examination  of  tissue  of  group  5,  24 
hours  after  exposure  (HE  stain).  Capillaries  of  dermis  dilated  and 
congested,  (c)  Histologic  examination  of  tissue  of  group  6,  24 
hours  after  exposure  (HE  stain).  Capillaries  of  dermis  dilated  and 
congested.  Edema  of  epidermal  cells,  (d)  Histologic  examination 
of  tissue  of  group  7,  24  hours  after  exposure  (HE  stain).  Capillaries 
of  dermis  dilated  and  congested.  Acervate  inflammation  cells.  Edema 
of  epidermal  cells,  (e)  Specimen  taken  5  minutes  after  level  IV 
exposure.  Gaps  between  epidermal  prickle  cells  widened  (electron 
microscope,  11,200X).  (f)  Specimen  taken  5  minutes  after  level  V 

exposure.  Gaps  between  epidermal  prickle  cells  widened  (electron 
microscope ,  11, 200X ) . 
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products  of  the  photochemical  reaction  that  cause  the  delayed  reaction.  To 
gain  a  further  understanding  of  the  mechanism  of  ultraviolet-produced 
erythema,  biochemical  studies  will  be  needed. 


Table  4.  Probability  analysis  of  skin  damage  to  white 
pigs  resulting  from  exposure  to  308-nm  laser  light 


Visual 

inspection 

Weighted  regression 
roethocf 

Regression  formula 

y-0.4759-f2.6184« 

MBDto 

(mJ/cm^) 

68.09 

68.84 

95%  confidence 
interval,  mJ/cm* 

48.86~65.00 

45.97-^68.05 

Probability 

units 


Fig.  2.  Probability  of  development  of  erythema  versus 
logarithm  of  exposure 


Solid  line  represents  results  of  visual  examination; 
dashed  line  is  weighted  regression  lines.  Open  circles 
represent  experimental  points,  dark  circles  represent 
values  calculated  from  weighted  regression  formula. 


b.  The  Minimum  Reddening  Dose  (MRD)  and  Laser  Protection  Standards 

This  experiment  used  white  pigs  to  determine  the  damage  threshold  for 
the  irradiation  of  skin  with  308-nm  pulsed  ultraviolet  laser  light.  MRDp..  was 
53.84  mJ/cm^,  i.e.  spots  receiving  this  dose  had  a  50%  probability  of 
erythema.  Similarly,  from  the  regression  equations  we  can  also  determine 
that  the  dose  at  which  there  will  be  a  10%  incidence  of  erythema  (MRD 
is  6.91  mJ/cm2.  in  existing  international  laser  protection  standards 
(lEC,  ANSI,  WHO  Environmental  Protection  Standard  23  and  the  like),  the 
maximum  permissible  irradiation  dose  for  308-nm  pulsed  laser  light  is  6.19 
mJ/cm^.  The  MRD^q  value  is  only  10%  higher  than  the  MPE  (maximum  permissible 

exposure),  indicating  that  the  boundary  between  nominal  safety  and  possible 
damage  is  unclear.  Consequently,  to  assure  safety  when  using  supershort-pulse 
ultraviolet  lasers,  reference  must  be  made  to  existing  maximum  permissible 
exposure  standards. 


References  omitted. 
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[Article  by  Xu  Zhongying  [1776  0112  5391],  Institute  of  Semiconductors, 

Chinese  Academy  of  Sciences:  "Picosecond  and  Femtosecond  Spectroscopy  in 
Semiconductors";  Paper  received  12  February  1985;  revised  manuscript  received 
29  May  1985] 

[Text]  Abstract:  This  paper  describes  various  time-resolved  nonlinear  optical 
techniques  which  are  applicable  to  semiconductors,  such  as  hot  luminescence 
correlation  technique,  transmission  correlation  technique,  population  mixing 
technique,  up-conversion  gate,  and  pump-probe  technique.  The  examples  of 
their  applications  to  semiconductors  are  given. 

I.  Introduction 

The  development  of  modern  picosecond  (ps)  technology  provided  an  extremely 
useful  tool  for  research  on  fast  nonequilibrium  phenomena.  These  fast  nonequi¬ 
librium  processes  are  generally  extremely  short,  from  several  hundred  pico¬ 
seconds  to  several  tens  of  femtoseconds  (fs).  This  time  domain  has  special 
significance  for  solid-state  physics.  Many  Important  physical  phenomena,  such 
as  thermoelectron  relaxation,  nonequilibrium  thermophonon  relaxation,  and  the 
electron-hole  collision  process,  occur  within  several  picoseconds  or  even 
shorter  times. 

When  semiconductor  material  is  excited  by  photons  with  energy  of  hVo,  a  large 
amount  of  nonequilibrium  hot  carrier^  is  produced.  These  hot  carriers  possess 
extra  energy  hVo  ~  which  the  electron  extra  energy  is  AEg  =  (hVo  -  Eg)  x 

(1  +  me/mh)"’-.  These  high  excitation  state  thermoelectrons  can  rapidly  lose 
their  extra  energy  and  relax  to  the  bottom  of  the  conduction  band.  Studying 
tHe  relaxation  process  of  these  thermoelectrons  can  provide  much  information 
on  the  electron-electron  and  electron-phonon  interaction  in  semiconductors. 

Generally  speaking,  under  the  action  of  an  extremely  short  light  pulse,  high 
excitation  state  nonequilibrium  carriers  first  of  all  go  through  carrier- 
carrier  scattering,  plasma  yuan  [0337]  emission  and  LO  phonon  scattering  to  a 
quasi-equilibrium  state.  At  this  time  the  electron  temperature  is  much  higher 
than  the  temperature  of  the  crystal  lattice,  the  relaxation  process  is  extreme¬ 
ly  short,  as  the  latest  experimental  measurements  have  shown,  in  GaAs, 
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Alo^32Gao,6BAs  and  GaAs  duoliangzij ing  [1122  6852  1311  7088]  structures  are 
35, *60,  and  50  fs,  respectively.  Under  the  action  of  a  picosecond  laser 
pulse,  experimentally,  there  is  no  way  to  observe  the  above -described  femto¬ 
second  relaxation  process.  The  thermoelectrons  produced  by  excitation  are 
already  in  a  quasi-equilibrium  state,  and  the  electron  temperature  is  much 
higher  than  the  crystal  lattice  temperature.  The  relaxation  process  after 
this  is  carried  out  primarily  through  emission  of  phonons,  and  the  relaxation 
time  is  from  several  picoseconds  to  several  tens  of  picoseconds  determined  by 
crystal  lattice  temperature,  excitation  intensity,  electron  energy, etc.  The 
above-described  relaxation  process  can  be  illustrated  by  Figure  1. 


Figure  1.  Schematic  Drawing  Showing  (a)  Various  Relaxation  Process 

Following  Ultrafast  Laser  Pulse;  (b)  Carrier  Distribution 
at  Different  Time  Delays 


The  aim  of  this  paper  is  to  introduce  the  application  of  picosecond  and  femto¬ 
second  laser  spectroscopy  in  research  on  semiconductor  physics,  and  to  give 
some  representative  recent  research  results. 


II.  Picosecond  and  Femtosecond  Lasers 

To  produce  an  ultrafast  laser  pulse  it  is  necessary  to  use  laser  mode-locked 
technology.  There  are  now  several  mode-locked  lasers  which  produce  ultrafast 
laser  pulses  of  several  hundred  ps  to  12  fs.  Table  1  sets  out  the  primary 
characteristics  of  various  primary  ultrafast  pulse  lasers. 

Synchronous  pump  dye  lasers  and  colliding  pulse  mode— locked  dye  lasers  can 
produce  continuous  pulse  sequence,  high  repetition  frequency  and  permit  use 
of  signal  averaging  technique  and  are  already  widely  used  in  measuring  fast 
relaxation  processes.  The  synchronous  pump  dye  laser  also  has  such  advantages 
as  a  large  tuning  range,  broad  wavelength  fugai  [1788  5556  covering?]  range, 
and  can  simultaneously  synchronously  pump  two  or  even  many  dye  lasers  to 
implement  synchronous  output  of  many  wavelengths .  When  an  unusually  high 
excitation  is  required  (such  as  when  studying  such  physical  processes  as 
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Table  1.  Characteristics  of  Various  Ultrafast  Laser  Systems 


Laser 

type 

ESHI* 

YAG 

Dye 

Dye 

1  Dye 

BH 

Argon  ion 

Excitation 

mode 

Flash- 

lamp 

Flash- 

lamp 

Flash- 

lamp 

Continuous 

laser 

Mode 

locked 

laser 

Continuous 

laser 

Discharge 

0.1/s 

lO/s 

10/s 

Continuous 

Continuous 

Continuous 

Mode 

locked 

mode 

Passive 

Passive 

Passive 

Passive 

Colliding 
pulse  mode 
locked 

Active 

Working 

wavelength 

(mm) 

1,060 

1,060 

mniiii 

580-620 

420-900 

620 

515,  488 

Pulsewidth 

(PS) 

6-8 

30 

2-5 

0.2-1. 5 

0.5-50 

EBBR 

150 

Peak  power 

10’ 

10’ 

10 

10^ 

10^ 

10’ 

10  ’ 

*12  fs  is  the  latest  result  of  Dr  J.M.  Halbout  of  IBM,  and  is  currently  the 
shortest  light  pulse.  The  system  uses  pulse  compression  technology.  55  fs  is 
a  research  result  reported  by  Cornell  University's  Professor  Tang  Zhongliang 
[3282  0112  5328]  who  has  not  yet  adopted  pulse  compression  technology. 


exciton  molecules  and  electron-hole  drop)  frequently  a  high  power  Neodymium 
glass  and  YAG  solid-state  laser  is  used. 

The  colliding  pulse  mode  locked  laser  was  proposed  in  1981  by  R.L.  Fork,  et  al., 
of  Bell  Laboratories.  This  is  a  ring  laser.  Within  the  chamber  there  are  two 
pulses  propagated  in  opposing  directions  precisely  synchronous  in  saturation 
absorber.  To  produce  a  transient  standing  wave  light  field,  a  transient  gate 
of  absorber  molecule  and  particle  distribution  thus  effectively  reduces  loss 
within  the  chamber  and  functions  as  a  sychronous,  stable,  compressed  laser 
pulse.  The  greatest  characteristic  of  this  laser  is  that  its  structure  is 
simple  (see  Figure  2)  consisting  of  a  continuously  operating  argon  ion  laser, 
gain  dye  jet  (R6G,  40-50  y®)  and  saturation  absorbing  dye  jet  (DODCI,  10  ym) . 
Its  output  pulse  width  is  100  fs,  but  using  pulse  compression  technique  the 
pulse  width  can  be  further  reduced. 

III.  Measuring  Technique  and  Its  Application  in  Semiconductor  Physics 

Up  to  now  there  has  not  been  a  picosecond-femtosecond  photoelectron  detector. 
The  only  electronic  device  which  can  record  picosecond  processes  is  the  streak 
camera,  but  this  device  is  high-priced,  has  low  resolution,  and  the  stability 
and  left  are  not  very  ideal. 
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ooDCi  absorber 


Figure  2.  Schematic  Representation  of  Colliding  Pulse  Mode-Locking 
(CPM)  Dye  Laser  System 

At  present  in  ultrafast  laser  pulse  measurement  technique,  nonlinear  optical 
measurement  technology,  common  correlation  measurement  technique,  frequency 
up-conversion  technique,  and  pump-probe  technique  are  mainly  used. 

1 .  Correlation  Measurement  Technique 

Correlation  measurement  technique  is  a  newly  developed  picosecond  spectral 
measurement  technique  which  is  very  similar  to  the  autocorrelation  measurement 
technique  used  t<?  measure  picosecond  laser  pulses  themselves .  Here  we  will 
introduce  photoluminescence  correlation  measurement  technique,  transmission 
correlation  measurement  technique^  and  population  mixture  technique,  and  pro¬ 
vide  the  experimental  results  of  their  use  in  measuring  thermoelectron  telaxa- 
tion  in  GaAs  and  GaAs  duoliangzijlng  and  heavy  doped  P-GaAs  shaozi  [1421  1311] 
life. 

(1)  The  photoluminescence  correlation  measurement  setup  is  as  illustrated  in 
Figure  3.  The  5  ps,  10®  Hz  laser  pulse  series  emitted  by  the  synchronous 
pulse  dye  laser  (DCM  dye,  650  nm)  is  divided  by  the  polarizing  prism  into  two 
beams  of  equal  intensity  and  perpendicular  polarized  direction  and  after  intro¬ 
ducing  a  time  delay  t,  the  two  beam  pulses  may  merge  through  a  focused  excita¬ 
tion  sample,  after  the  sample  is  excited  it  produces  a  large  volume  of  thermo¬ 
electrons,  and  the  thermoelectrons  quickly  relax  to  the  bottom  of  the  conduc¬ 
tive  band.  Although  the  relaxation  velocity  far  exceeds  the  radiative  recom¬ 
bination  velocity,  there  still  exists  a  transition  in  the  high  energy  band 
which  forms  a  high  energy  tail  in  the  photoluminescence  spectrum  (illustrated 
in  the  inset  in  Figure  3) .  Due  to  the  dual  molecular  characteristic  of  the 
recombination  luminescence,  the  luminescence  intensity  displays  a  nonlinear 
relationship  with  the  excitation  intensity,  thus  when  the  two  pulses  coincide 
(t  =  0),  the  luminescence  intensity  intensifies  the  nonlinearity  forming  a  hot 
luminescence  correlation  peak.  Any  change  in  the  thermoelectron  distribution 
is  reflected  in  the  luminescence  correlation  signal  intensity.  When  t  greatly 
exceeds  the  thermoelectron  relaxation  time,  the  correlation  peak  disappears, 
and  in  experiments  only  DC  background  is  seen.  The  relationship  between 
luminescence  correlation  peak  and  the  delety  time  T  between  the  two  pulses  can 
be  expressed 

PL(T)'^|_^diBexp(  -  \x\/Tt)AO('t-a>), 


(1) 


Polarizing  light 
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Figure  3.  Experimental  Setup  of  Hot  Luminescence  Correlation 
Technique 


in  which  x  is  any  variable,  Tg  is  the  thermoelectron  relaxation  time  constant, 
J  dtl{t+v)l(li)  is  the  correlation  function  of  the  laser  pulse  in  which 
I(t)  and  I(t  +  t)  is  the  intensity  of  the  two  pulses. 


Equation  (1)  was  derived  on  the  basis  of  the  changes  in  the  time  of  carrier 
occupation  density,  and  shows  that  the  luminescence  correlation  signal  is  the 
convolution  of  the  laser  pulse  autocorrelation  function  and  the  exponential 
function  of  time  constant  as  Tg.  When  the  laser  pulse  is  far  smaller  than  the 
thermoelectron  relaxation  time,  PL(t)  is  the  time  constant  equivalent  to  the 
exponential  function  of  electron  relaxation  time  Tg;  yet  when  the  electron 
relaxation  function  and  the  laser  pulse  width  are  analogous,  it  is  necessary 
to  solve  the  correlation  function  to  find  the  relaxation  time.  The  results  of 
relaxation  process  of  thetmoelectrons  experimentally  measured  in  GaAs  and 
GaAs-GaAlAs  duoliangzljing  showed  that: 

a)  Under  the  same  excitation  conditions,  the  cooling  velocity  of  thermoelec¬ 
trons  in  liangzijlng  [6852  1311  7088  literally  "quantum  trap"]  structures  is 
slow  and  the  relaxation  time  is  long.  This  is  because  in  the  cooling  process, 
the  thermoelectrons  emit  a  large  amount  of  LO  phonons,  and  the  multilayered 
boundaries  of  the  liangzijlng  material  causes  the  thermoelectron  system  to 
encounter  obstacles  in  crystal  lattice  energy  conversion.  Figure  4  is  a  typi¬ 
cal  hot  luminescence  correlation  curve.  The  black  dots  in  the  figure  repre¬ 
sent  the  computer  results,  from  the  figure  we  obtain  a  4.8  ps  relaxation  time 
for  the  GaAs  semiconductor  material  thermoelectrons  and  35  ps  for  the  llangzi- 
jing  structure. 

(b)  The  electron  relaxation  relocity  of  different  energy  positions  is  differ¬ 
ent,  thus  this  method  can  be  used  to  conveniently  study  the  relaxation  pro¬ 
cesses  of  electrons  in  different  energy  states.  Figure  5  illustrates  the 
relaxation  process  of  thermoelectrons  of  different  energies  in  a  GaAs 
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Figure  4.  Typical  Hot  Luminescence  Correlation  Peaks  of  QW  GaAs  and  Bulk 
GaAs 


Figure  5.  Comparison  of  Hot-Luminescence  Correlation  Peaks  at  Different 
Photon  Energies  for  GaAs  QW  Under  the  Excitation  Intensity  of 
-3.8  X  lO^Vcm® 

(a)  «i,-U7ineV,  T,-84ps;  (b)  Bj,-15Om0V,  T,-51p8; 

(c)  ®i.-190ineV,  T,-35p8;  (d)  »i.-22Om0V,  T,-23ps 
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llangzijlng,  the  black  ddts  in  the  illustration  represent  the  computer  results, 
the  theoretical  computation  of  polarization  of  the  experimental  curve  near  the 
low  energy  state  t  =  0  is  created  by  the  filling  effect  of  the  electrons 
going  from  high  energy  state  to  low  energy  state. 

(c)  The  thermoelectron  relaxation  time  is  related  to  the  excitation  level: 
under  high  excitation,  the  relaxation  process  is  slow.  This  is  because  the 
electron-phonon  system  screen  greatly  reduces  the  energy  exchange  between  the 
carrier  and  the  crystal  lattice.  Figure  6  summarizes  the  experimental  results. 
In  the  figure,  the  numbers  represent  the  energy  of  the  electrons  under  study 
and  it  can  be  seen  that  at  low  excitation  (•'^lO^’cm"^)  the  relaxation  time  of 
the  semiconductor  materials  and  the  llangzijlng  structure  tend  to  be  the  same. 

If  we  use  a  CPM  laser  with  a  pulse  width  of  0.1  ps  to  replace  the  synchronous 
pump  dye  laser  with  a  pulse  width  of  5  ps,  we  can  observe  the  correlation  peaks 
illustrated  in  Figure  7.  Based  on  the  above  principles,  we  can  obtain  a  relax¬ 
ation  time  of  the  electron  from  high  excitation  state  to  quasi-equilibrium 
state  of  approximately  ~100  fs.  In  addition,  we  also  observe  that  when  the 
excitation  intensity  increases,  wings  appear  on  the  two  sides  of  the  correlation 
peak,  the  time  constant  is  approximately  1  ps,  the  relative  intensity  increases 
as  the  excitation  intensity  Increases,  and  it  is  the  function  of  electron 
energy.  When  energy  is  low,  the  relative  intensity  of  the  wings  increases. 
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Figure  6.  Hot  Carrier  Relaxation  Time  vs.  Excitation  Intensity  for  Both 
QW  and  Bulk  GaAs  at  Different  Energies 


Figure  7.  Rising  Wings  Due  to  Band-Filling  Effect 

The  wings  are  formed  because  in  the  central  energy  valley  electron  relaxation, 
the  electrons  initially  scattered  to  the  energy  valley  jointly  occupied  by  L 
and  X  largely  leak  towards  T  energy  valley  and  produce  the  energy  band  filling 
effect.  As  it  draws  close  to  the  bottom  of  the  conduction  band,  the  energy 
band  filling  effect  becomes  clear  and  the  resultant  relative  intensity  of  the 
wings  increases. 
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Figure  8.  Experimental  Setup  for  Transmission  Correlation  Technique 
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(2)  Transmission  correlation  measurement  technique  is  based  on  the  absorption 
saturation  effect,^  see  Figure  8.  The  coaxial  focused  excitation  sample  after 
the  -100  fs,  80  MHz  laser  pulse  emitted  by  the  CPM  laser  has  undergone  polariz¬ 
ing  beam-splitting  and  time  delay,  the  DC  probe  system  probes  the  overall  trans¬ 
mission  intensity  of  the  two  pulses  through  samples.  To  improve  the  sensitiv¬ 
ity  of  the  probe  and  reduce  noise  the  differential  reception  technique  was 
adopted.  The  differential  signal  is  as  Illustrated  in  Figure  9(a).  Figure  9(b) 
is  the  transmission  correlation  peak  after  integration.  The  broken  line  is  the 
autocorrelation  signal  of  the  laser  pulse.  The  overall  transmission  intensity 
measured  with  changes  in  the  delay  time  x  can  similarly  be  expressed  in  the 
form  of  Equation  (1) . 


Figure  9.  Transmission  Correlation  Curve 

(a)  The  derivative  dither  trace  of  transmission  correla 
tion  peak; 

(b)  The  corresponding  integrated  curve  of  transmission 
correlation  peak 
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The  femosecond  relaxation  processes  measured  in  GaAs,  Gao ^ esAlo. 32AS  liangzi- 
j ing  structures  were  35,  60,  and  50  fs,  respectively.  Note  that  these  values 
are  much  smaller  than  the  width  of  the  100  fs  laser  pulse  itself.  The  energy 
of  the  photon  emitted  in  the  experiment  was  2.02  eV,  the  emission  intensity 
was  ~10^’cm“^.-  This  is  the  fastest  relaxation  time  measurement  at  present. 

Wings  on  both  sides  of  the  correlation  peaks  were  also  observed  in  transmis¬ 
sion  correlation  measurement.  Their  time  constant  was  1  ps. 

(3)  Population,  mixing  technique  is  another  correlation  measurement  technique 
used  for  measuring  the  carrier  dynamics  in  semiconductors.  As  illustrated  in 
Figure  10,  after  two  laser  pulses  of  equal  intensity  go  through  time  delay  t, 
an  excitation  semiconductor  sample  is  focused,  producing  photoluminescence. 

The  two  laser  beams  are  modulated  by  a  chopper  which  uses  different  frequencies, 
when  the  phase  locked  amplifier  is  adjusted  to  slip  frequency  (sum  frequency) 
frequency,  one  can  probe  the  transition  which  occurs  between  the  electrons  pro¬ 
duced  by  the  excitation  of  one  beam  and  the  hole  produced  by  the  excitation  of 
the  other  beam  and  the  electron-hole  combined  luminescence  signal  excited  by 
the  samp  laser  pulse  is  isolated  by  heterodyne  receiving  technique.  Supposing 
the  background  electron  and  hole  concentration  of  the  material  is  no  and  po» 
the  concentrations  of  the  carrier  produced  by  the  first  laser  pulse  are  n(t) 
and  p(t),  the  concentrations  of  the  carrier  produced  by  the  second  laser  pulse 
are  n(t  +  t)  and  p(t  +  x) ,  the  overall  combined  luminescence  intensity  is 


I(t,  r)-^[n(<+T)+«(0+«o] 


(2) 


Since  heterodyne  reception  was  adopted,  the  signal  intensity  obtained  was 


n(^t+T)p(t)dt+^  n(t)p(t-\-T)dt, 

1. e.,  the  luminescence  signal  measured  is  in  direct  proportion  to  the  correla¬ 
tion  function  of  the  carrier  density,  and  clearly  also  is  the  correlation 
function  of  the  delay  time  of  the  two  laser  pulses. 

Rosen  et  al.^  used  this  technique  to  measure  the  life  of  a  few  carriers  in 
highly  doped  P-GaAs(Zn:  6  x  10^®/cm^).  The  laser  usedwas  a  ring  type  dye 
laser  (<0.5  ps,  610  nm)  .  The  results  of  the  experiment  are  given  in  Figure 
11.  The  electron  life  obtained  was  39  ±  7  ps. 

2.  Frequency  Up -Conversion  Technique 

The  frequency  up— conversion  technique  is  also  called  the  optical  sum  frequency 
technique.  It  uses  ultrafast  pulse  as  a  time  gate  to  measure  the  time  change 
process  of  transient  processes.  Compared  with  the  optical  Kerr— cell,  the 
frequency  up— conversion  technique  is  suitable  for  use  in  medium  power ,  highly 
repetitive  frequency  continuous  output  pulse  lasers  (such  as  the  synchronous 
pump  dye  laser)  where  the  measurement  signal  may  be  very  weak,  the  wavelength 
partition  capability  is  strong,  and  the  time-resolved  precision  is  limited  by 
the  width  of  the  laser  pulse,  thus  it  is  especially  suited  to  time-resolved 
spectrum  measurement . 
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Figure  10.  Schematic  Diagram  of  the  Population  Mixing  Technique 


Figure  11.  Measured  Correlation  Function  for  P-GaAs  (6.4  x  10^®Zn  cm 
at  115  K  and  840  nm 


The  basic  experimental  apparatus  for  frequency  up-conversion  is  Illustrated  in 
Figure  12.  To  measure  a  transient  spectrum,  the  laser  pulse  is  divided  into 
two  beams:  one  is  the  excitation  sample  which  produces  the  transient  lumines¬ 
cence;  the  other  beam  is  reflected  back  by  the  cube  comer  reflector  and 
introduces  time  delay  T.  The  two  beams  operate  coaxially  in  nonlinear 
crystals  to  produce  a  frequency  which  is  a  sum  frequency  the  light  of 
(03  =  Ci)i  +  (02  which  i^en  expressed  as  a  wavelength  can  be  written 


Xa 


(4) 
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Figure  12.  Schematic  of  the  Up -Conversion  Gate  Technique 


in  which  Xi  is  the  laser  wavelength;  Xz  is  the  wavelength  of  the  luminescence 
measured;  X3  is  the  wavelength  actually  measured  and  we  can  artificially  set 
up  a  monochromator  wavelength  (X3)  to  achieve  the  aim  of  measuring  different 
X2.  Clearly  this  can  be  used  to  obtain  a  group  of  spectrum  curves  of  differ¬ 
ent  delay  times.  On  the  other  hand,  if  a  monochromator  is  not  used  in  the 
above  apparatus  to  split  the  light  but  to  measure  the  overall  intensity  of 
the  luminescence,  the  process  of  changes  in  light  intensity  can  be  obtained 
for  different  time  delays  (luminescence  decay  process) . 

Daly  and  Mahr  used  this  technique  to  measure  the  time-resolved  spectrum  of 
CdSe’  using  a  synchronous  dye  laser  (10  ps)  with  excitation  intensity  of 
~5  X  10  cm“^.  Figure  13  is  a  graph  of  the  time-resolved  spectrum.  When 

T  =  20  ps  the  experiment  observed  that  the  central  wavelength  was  6280  A  wide 
luminescence  peak,  and  as  the  delay  time  increased  the  luminescence  peak 
became  narrower  and  shifted  towards  the  wavelength.  After  90  ps,  the  position 
of  the  peak  width  and  the  peak  value  were  basically  unchanged,  because  in  the 
process  of  cooling  the  plasma,  the  electron-hole  formed  liquid  drops.  Through 
this  kind  of  time-resolved  spectrum  we  could  also  obtain  important  information 
on  the  thermoelectron  relaxation  process.  For  example,  we  could  obtain  tran¬ 
sient  electron  temperatures  directly  on  the  basis  of  the  changes  in  intensity 
of  the  spectrum’s  high  energy  terminal  and  thus  obtain  the  time  evolution 
process  of  thermoelectron  distribution. 

3.  Extitation-Probe  Technique 

Under  a  great  many  conditions,  transient  relaxation  process  is  not  accompanied 
by  light  emission  (or  there  are  difficulties  in  measuring  light  emission)  and 
often  accompanying  the  changes  in  some  optical  properties  of  the  sample,  such 
as  changes  in  reflectivity  and  transmissivity,  under  these  circumstances  the 
excitation-probe  technique  is  often  used  for  measurements.  In  the  excitation- 
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Figure  13.  Luminescent  Intensity  Spectra  From  CdSe  at  1.8  K  for 
Fixed  Time  Delays  After  the  Excitation 


Figure  14.  Experimental  Arrangement  of  the  Time  Resolved  Nearband  Gap 
Absorption  Spectroscopy 
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probe  technique  the  first  beam  pulse  puts  the  sample  being  studied  into  the 
excitation  state,  then  the  second  pulse  (time  delay  t)  probes  the  time  suc¬ 
cession  process  of  the  excitation  state. 

C.V.  Shank,  et  al.,  used  the  excitation-probe  technique  to  measure  the  pro¬ 
cess  of  absorption  dynamics  of  the  light  excitation  carrier  in  GaAs,®  and 
observed  that  the  time  constant  of  thermoelectron  relaxation  to  the  crystal 
lattice  temperature  was  4  ps,  the  exciton  shield  and  band  gap  renormalization 
occurred  within  0.5  ps.  The  experimental  apparatus  is  illustrated  in 
Figure  14.  The  light  pulse  produced  by  the  mode  locked  dye  laser  goes 
through  amplification,  the  synchronous  alcohol  pool  and  water  pool,  the  former 
producing  jilaman  [3423  0812  2581  Gilman?]  scattering  producing  a  pulse  light 
excitation  sample  of  750oX,  and  the  latter  producing  a  white  light  pulse 
through  interaction  of  the  parameter's  four  photons  to  serve  as  a  probe  light 
with  a  sample  of  double  heterojunction  structure.  From  the  results  of  analy¬ 
tical  experiments  we  know  that  before  excitation  there  is  a  definite  exciton 
absorption  line;  when  x  =  1  ps  the  exciton  peak  disappears.  This  is  because 
the  great  number  of  excitation  carries  creates  an  exciton  shield;  we  observed 
that  the  absorption  boundary  shifted  to  undisturbed  free  excitons,  this  was 
caused  by  the  carrier  exchange  interaction  energy  band  compression;  when 
T  =  12  ps,  energy  greater  than  the  transmission  of  exciton  energy  increases, 
this  is  caused  by  the  energy  band  filling  effect  (see  Figure  15). 


Figure  15. 


.7900  8100  8300  7900  8100  ^300  7900o  8100  8309 

Wavelength  (A) 


Time-Resolved  Transmission  Spectra  and  Difference  Spectra 
of  GaAs  Near  the  Exciton  Resonance 


D.  Von  der  Llndle,  et  al.,  used  the  excitation-probe  technique  to  study  the 
time-resolved  phonon  spectrum  in  GaAs  and  InP  materials  (77K),^  but  the 
excitation  light  pulse  and  probe  light  pulse  used  came  from  the  same  light 
source  (synchronous  pump  dye  laser)  and  the  excitation  and  the  probe  wave¬ 
lengths  were  the  same.  The  experimental  apparatus  is  illustrated  in  Figure  16. 

GaAs  and  InP  belong  to  the  T^  point  group  crystals.  The  polarization  direc¬ 
tion  of  the  LO  phonon  scattering  light  which  is  scattered  from  (100)  faces  is 
perpendicular  to  the  polarization  direction  of  Incidence  light,  therefore 
using  a  suitable  polarization  plate  the  scattering  introduced  by  the 
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Figure  16.  Schematic  of  Time-Resolved  Phonon  Spectroscopy 


excitation  light  can  be  filtered  out.  The  excitation  light  Is  only  an  excita¬ 
tion  sample.  It  produces  LO  phonons  and  the  probe  light  at  the  same  time 
functions  as  excitation  and  probe,  but  the  LO  scattering  light  Introduced  by 
the  probe  light  Is  unrelated  to  At, but  Is  only  structural  background  and  the 
changes  In  LO  phonons  Introduced  by  excitation  light  creates  temporal  changes 
In  the  scattering  light  Intensity,  thus  Information  about  LO  phonon  relaxation 
Is  obtained.  In  the  experiment  we  obtained  quasl-equlllbrlum  LO  phonon  relax¬ 
ation  times  In  GaAs  and  InP  of  7  +  1  and  5  ±  1  ps,  respectively. 

IV.  Conclusion 

The  application  of  ultrafast  pulse  lasers  In  semiconductors  and  other  research 
areas  Is  very  lively,  for  example  picosecond  optical  waveguide  switches, 
physics  research  on  amorphous  materials, research  on  semiconductor  surface 
and  boundary  transient  state  processes, etc.  While  Introducing  the  dynamic 
state  of  research  In  this  area,  the  key  emphasis  In  this  paper  Is  on  the 
experimental  techniques  and  research  content  which  the  author  has  touched 
upon  In  past  work. 
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ULTRAFAST  IMAGE -CONVERTER  CAMERAS:  POWERFUL  TOOLS  FOR  TRANSIENT  RECORDING  OF 
OPTICAL  PHENOMENA 

Beijing  WULI  [PHYSICS]  in  Chinese  Vol  14  No  3,  Mar  86  pp  168-172,  164 

[Article  by  Zhou  Xuan  [0719  2467],  Institute  of  Semiconductors,  Chinese  Acad¬ 
emy  of  Sciences:  "A  Powerful  Tool  for  Capturing  Transient  Optical  Phenomena; 
Ultrafast  Image-Converter  Cameras"] 

[Text]  For  some  time  many  people  have  been  familiar  with  the  theory,  perform¬ 
ance  and  operation  of  ordinary  camera,  but  still  may  be  rather  unfamiliar  with 
ultrafast  photographic  technology.  As  the  term  suggests,  ultrafast  photo¬ 
graphic  technology  is  primarily  for  studying  ultrafast  moving  targets  or  tran¬ 
sient  optical  phenomena  and  complex  (photographic)  problems  of  elapsing  pro¬ 
cesses.  As  photographic  system,  "ultrafast"  refers  to  two  aspects:  one  means 
that  there  is  a  very  high  shooting  frequency  (described  as  the  number  of  Images 
that  are  shot  within  a  unit  of  time) ;  the  other  meaning  is  the  time  resolution 
ability,  i.e.,  how  fast  a  transient  process  the  camera  can  capture  (for  a 
shutter  camera,  described  as  the  shortest  exposure  time).  When  higher  shooting 
frequencies  are  required,  it  actually  also  demands  higher  time  resolution. 

For  ordinary  cameras  with  mechanical  shutters,  the  shortest  exposure  time  is 
several  thousandths  of  a  second,  and  this  is  sufficient  time  resolution  for 
shooting  most  moving  targets  in  dally  life.  But  the  many  optical  phenomena 
and  elapsed  processes  involved  in  scientific  research  are  much  shorter,  so 
that  using  an  ordinary  camera  to  shoot  such  high-speed  phenomena  is  completely 
impossible.  In  addition  to  demanding  extremely  short  exposure  times,  another 
difficulty  encountered  by  high-speed  photography  is  the  problem  of  the  light 
intensity  of  the  target  object.  The  essence  of  the  photographic  process  is  to 
cause  the  recording  photographic  plate  to  receive  a  certain  exposure  to  light. 
When  the  exposure  time  is  extremely  short  (which  is  necessary  in  shooting  high¬ 
speed  phenomena) ,  to  obtain  a  certain  light  exposure  requires  extremely  high 
intensity  of  illumination.  To  adapt  to  the  development  of  modern  science  and 
technology  various  kinds  of  high-speed  photographic  equipment  has  appeared, 
such  as  focusing  cameras,  Kerr  cell  (or  Pockel's  cell)  camera  and  holocameras 
which  use  laser  light  sources.  These  cameras  are  all  applied  to  different 
degrees.  However,  until  the  appearance  of  high-speed  cameras  using  modern 
Imaging  tubes  which  have  microchannel  plates  as  information  conversion  com¬ 
ponents,  no  one  camera  (such  as  the  image— converter  camera)  could  resolve 
sufficiently  the  contradiction  between  high-speed  response  and  high 
sensitivity. 
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Here  we  describe  two  ultrafast  image-converter  cameras— a  one-time  shutter  type 
and  a  scanning  streak  camera — ^vdiich  are  the  main  trends  in  current  ultrafast 
camera  development.  The  time  resolution  of  the  scanning  streak  camera  in 
particular  is  smaller  than  1  ps,  and  it  is  almost  the  only  Instrument  that  can 
be  used  for  studying  high-speed  phenomena  within  the  picosecond  range.  It  can 
be  anticipated  that  high-speed  image-converter  cameras  will  be  more  widely 
used  in  scientific  research  as  a  new  optical  measurement  instrument  the  more 
they  develop. 

I.  Basic  Problems  of  High-Speed  Photography  Capturing  One-Time  Optical 
Phenomena 

We  will  explain  some  of  the  basic  problems  which  must  be  resolved  in  high¬ 
speed  photographic  technology  using  photoelectric  study  of  gas  electric  dis¬ 
charge  as  an  example. 

Using  a  spark  gap  to  study  the  gas  electric  discharge  process  is  one  of  the 
basic  methods  of  understanding  the  characteristics  of  gas  breakdown  under 
high  pressure.  Gas  electric  discharge  has  the  following  characteristics: 

(1)  Gas  breakdown  is  determined  by  many  factors  such  as  the  voltage  between 
spark  gap  electrodes,  degree  of  gas  lihua  [4418  0553]  and  thus  it  is  diffi¬ 
cult  to  determine  accurately  the  time  at  which  breakdown  will  occur  before 
the  fact.  The  gas  breakdown  process  can  be  divided  into  three  stages.  In 
the  initial  stage,  even  though  a  DC  voltage  is  applied  to  the  spark  gap 
electrodes,  there  is  still  a  delay  before  it  produces  the  initial  electrons 
which  are  used  for  avalanche  multiplication.  When  the  avalanche  multiplica¬ 
tion  begins,  the  avalanche  current  rapidly  rises  to  its  critical  value.  Then 
the  gas  in  the  spark  gap  is  completely  broken  down  and  forms  a  conduction 
channel.  There  is  no  doubt  that  if  we  can  record  the  entire  breakdown  elapse 
process  we  can  obtain  a  great  deal  of  useful  Information  from  it.  But,  before 
gas  electric  discharge  reaches  its  critical  point,  it  is  only  accompanied  by 
production  of  a  very  weak  light  radiation,  and  thus  shooting  the  gas  pre¬ 
discharge  process  required  a  camera  with  extremely  high  sensitivity. 

(2)  The  duration  of  the  gas  discharge  process  is  very  short,  especially  in 

a  compressed  gas  with  a  very  high  dielectric  constant  (such  as  SFe)  the  dura¬ 
tion  has  already  entered  the  range  of  nanoseconds.  Thus,  to  capture  this 
transient  process,  a  camera  must  have  a  correspondingly  short  possible  expo¬ 
sure  time.  Gas  breakdown  is  actually  a  single  phenomenon  and  the  time  when 
the  breakdown  occurs  is  also  hard  to  predict  in  advance.  Thus  to  record  such 
a  process,  the  camera's  shutter  must  be  activated  by  the  measured  event  itself 
For  this  reason,  there  must  be  a  delay  time  from  the  time  the  event  occurs 
until  the  camera  shutter  opens.  And  this  time  must  be  as  short  as  possible 
so  that  it  can  record  the  initial  process  of  the  measured  phenomenon.  Of 
course,  by  lengthening  the  optical  path  between  the  measured  phenomenon  and 
the  camera  lens  the  influence  of  the  above  mentioned  delay  time  can  be  aug¬ 
mented  (in  air,  1  ns  corresponds  to  30  cm  of  optical  path).  However,  techno¬ 
logically  this  is  not  always  feasible,  because  it  makes  the  already  weak  light 
j^udiation  even  weaker .  The  best  method  of  solution  is  for  the  camera  itself 
to  have  the  ability  to  intensify  the  light. 
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(3)  Along  with  the  light  radiation  of  gas  discharge  there  may  be  a  very 
broad  optical  spectrum  range.  Thus,  the  camera  also  should  have  a  corre¬ 
sponding  light  spectrum  response  range,  and  it  is  best  if  it  can  convert 
incoming  light  of  a  very  broad  range  directly  into  the  visible  light  spec¬ 
trum  suited  to  observation  by  the  human  eye  or  recording  ordinary  film. 

(4)  The  aim  of  photography  is  to  reproduce  accurately  the  measured  target. 

To  study  the  spark  gap  discharge  channel  in  detail,  the  camera  also  should 
have  a  high  enough  spatial  resolution  and  small  enough  distortion. 

On  the  basis  of  the  above  characteristics,  the  camera  system  for  researching 
the  gas  discharge  process  can  be  reduced  to  the  following  basic  demands: 

1)  have  a  very  high  time  resolution;  2)  have  a  very  high,  adjustable  optical 
gain;  3)  can  be  triggered  by  the  measured  event  itself  and  have  the  smallest 
possible  trigger  delay  and  wobble;  4)  have  a  broad  optical  spectrum  response 
range;  5)  have  a  hl^ spatial  resolution;  6)  have  abundant  information  pro¬ 
cessing  functions;  and  7)  have  a  structure  which  is  compact,  small  and 
light. 

The  above  basic  demands  are  important  issues  which  high-speed  photographic 
technology  has  been  trying  to  resolve  for  a  long  time.  The  image-converter 
camera  can  take  into  account  the  above  basic  demands  and  has  already  become 
a  powerful  tool  for  capturing,  recording,  and  analyzing  transient  optical 
phenomena,  and  is  playing  a  very  Important  role  in  research  on  high-speed 
phenomena . 

II.  Comparison  of  Several  Common  High-Speed  Optical  Shutters 

Light  shutters  (optical  shutters)  are  the  central  part  of  any  type  of  shutter 
camera.  The  shortest  transient  optical  time  provided  by  mechanical  shutters 
is  generally  about  50  ns.^  The  Kerr  cell  and  Pockel's  cell  were  the  earliest 
electrooptical  components  to  use  an  optical  shutter,  but  they  both  have  the 
shortcomings  of  narrow  optical  response  range,  high  opening  voltage,  and 
considerable  loss  in  transmitting  light  energy.  The  appearance  of  the  new 
image-converter  based  on  the  principles  of  electronoptlcs  added  new  driving 
force  to  the  development  of  high-speed  photographic  technology. 

The  types  and  structure  of  image-converters  are  very  varied.  However,  all 
types  of  image-converters  encompass  one  basic  principle:  First  of  all,  an 
electron  analog  of  the  target  image  to  be  measured  is  formed  with  the  help 
of  a  photocathode ,  then  the  photoelectrons  released  by  the  photocathode  are 
controlled  so  that  they  are  accelerated  by  the  action  of  a  strong  electron 
field.  Finally,  these  high-energy  electrons  collide  with  a  fluorescent 
screen  (anode)  exciting  the  fluorescence  and  forming  a  secondary  optical 
image.  This  image  can  be  observed  by  the  human  eye  or  recorded  using  con¬ 
ventional  photographic  technology. 

Since  it  is  relatively  easy  to  control  and  amplify  the  electron  beam,  the 
image  converter  can  be  used  to  approximate  a  high-speed  inertia-less  optical 
shutter.  Image-converter  t3rpes  are  differentiated  mainly  by  the  different 
method  adopted  for  controlling  the  electron  beam.  One  worth  mentioning  is 
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a  recently  developed  method  called  the  close-frame  focusing  method,  which  in 
essence  makes  the  distance  between  the  photocathode  and  the  optical  screen 
electrode  surface  extremely  close  (only  several  mm)  so  that  there  is  formed 
iu  the  vicinity  of  the  photocathode  an  electron  cloud  which  can  be  accele¬ 
rated  by  a  strong  electrical  field  and  travel  directly  toward  the  optical 
screen  electrode  along  the  electrical  field  gradient  and  does  not  cause 
obvious  dispersion.  This  image-converter,  which  contains  only  a  photocathode 
and  an  optical  screen  electrode,  is  called  a  diode  image-converter.  When  it 
is  used  as  an  optical  shutter,  a  5-10  kV  pulse  voltage  is  applied  between  the 
two  electrodes.  If  a  third  electrode  (control  grid)  is  added  in  this  image- 
converter.',  then  it  forms  a  triode  image— converter .  The  control  grid  requires 
only  several  dozen  volts  to  be  able  to  control  the  electron  beam  effectively. 
In  the  past  dozen  years  a  new  close-frame  focus  image-converter  which  has  a 
microchannel  board  has  also  appeared.  Since  it  has  an  extremely  fast  reac¬ 
tion  time,  a  very  high  optical  gain  and  flexible  control  ability  it  has  been 
applied  extensively  in  high-speed  photography  technology.  See  Table  1  for  a 
comparison' '.of  the  typical  performance  of  the  different  electrooptical  high¬ 
speed  optical  shutters. 


Table  1.  Comparison  of  Typical  Performance  of  Different  Optical  Shutters 


Item 

Optical 

gain 

Response 

time 

Control 

voltage 

Mode  of 
application 

Kerr  cell 

0.1 

100  ps 

15-30  kV 

single  frame 

Pockel  cell 

0.5 

100  ps 

1-5  kV 

single  frame 

Diode  image-converter  tube 

10 

<1  ns 

5-10  kV 

single  frame 

Triode  image-converter  tube 

100 

<1  ns 

200  V 

multiframe 
or  streak 
recording 

Clo se-f ocus  image-converter 
with  MCP 

10 '*-10® 

<1  ns 

200  V 

single  frame 

III.  The  Simplest  Diode  Image-Converter 

The  basic  structure  of  the  diode  image— converter  is  illustrated  in  Figure  1. 
The  incident  light  shines  on  the  photocathode  P  and  through  it  releases  photo- 
electrons  which  form  an  electron  cloud  in  the  vicinity  of  the  photocathode 
surface.  If  a  DC  high  voltage  is  applied  between  the  photocathode  P  and  the 
anode  A,  then  the  photoelectrons  will  be  accelerated  by  the  strong  electron 
field  and  through  the  appropriate  electron  optical  system  0  also  be  focused  on 
the  optical  screen  electrode  (anode)  A,  and  thus  release  electrons  a  second 
time,  with  the  result  being  that  an  optical  image  is  formed  on  the  optical 
screen  electrode.  From  the  above-described  optical-electron  and  electron- 
optical  conversions,  the  incident  optical  image  is  enhanced.^  This  image- 
converter  which  has  image— enhancing  capability  is  called  an  image  enhancer. 
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It  should  be  pointed  out  that  there  is  no  added  electron  optical  system  0  fisr 
the  close-frame  focus  image-converter  with  an  extremely  close  distance  between 
the  photocathode  and  the  anode.  It  is  only  because  the  distance  between  the 
photocathode  and  the  anode  is  extremely  close  (several  millimeters) ,  when  a 
voltage  of  several  kV  is  applied  between  the  electrodes*  the  field  intensifi¬ 
cation  is  extremely  high.  Thus  the  photoelectrons  released  by  the  photo¬ 
cathode  will  move  in  a  uniform  electron  field  which  is  similar  to  a  panel 
capacitor*  and  the  accelerated  electron  field  reaches  the  optical  screen 
electrode  almost  parallel  and  thus  the  electron  beam  is  focused. 


Figure  1.  Diode  Image-Converter 


The  principles  of  one-time  high-speed  photography  using  the  close-frame  focus 
diode  image-converter  as  an  optical  shutter  are  illustrated  in  Figure  2. 

The  light  emitted  by  the  measured  object  G  is  focused  through  the  objective 
lens  0  on  the  surface  of  the  photocathode  P.  On  the  other  hand*  the  optical 
signal  synchronous  with  the  measured  optical  phenomenon  is  changed  through 
photoelectric  conversion  into  an  electron  trigger  signal  which  initiates  the 
high-voltage  narrow-pulse  exciter.  The  high-voltage  narrow-pulse  output  by 
the  exciter  is  applied  between  the  photocathode  and  optical  screen  electrode 
of  the  image-converter.  Within  the  pulse  width  action  time  (i.e.*  the  expo¬ 
sure  time)  the  photoelectrons  emitted  by  P  are  accelerated  by  the  strong  elec¬ 
tron  field  and  then  new  electrons  are  excited  by  collision  with  the  optical 
screen  electrode  forming  a  new  visible  image.  Through  lens  O2*  the  optical 
image  formed  on  the  optical  screen  electrode  is  recorded  as  an  image  on  film  F . 

Each  incident  high-energy  electron  which  collides  with  the  optical  screen 
electrode  surface  can  produce  many  photons  (typical  value  is  100-200)  and  thus 
the  image-converter  itself  is  a  component  with  a  very  high  optical  gain.  How¬ 
ever*  this  light  amplification  which  comes  from  electron  excitation  is  almost 
completely  canceled  out  by  the  two  following  factors .  First  *  the  typical 
quantum  efficiency  of  the  photocathode  is  only  about  20  percent;  second*  the 
photon  produced  by  collision  of  high-energy  electrons  with  the  optical  screen 
electrode  scatters  in  all  directions  completely  Irregularly*  thus  the  inten¬ 
sity  distribution  of  the  output  light  created  is  also  dispersed.  After  taking 
these  factors  into  consideration*  the  overall  light  gain  of  the  photographic 
system  illustrated  in  Figure  2  is  close  to  1*  thus  performance  cannot  be 
satisfactory*  and  it  is  especially  not  suited  to  weak  light  targets. 
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Figure  2.  Principles  of  Single  Event  High-Speed  Photography  Using 
an  Image-Converter  as  Optical  Shutter 


IV.  New  Image-Converter  Tubes  With  MicroChannel  Plate  and  Optical  Fiber 
Coupling  Window 

Two  major  technological  advances  which  have  been  made  in  recent  years  have 
greatly  improved  the  performance  of  the  diode  image-converter  tube  described 
above.  The  first  advance  benefited  from  achievements  in  optical  fiber  tech¬ 
nology.  Using  an  optical  fiber  window  to  replace  the  ordinary  output  rear 
window  of  the  diode  image-converter  tube  permitted  transmitting  the  optical 
image  produced  on  one  plane  to  another  plane  almost  without  distortion.  Since 
the  optical  image  transmission  efficiency  was  improved,  the  overall  optical 
gain  of  the  photographic  system  improved  from  nearly  1  to  15-30.  The  second 
advance  in  technology  was  the  adoption  of  the  so-called  mlcrochannel  plate 
(MCP)  structure . ^ 


Optical  screen  electrode  i 


1  kV  10— 20kV 


Figure  3.  Shutter  Tube  With  MCP  and  Optical  Fiber  Output  Window 

The  structure  of  a  new  close  focus  shutter  tube  with  an  MCP  and  optical  fiber 
output  window  is  illustrated  in  Figure  3.  The  distance  between  the  photb- 
cathode  and  the  MCP  input  plane  and  between  the  optical  screen  electrode  and 
the  MCP  output  plane  are  several  millimeters  (i.e.,  close  frame  focus).  The 
photocathode  potential  should  be  lower  than  the  MCP  input  terminal  potential, 
the  photoelectrons  which  come  from  the  photocathode  and  carry  the  image  infor¬ 
mation  project  onto  the  input  terminal  surface  of  the  MCP,  and  enter  the  micro- 
channel.  An  acceleration  voltage  typically  about  1  kV  is  applied  between  the 
MCP's  input  and  output  terminal  surfaces.  After  the  electrons  entering  the 
mlcrochannel  have  been  multiplied  they  appear  again  at  the  MCP's  output 


65 


io“of  Tf\Zr:i\TZ  a  high 

optical  screen  electrode  where  they  excite  fluorescencL^^Jie°f?^^^^ 

graphlc\ysSrt"iifcoLtrLSrobL"Ss\foverIiro\^^ 

hundred  thousand-fold  between  the  nh^fnf.^.^^,r.^o  i  ^  ^  ^  several 

It  must  be  pointed  out  photographic  film! 

f^S?convertefiuh°MS”anfopUca^fT®* 

to  tooting  hlgh-opoL  SoaroptLl  JatgotL"'''’ 


V. 


streak  Image-Converter  Tube  With  MicroChannel  Plate 


Tti  srr  ^  « 

main  difference  between  the  etreat  ^  scanning  streak  camera  system. ^  The 
trated  In  pSure  3  ifthaj  tf rf  “"■*  “'’“tter  tube  lllus- 

plates  and  tL  necesLJJ^ddJtSn^f^nrtr'd^f 'liJf  daf 

c'^ttpif^L^e'fp'i-n:  s:i-d€rF' 

Jr^ranf^u^ll'LS^P^lLf^y^tS  and 

rXt“=drrni 

Ing  the  deflecting  plates  is  develoLd  the  photoelectric  beam  enter- 

the  deflecting  platrSaL  It  t  I  t  ^  direction  perpendicular  to 

ss£rL=:5-s^^ 

s«eak'’“J™;ySrbrt8SneJr^„rcS:g:ri‘rJbir:“e:iEn^ 

oS^sLirby^Lr 

r%TS!  ?baJ-lt“»  “ariSt.:^=°^“““"  tbar/J:  “Jd^f la 


66 


Input 

laser 

pulse 


Optical 

delay 


High- 
speed 
photo- J 
electri 


P/C  MCA  S/P  P 
Slit 


IT 


Streak  imaged  mcp 
converter  ' 


detector 

Adjustable-] 


Igh-voltage  scan 
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Figure  4.  Simple  Scanning  Streak  Camera  System 


L  is  lens;  P/C  is  photocathode;  M  is  aperture  plate;  C  is  focusing 
electrode;  A  is  anode;  S/P  is  deflecting  plate;  MCP  is  micro- 
channel  plate;  P  is  optical  screen  electrode 


VI.  Image-Converter  Tube  Pulse  Control 

The  shutter-type  image-converter  tube  and  the  gating  pulse  generator  are 
combined  to  form  the  simplest  shutter-type  image-converter  camera  with  high 
optical  gain  and  rapid  response.  Similarly >  the  scanning  streak  image- 
converter  tube  and  the  scan  generator  which  controls  the  photoelectron  beam 
deflection  are  combined.  There  is  no  doubt  whatsoever  that  the  performance 
of  these  two  image-converter  cameras  is  primarily  determined  by  the  image- 
converter  used  itself,  but  here  there  is  also  the  problem  of  how  to  use  the 
image-converter  performance  effectively.  In  particular,  the  two  important 
performance  norms  of  the  shutter-type  image-converter  camera — optical  gain 
and  exposure  time — are  directly  determined  by  the  gating  pulse  amplitude  and 
width;  and  the  most  important  performance  norms  of  the  scanning  streak 
camera — streak  speed — is  directly  determined  by  the  scan  voltage  pulse  ampli¬ 
tude  and  ascent  time.  In  sum,  at  bottom  the  two  are  summed  up  in  the  demands 
of  high  voltage-high-speed  pulse,  and  after  deciding  to  use  the  image- 
converter,  image-converter  pulse  control  becomes  one  of  the  key  technical 
questions  in  image-converter  research  and  development. 

The  basic  demands  for  an  image-converter  controlled  pulse  generator  can  be 
summarized  as  follows:  1)  have  a  very  high  output  pulse  amplitude  and  the 
polarity  demanded;  2)  have  the  pulse  shape  demanded  (for  shutter-type  cameras 
it  must  be  a  slope  wave  of  excellent  linearity)  and  narrow  pulse  width  (or 
rapid  ascent) ;  3)  the  generator  is  triggerable  (electronic  or  optical)  and 

have  the  smallest  possible  trigger  delay  and  trigger  wobble;  4)  have  strong 
tolerance  negative  load  ability;  and  5)  have  small  mass,  light  weight,  and 
long-range  reliable  application. 

For  a  fairly  long  time  athyratron  or  a  spark  gap  was  used  as  the  control  pulse 
source  of  the  image-converter.  The  trigger  delay  and  wobble  Involved  with 
this  are  comparatively  large.  With  the  appearance  of  several  new  types  of 
high-speed  semiconductor  components,  the  image-converter’s  control  pulse 
source  has  gradually  been  replaced  by  a  solid-state  pulse  source  with  excel¬ 
lent  performance.^”^ 
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VII.  Basic  Organization  and  Application  of  Ultrafast  Image-Converter  Cameras 

As  we  said  above,  there  are  a  variety  of  ultrafast  photographic  systems.  The 
image-converter  and  the  corresponding  control  pulse  trigger  are  two  indispens¬ 
able  basic  constituent  parts  of  a  high-speed  camera  which  uses  an  image- 
converter  as  the  information  conversion  component.  Figure  5  illustrates  the 
basic  block  diagram  of  an  ultrafast  image-converter  camera  system.  The  design 
of  the  image  information  processing  element  in  the  system  is  very  flexible  and 
can  be  determined  depending  on  the  specific  situation. 


Figure  5.  Composition  of  Ultrafast  Image-Converter  Camera 


The  high-speed  camera  Introduced  so  far  is  called  the  "single  frame"  camera, 
i.e.,  corresponding  to  one  trigger  signal  (synchronous  with  the  measured  opti¬ 
cal  phenomenon)  there  is  produced  a  gating  narrow  pulse  to  be  used  by  the 
shutter-type  image-converter  or  a  high-voltage  fast  slope  pulse  to  be  used  for 
a  streak  image-converter  scan,  causing  the  high-speed  camera  within  the  corre¬ 
sponding  time  to  photograph  an  image  of  the  target  and  obtain  one  photographic 
frame.  Although  the  exposure  time  (or  scan  time)  is  very  short,  the  time 
interval  between  the  two  neighboring  exposures  (or  two  scans)  is  very  long, 
thus  in  reality  it  can  be  accepted  that  the  camera  was  operated  only  once. 

In  practical  applications,  sometimes  it  is  necessary  to  carry  out  continuous 
photographing  of  the  measured  target,  i.e.,  one  hopes  to  obtain  many  photo¬ 
graphic  frames.  At  such  times  it  is  necessary  to  add  additional  control  elec¬ 
trodes  and  corresponding  control  circuits  to  the  foundation  of  the  image- 
converter  described  above.  For  example,  for  a  shutter  type  image— converter , 
after  each  exposure,  a  step  voltage  deflecting  electron  ray  can  go  to  another 
position  on  the  fluorescent  screen  to  make  the  next  exposure .  In  this  way , 
each  exposure  not  only  is  staggered  in  time  but  in  space  is  also  located  in 
a  different  place  on  the  fluorescent  screen.  After  many  exposures,  a  photo¬ 
graph  of  a  certain  size  can  be  obtained.  This  multiple  frame  operation  method 
has  very  important  significance  in  actual  applications.  For  example,  the 
changes  in  attitude  of  a  flight  body  (such  as  a  guided  missile)  and  the  changes 
in  spark  gap  discharge  channels,  etc.,  can  be  photographed.  Of  course,  the 
more  pictures  a  fluorescent  screen  of  fixed  dimensions  is  divided  into,  the 
smaller  the  images.  Reference  9  also  describes  a  new  image  dissection  type 
special  image  division  image-converter,  which  can  record  three  frames  in  one 
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image,  with  a  division  cycle  smaller  than  100  ps.  The  picture  division  camera 
using  this  image-converter  has  been  used  in  laser  fusion  research.  To  use 
scanning  streak  cameras  to  study  highly  repetitive  high-speed  optical  phenom¬ 
ena  (such  as  checking  continuous  mode -locked  laser  pulse  sequence)  high 
frequency  sine  wave  voltage  can  be  used  to  replace  linear  slope  wave  voltage 
as  the  streak  image-converter's  scan  signal.  This  scanning  mode  is  called 
"synchroscan. Present  standard  commercial  scanning  streak  cameras  (such 
as  the  IMA.C0N  675)  is  already  equipped  with  this  type  of  S5mchronous  scanning 
device . 

Since  ultrafast  image-converter  cameras  with  MCP  have  extremely  high  optical 
gain  and  extremely  short  time  resolution  they  have  become  indispensable  tools 
for  high-speed  phenomena  research  and  their  range  of  applications  involves 
many  Important  technical  areas.  In  addition  to  the  above  described  applica¬ 
tions  of  ultrafast  image-converter  cameras  in  gas  breakdown,  laser  fusion  and 
ultrashort  light  pulse  measurement,  they  have  also  been  applied  in  such  areas 
as  measurement  of  laser  produced  plasma,  research  on  molecular  activated 
states,  X-ray  diagnosis,  nonlinear  optics,  picosecond  spectroscopy,  fluorescent 
spectroscopy,  synchronous  accelerator  radiation,  explosions,  optical  communica¬ 
tions  and  optical  information  processing  and  photobiology,  photochemistry,  and 
photophysics.  Of  course,  to  adapt  to  different  application  goals,  it  is  neces¬ 
sary  to  adapt  and  expand  the  basic  photographic  system  illustrated  in  Figure  5 
and  the  image-converter  type  used  should  also  be  selected  with  a  view  to  the 
research  object,  but  as  far  as  the  basic  problems  involved  with  image-converter 
cameras  is  concerned,  generally  they  are  as  described  in  this  paper. 
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BEIJING  ELECTRON-POSITRON  COLLIDER  PROJECT  DETAILED 

Hong  Kong  LIAOWANG  OVERSEAS  EDITION  in  Chinese  No  14,  7  Apr  86  pp  4-5 

[Article  by  Gu  Hainan]  [7357  6701  0589]] 

[Text]  The  Beijing  Electron-Positron  Collider  (BEPC)  project  which  attracted 
a  great  deal  of  attention  from  the  international  scientific  community  has  been 
under  construction  for  over  2  years.  Recently,  this  reporter  traveled  to  the 
construction  site  and  had  an  interview  with  the  three  gentlemen  in  charge  of 
the  project:  Gu  Yu,  Xie  Jialin,  and  Fang  Shouxian. 

In  the  project  coordination  room,  group  leader  Gu  Yu  briefly  summarized  the 
current  project  status.  For  detailed  questions,  he  referred  the  author  to 
professors  Xie  Jialin  and  Fang  Shouxian,  who  are  experts  in  particle 
accelerators . 

Pointing  to  a  large  BEPC  schedule  and  progress  chart  on  the  wall,  the  two 
experts  said:  "Since  October  1984,  when  director  Deng  Xiaoping  personally 
endorsed  the  project,  dozens  of  organizations  and  thousands  of  engineering  and 
technical  personnel  had  dedicated  their  efforts  to  the  project;  during  this 
time,  the  prototypes  of  six  key  facilities  used  for  the  collider  had  been 
developed  and  certified;  other  equipment  is  also  being  developed.  The  near 
1-km  tunnel  is  partially  finished.  In  short,  the  completion  date  of  the 
entire  project  is  nearly  in  sight.  A  large  national  laboratory  which  will 
attract  physicists  around  the  world  to  explore  the  secrets  of  the  microscopic 
world  will  soon  be  erected  on  Chinese  soil! 

In  the  office  of  the  front-line  project  director  Fang  Shouxian,  the  reporter 
noticed  a  lunch  box  on  the  desk  and  blankets  on  the  sofa. 

"Do  you  use  your  office  as  your  living  quarters?" 

"To  build  such  a  high-technology  facility  is  not  only  unprecedented  in  China 
but  also  a  very  challenging  task  Internationally.  Director  Deng  Xiaoping  had 
requested  that  our  accelerator  must  be  completed  on  or  ahead  of  schedule.  Now^ 
that  the  completion  date  is  in  sight  and  this  gold  medal  for  the  state  is 
nearly  at  hand,  the  entire  staff  on  this  project  are  working  around  the  clock. 
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"What  is  an  electron— positron  collider?  What  Is  It  used  for?"  These  are 
probably  the  questions  that  most  readers  would  like  to  know  the  answers  to. 

Prof.  Xle  and  Prof.  Fang  gave  the  following  explanation:  "A  collider  Is  a 
large  research  apparatus  which  consists  of  thousands  of  precision  parts 
connected  together  In  series.  When  the  power  Is  turned  on,  each  component 
must  perform  In  exactly  the  prescribed  manner.  Two  beams  of  electrons  and 
positrons  each  0.1  mm  thick,  1  mm  wide,  and  10  cm  long  revolve  In  the 
accelerator  along  a  240-m  acceleration  tube.  At  every  revolution  a  collision 
occurs,  and  there  are  more  than  1  million  collisions  per  second.  In  order  to 
achieve  collision,  the  collider  design  must  meet  very  stringent  requirements. 

This  collider  Is  expected  to  be  ready  for  beam  generation  by  the  end  of  next 
year,  and  the  first  electron-positron  collisions  are  expected  to  take  place  In 
1988.  The  collider  consists  of  the  following  components:  the  Injector,  the 
storage  ring,  the  detector,  the  synchronous  radiation  test  region,  and  the 
computing  center.  The  basic  objectives  for  building  the  collider  were:  1) 
to  promote  the  research  of  high-energy  physics  In  China,  and  2)  to  promote 
the  application  of  synchronous  radiation. 

Sketching  a  schematic  diagram  on  a  piece  of  paper,  Xle  Jlalln  proceeded  to 
explain  the  extreme  Importance  of  basic  research.  The  discovery  of 
television,  radar.  X-ray's,  the  atomic  bomb  and  atomic  energy  all  hinged  on  a 
good  understanding  of  the  structure  of  atoms  and  nuclei.  The  appearance  of 
these  things  unquestionably  had  profound  effect  on  human  lives.  The  particles 
that  make  up  the  microscopic  structure  of  matter  are  so  small  that  they  are 
not  visible  to  human  eyes  or  microscopes.  Therefore,  scientists  Invented 
different  types  of  accelerators  to  study  these  "fundamental  particles" — the 
smallest  units  of  matter.  The  collider  Is  a  new  type  of  accelerator  which  was 
developed  In  the  1970 's.  Everytlme  a  breakthrough  occurs  In  these  studies, 
human  knowledge  about  nature  takes  a  giant  step  forward.  In  the  early  1970' s, 
the  well-known  Chlnese-Amerlcan  physicist  Ding  Zhaozhong  discovered  the  J 
particle  using  a  proton  cyclotron  at  the  Brooke  Haven  Laboratory.  At  about 
the  same  time,  American  scientist  (Llchter)  also  found  this  particle  In  a 
positron-electron  collider  at  the  Stanford  Linear  Accelerator  Center.  In  the 
latter  case,  the  signal  was  loud  and  clear,  and  the  noise  background  was  low. 
This  Illustrates  the  superiority  of  the  electron  collider. 

The  collider  can  also  be  used  In  the  study  of  the  application  of  synchronous 
radiation.  Synchronous  radiation  Is  the  light  emitted  by  high-energy 
electrons  traveling  along  a  circular  orbit  in  a  synchronous  cyclotron.  This 
light  source  provides  the  basic  tool  for  studying  a  variety  of  scientific 
topics  such  as  surface  physics,  atomic-molecular  physics,  solid-state  physics, 
material  science,  biochemistry,  biophysics,  genetic  engineering,  catalytic 
chemistry,  environmental  science,  medicine.  X-ray  microscopy  and 
photo-etching . 

\ 

At  this  point.  Prof.  Xle  showed  us  a  silicon  chip  covered  with  crisscrossing 
Integrated  circuits  and  said:  ”The  connections  of  today’s  Integrated  circuits 
are  several  microns  wide;  by  using  synchronous  radiation,  the  width  of  the 
connections  can  be  reduced  to  less  than  1  micron,  thus  greatly  Increasing  the 
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level  of  integration.  Synchronous  radiation  can  also  be  used  in  medicine.  To 
detect  the  blockage  of  coronary  arteries  of  elderly  persons  because  of 
cholestral  accumulation,  one  of  the  methods  is  to  take  two  pictures  of  the 
heart  with  adjacent  wavelengths  using  synchronous  radiation,  and  process  the 
information  using  computers.  Thus,  the  patient  runs  no  risk  of  surgery,  and 
is  exposed  to  very  little  radiation  dosage.  Therefore,  it  could  very  well 
become  a  widely  used  tool  for  treating  arteriosclerosis. 

Today  there  are  more  than  10  positron-electron  colliders  in  the  world.  The 
United  States  is  planning  to  spend  6  billion  dollars  over  a  10— year  period  to 
build  a  100— km  super-conducting  super  collider  (SSC).  The  European  Nuclear 
Research  Center  is  in  the  process  of  constructing  a  27-km  positron-electron 
collider.  China's  current  economic  resources  do  not  allow  expenditure  of  this 
magnitude.  But  while  Its  energy  level  Is  rather  moderate,  our  collider  will 
have  a  unique  place  In  the  International  community  because  the  energy  level  of 
this  collider  (2.2  billion  eV  to  2.8  billion  eV)  allows  the  study  of  certain 
problems  in  physics  that  cannot  be  investigated  with  other  types  of 
accelecators.  At  present,  our  collider  has  the  highest  design  brightness  In 
the  world  in  this  energy  region. 

The  director  of  the  project  office,  Liu  Hualzu  explained  that  the  project 
guidance  group  organized  by  the  government  is  composed  of  members  from  the 
Chinese  Academy  of  Sciences,  the  National  Planning  Commission,  the  National 
Economic  Commission,  and  officials  of  Beijing  City.  This  is  a  group  with 
clearly  defined  responsibilities  and  high  operating  efficiency;  It  not  only 
provides  the  day-to-day  project  guidance  from  the  Chinese  Academy  of  Sciences, 
but  is  also  responsible  for  coordinating  the  activities  of  nearly  100  units 
from  more  than  10  ministries  and  committees  around  the  country.  Because  of 
the  good  coordination  effort,  the  project  is  progressing  rapidly.  He  pointed 
out  that  of  the  nearly  10,000  parts  and  components  of  the  collider,  70  percent 
are  non-standard,  special-purpose  devices  which  must  meet  high  precision 
requirements,  and  are  technically  very  different  to  build.  But  the  prototypes 
of  most  of  these  special-purpose  devices  had  been  developed  in  about  a  year, 
and  they  are  already  being  mass-produced.  The  civil  construction  work  Is  also 
moving  at  a  rapid  pace.  The  200— m  tunnel  of  the  electron  linear  accelerator 
and  the  large  ^11  containing  the  accelerator  tube  have  already  been 
completed.  Such  speed  of  construction  is  unprecedented  in  China  and 
considered  above  average  based  on  international  standards.  On  this  project, 
we  not  only  took  full  advantage  of  the  most  advanced  technologies  available  to 
us,  but  also  saved  considerable  amount  of  money  for  the  state. 

After  the  interview,  the  reporter  was  given  a  tour  of  the  construction  site. 
Directors  of  Research  Yan  Talxuan  and  Luo  Tlngxlung  pointed  to  a  large  curved 
magnet  being  tested  and  said:  "Each  of  these  magnets  is  made  of  3,000  pieces 
of  0.5-mm  thick  silicon  steel  plates  pressed  together.  The  collider  requires 
40  of  such  magnets.  Technicians  at  the  Xian  Feng  Electric  Motor  Co.  in 
Shanghai  can  produce  3  to  4  high-quality  magnets  per  month,  and  they  are 
working  around  the  clock.  In  the  past,  we  never  attempted  to  make  this  t3rpe 
of  magnet.  In  order  to  guide  the  electrons  along  a  prescribed  track,  the 
collider  also  needs  60  focusing  magnets.  These  magnets  must  satisfy  even 
higher  precision  requirements,  and  each  piece  is  pressed  together  with  1,00 
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silicon  steel  plates.  They  are  being  built  by  workers  of  the  High  Energy 
Institute,  and  are  scheduled  to  be  Installed  in  the  tunnel  by  the  end  of  this 
year.  In  addition,  56  precision  acceleration  tubes  which  are  3-m  long  and 
made  of  oxygen-free  copper,  16  high-power  speed  regulation  tubes,  modulators, 
and  banks  of  linear  accelerators  are  also  needed.  Therefore,  this  is  truly  a 
rare  scientific  and  technological  challenge  for  this  country.  To  ensure  that 
the  project  will  be  on  schedule,  thousands  of  technical  personnel  and  workers 
have  devoted  their  lives  to  the  project;  some  scientists  give  up  the 
opportunity  to  study  abroad  and  spend  their  holidays  working  on  the  project. 

Speaking  of  the  enthusiasm  exhibited  on  this  project.  Fang  Shouxlan  said: 

"Many  of  us  have  worked  abroad  because  we  did  not  have  our  own  accelerator  and 
had  to  use  other  country's  accelerators  for  research.  Since  returning  home 
from  the  Dubna  Joint  Institute  for  Nuclear  Research  of  the  Soviet  Union,  I 
have  participated  in  the  design  of  several  accelerators  with  my  colleagues 
over  the  past  20  years,  but  all  of  them  were  abandoned  because  of  lack  of 
adequate  preparation.  Therefore,  the  colleagues  of  my  generation  are 
determined  to  make  our  dream  come  true,  l.e.,  to  successfully  build  an 
accelerator  within  our  lifetime  so  as  to  boost  the  pride  of  the  Chinese  people 
and  to  reserve  a  place  for  China's  high-energy  physics  on  the  international 
stage." 
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THE  TIME  PROJECTION  CHAMBER:  NEW  CONCEPT  IN  PARTICLE  DETECTION 

Shanghai  HE  JISHU  [NUCLEAR  TECHNIQUES]  In  Chinese  No  5,  May  86  pp  19-23,  37 

[Article  by  Chen  Chaoqing  [7115  2600  3237]] 

[Text]  I.  Introduction 

The  time  projection  chamber  (TPC)^  is  a  new  detection  device  that  can  simul¬ 
taneously  measure  the  three-dimensional  position,  momentum  of  charged  parti¬ 
cles  and  perform  discrimination  of  particle  types.  It  differs  from  a 
conventional  cylindrical  drift  chamber  in  that  inside  a  TPC,  the  electric 
field  is  parallel  to  the  magnetic  field,  thus  the  force  E  x  B  has  little 
effect  on  electron  drift;  the  detection  surfaces  are  located  on  the  two  end 
covers,  and  because  there  are  no  metallic  wires  along  the  axial  direction  of 
the  cylinder,  the  sensitivity  volume  can  be  made  very  large,  which  can 
accommodate  several  hundred  sample  measurements  of  the  energy  loss  (dE/dx), 
thus  providing  good  discrimination  capability.  It  has  extremely  good  spatial 
and  momentum  resolution,  dE/dx  resolution,  and  facilitates  pattern  recogni¬ 
tion;  it  is  a  relatively  advanced  high-energy  detector  which  can  be  used  in 
collision  physics  experiment  or  stationary  target  experiment.  The  same  prin¬ 
ciple  can  also  be  used  to  build  a  simpler  TPC  which  would  have  sufficiently 
good  performance  to  be  used  as  a  powerful  tool  in  the  study  of  medium-energy 
physics,  nuclear  physics,  and  heavy  ions. 

The  basic  concept  of  TPC  was  proposed  in  1974  by  Dr  D.R.  Nygren^  of  the 
Lawrence  Berkeley  Laboratory  (LBL) ,  and  it  was  first  developed  and  built  at 
the  LBL.^“®  Over  the  last  10  years,  many  laboratories  around  the  world  had 
designed  and  built  various  types  of  TPC’s  based  on  the  requirements  of  the 
individual  physics  experiments. 

Under  the  guidance  of  Dr  Nygren,  the  author  participated  in  part  of  the 
construction  work  and  testing  of  the  LBL-TPC.  This  article  briefly  describes 
the  structure,  operating  principle,  and  performance  test  results  of  the  LBL- 
TPC. 


II.  The  Structure  of  LBL-TPC 

Figure  1  shows  a  structural  diagram  of  the  time  projection  chamber.  It  is  a 
special  cylindrical  drift  chamber  which  Is  2m  long  and  has  an  inside  diameter 
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of  40cm  and  an  outside  diamter  of  2m,  The  cylinder  is  divided  into  two  com¬ 
pletely  symmetric  sections  by  a  thin-film  electrode.  The  thin-film  electrode 
is  a  mesh  electrode  made  of  stainless  steel  wires}  it  offers  low  resistance 
to  airflow  and  is  transparent  to  particles.  A  high  negative  voltage  is 
applied  to  the  electrode;  the  maximum  design  value  is  -200  kV,  The  inside 
and  outside  walls  of  the  cylinder  are  covered  with  copper— plated  annular 
electrodes,  called  the  electric  field  cage.  A  resistance  voltage  divider  is 
used  to  form  a  uniform  drift  field. 
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Figure  1.  Structural  Diagram  of  a  Time  Projection  Chamber 


The  two  end  covers  of  the  cylinder  each  has  six  pie-shaped  multiwire  positive- 
ratio  chambers  (MWPC).  Each  MWPC  contains  four  planar  electrodes:  the  cathode, 
the  anode,  the  screen  grid,  and  the  gate  grid,  as  shown  in  Figure  2.  The  dis¬ 
tance  between  the  anode  and  the  cathode  are  both  4mm;  the  distance  between  the 
gate  grid  and  the  screen  grid  is  8mm.  The  grid  mesh  is  made  of  76  ym 
beryllium  copper  wire  which  is  separated  by  1  millimeter.  The  screen  grid 
wires  perform  two  functions:  one  is  to  separate  the  drift  region  from  the 
amplification  region;  the  other  is  a  funnel-like  function,  which  "focuses"  the 
electric  field  lines  onto  the  signal  wires.  Thus,  the  ionizing  electrons  pro¬ 
duced  by  each  of  the  4mm-long  particle  tracks  are  drifted  to  the  corresponding 
signal  wire.  The  gate  grid^*'  is  used  to  prevent  the  feedback  of  positive  ions 
and  to  reduce  the  effect  of  space  charge.  It  performs  the  function  of  an  elec¬ 
tronic  switch.  Normally,  this  gate  is  always  closed,  which  blocks  the  drift 
electrons;  the  gate  is  opened  only  when  a  favorable  event  occurs  to  let  the 
drift  electrons  reach  the  anode  mesh  in  order  to_induce  an  avalanche.  To 
reach  this  goal,  a  common  negative  bias  voltage  Vq  is  applied  to  all  the 
gate  grid  wires,  so  that  the  drift  field  lines  can  freely  pass  through  the  gate 
grid,  as  shown  in  Figure  3a;  to  close  the  gate  grid,  two  low  voltages  with 
opposite  gate  grid  wires  to  prevent  the  electric  field  lines  and  hence  the 
drift  electrons  from  passing  through  the  gate  grid,  as  shown  in  Figure  3b.  The 
opening  and  closing  of  the  gate  grid  has  183  signal  wires  and  an  equal  number 
of  electric  field  wires.  The  signal  wires  are  made  of  20  ym  gold-plated 
tungsten  wires, and  are  used  to  measure  the  energy  loss  dE/dx.  The  distance 
between  signal  wires  is  4  millimeters.  The  electric  field  wires  are  made  of 
76  ym  beryllium-copper  wires;  they  are  placed  between  neighboring  signal  wires 
to  reduce  crosstalk.  The  cathode  plane  is  made  of  1  cm-thick  plates  with 
copper  coatings  on  both  sides.  Fifteen  rows  of  cathode  pads  are  etched  on 
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Figure  3.  Switching  States  and  Electric  Field  Distribution  of  the  Gate  Grid 

the  copper-coated  plates,  each  pad  has  an  area  of  7.5inm  x  7.5inm  and  the  insu¬ 
lated  gap  is  0.5  millimeters.  These  15  rows  of  pads  are  located  directly  be¬ 
low  the  corresponding  15  signal  wires  to  measure  the  coordinates  along  the 
wires . 

Each  pie-shaped  MWPC  carries  1,335  low-noise  preamplifiers  which  are  installed 
on  the  back  side  of  the  cathode  plate.  In  order  to  lower  the  temperature  rise 
caused  by  preamplification,  each  MWPC  is  equipped  with  a  water  cooling  system 
on  the  back  side. 

The  entire  TPC  is  housed  inside  a  high-pressure  container  filled  with  a  mix¬ 
ture  of  argon  (80  percent)  and  methane  (20  percent) .  The  maximum  design 
pressure  is  10  atm.  The  high-pressure  container  is  enclosed  by  a  magnetic 
coil. 
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III.  Operating  Principle 

When  charged  particles  pass  through  the  chamber,  an  ionized  track  is  formed 
in  the  gas.  Under  the  action  of  the  electric  field,  the  electrons  are  drifted 
toward  the  MWPC  of  the  end  cover.  During  the  drift  process,  lateral  disper¬ 
sion  of  the  electrons  is  constrained  by  the  axial  magnetic  field  which  is 
parallel  to  the  electric  field;  hence  the  electrons  can  drift  over  a  long 
distance  without  significant  dispersion.  The  amount  of  dispersion  is  deter¬ 
mined  by  the  magnetic  field  strength  B  and  the  gas  pressure  P; 

ctdCB)  =<7d(0)/V1  +  (ot)' 

where  OpCB)  and  c^(0)  are,  respectively,  the  RMS  lateral  dispersion  distance 
with  and  without  magnetic  field,  w  =  eB^m  is  the  gyromagnetlc  ratio  of  the 
electrons,  T  is  the  mean  collision  time  of  the  electrons,  which  is  inversely 
proportional  to  the  gas  pressure  P.  Clearly,  when  the  magnetic  field  strength 
B  is  sufficiently  large  so  that  ut>1,  then  lateral  dispersion  can  be  effec¬ 
tively  compressed. 

When  the  drift  electrons  reach  the  MWPC  of  the  end  covers,  ionization  avalanche 
occurs  in  the  vicinity  of  the  anode  wire,  thus  generating  a  negative  pulse 
signal.  At  the  same  time,  a  positive  pulse  signal  is  induced  on  the  pad  below 
the  anode  wire,  as  shown  in  Figure  4.  Each  avalanche  on  the  signal  wire  will 
produce  an  induced  signal  above  the  threshold  on  2  to  3  pads  below  the  signal 
wire.  The  magnitude  of  the  signal  is  proportional  to  the  solid  angle  subtended 
by  the  pad  relative  to  the  position  of  avalanche.  The  barycenter  of  the 
Induced  signal  determines  the  position  of  the  avalanche  on  the  signal  wire. 


Figure  4. 
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After  amplification  by  the  low-noise  FET  preamplifier,  the  negative  pulse 
generated  on  the  signal  wire  and  the  positive  pulse  generated  on  the  pad  are 
transmitted  to  the  remote  shaping  amplifier®  where  they  are  subject  to  further 
amplification  and  discrimination,  and  then  stored  into  a  charge  couple  device 
(CCD).®  CCD  is  an  analog  shift  storage  device  that  contains  455  charged 
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storage  elements.  Under  the  control  of  external  timing  pulses,  the  charge 
(pulse  height)  signal  moves  from  one  storage  element  to  another .  The  fre¬ 
quency  of  the  external  timing  pulse  at  write-in  time  is  10  MHz,  hence  the 
time  for  a  signal  to  move  from  one  element  to  another  is  100  ns;  a  CCD  can 
store  an  analog  signal  45.5  ys  in  duration.  The  maximum  drift  distance  in 
a  TPC  is  1  m,  the  electron  drift  speed  is  5cm/ys,  hence  the  maximum  drift 
time  is  20  ys.  The  delay  time  of  a  CCD  is  capable  of  storing  all  the  analog 
signals  within  the  drift  time  of  the  TPC.  At  read-out  time,  the  frequency 
of  the  external  timing  pulse  drops  to  110  KHz  to  facilitate  digital  conver¬ 
sion  of  the  analog  signals.  The  converted  data  are  stored  in  the  buffer 
storage  unit  for  computer  processing. 

1 .  Coordinate  Measurement 

We  define  the  X  coordinate  to  be  along  the  direction  of  the  signal  wire,  the 
Y  coordinate  to  be  in  the  anode  plane  and  perpendicular  to  the  signal  wire, and 
the  Z  coordinate  to  be  along  the  axial  direction  of  the  chamber.  The  X  and  Y 
coordinates  of  the  particle  track  are  determined  by  the  position  of  the  bary- 
center  of  the  signal  on  the  cathode  pad,  and  the  predicted  measurement  accuracy 
(RMS  deviation)  is  150  ym.  The  Z  coordinate  is  determined  by  the  drift  time 
of  the  ionizing  electrons  from  the  particle  track  to  the  MWPC  of  the  end  cover; 
the  predicted  measurement  accuracy  is  350  ym. 

2.  Momentum  Measurement 

Having  determined  the  three-dimensional  coordinates  of  the  charged  particle 
track,  one  derives  its  radius  of  curvature  R  on  a  plane  perpendicular  to  the 
magnetic  field,  and  determines  the  lateral  momentum  of  the  particles  P  =  0.03 
BR  (where  the  unit  of  P  is  GeV/C,  the  unit  of  the  magnetic  field  is  T,  and  the 
unit  of  the  radius  of  curvature  is  m) .  The  range  of  momentum  that  can  be 
measured  by  the  LBL-TPC  is  from  100  MeV/C  to  15  GeV/C.  When  the  magnetic  field 
strength  is  1.5  T,  the  predicted  momentum  resolution  is  AP/P^  =0.4  percent 
(GeV/C)"^. 

3.  Particle  Discrimination 

The  output  pulse  height  from  the  signal  wire  indicates  the  magnitude  of  the 
energy  loss  dE/dx  suffered  by  the  particles  in  passing  through  the  4mm-high 
pressure  gas.  The  energy  loss  dE/dx  of  a  charge  particle  is  related  to  its 
speed.  Once  the  particle  momentum  is  determined,  the  charge-to-mass  ratio 
of  the  particle  can  be  determined  by  measuring  its  energy  loss  dE/dx,  hence 
the  type  of  particle  can  be  identified. 

Figure  5  shows  the  relationship  between  the  most  probable  energy  loss  and 
momentum  for  several  different  types  of  particles  over  a  distance  of  4mm  in 
an  8.5-atm  gas  mixture  of  argon  (80  percent)  and  methane  (20  percent).  It 
can  be  seen  that  as  the  momentum  increases,  the  particle  energy  loss  can  be 
divided  into  three  different  regions, i.e.,  the  nonrelativist ic  region 
(BY<2)>  the  relativistic  region  (2<3y<100) ,  and  the  super  relativistic 
region  (By»100).  Here  3  =  v/c,  y  =  (1  ~  V  is  the  particle  velocity, 

and  c  is  the  speed  of  light  in  vacuum.  In  the  nonrelativistic  region,  dE/dx 
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Figure  5.  Relationship  Between  the  Most  Probable  Energy  Loss  of 
Charged  Particle  and  Its  Momentum 
Mote:  The  gas  mixture  contains  80  percent  argon  and 
20  percent  methane,  the  pressure  is  8.5  atm, 
and  the  sample  thickness  is  4mm 

varies  with  3”^,  and  particle  discrimination  is  relatively  easy.  In  the 
super  relativistic  region,  dE/dx  does  not  vary  with  momentum,  hence  one  can¬ 
not  use  dE/dx  for  particle  discrimination.  In  the  relativistic  region,  dE/dx 
Increases  logarithmically  with  3y*  Figure  5  shows  that  in  this  region  (where 
momentum  varies  from  1  GeV/c  to  several  hundred  GeV/c),  the  rise  in  energy 
loss  is  only  around  40  percent.  Although  the  energy  losses  for  different 
types  of  particles  at  the  same  momentum  value  are  not  the  same,  the  differ¬ 
ences  are  quite  small.  Therefore,  if  the  particles  are  to  be  distinguished 
by  the  relativistic  rise  in  ionization  energy  loss,  the  detector  must  have 
extremely  high  dE/dx  resolution. 

But  in  a  thin  gas  medium,  there  are  considerable  fluctuations  in  the  particle 
energy  loss,  which  is  referred  to  as  the  "Landau  tail."  It  is  formed  by  the 
energy  spectrum  of  the  "6  electrons"  produced  by  the  particles  passing  through 
the  medium.  Therefore,  the  dE/dx  resolution  of  a  single  sample  is  very  poor. 
The  key  to  improving  dE/dx  resolution  is  to  perform  multiple  independent 
measurements  of  the  energy  loss  of  each  particle.  In  general,  if  the  number 
of  samples  is  N,  then  the  resolution  can  be  improved  by  a  factor  of  ✓».  Each 
pie-shaped  MWPC  on  the  end  cover  has  183  signal  wires,  hence  as  many  as  183 
dE/dx  samples  can  be  taken  along  each  track.  Furthermore,  the  high-pressure 
gas  of  a  TPC  ensures  that  each  sample  contains  a  sufficient  number  of  ionizing 
electrons,  and  that  the  dE/dx  resolution  of  a  TPC  can  exceed  the  3  percent 
requirement  so  that  within  the  momentum  range  from  100  MeV/c  to  15  GeV/c,  it 
can  discriminate  such  particles  as  fr  mesons,  K  mesons  and  protons  as  well  as 
the  low-momentum  electrons  and  y  mesons. 
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IV.  Performance  Results 

1.  Spatial  Resolution^** 

The  spatial  resolution  Ojjy  on  the  XY  plane  depends  on  many  factors.  Based  on 
the  results  of  cosmic  ray  measurements,  this  relationship  can  be  expressed  in 
the  following  form; 

Oxr  =<^A  +  Og-^^^—e'^-seca  +  ale'^-cosatan^a 
I'M 

where  L  is  the  drift  distance  (m) ;  is  the  maximum  drift  distance  (Im) ; 
a  is  the  subtended  angle  between  the  projected  track  on  the  MWPC  and  the  Y 
axis;  y  is  the  electron  capture  rate  per  meter.  The  first  term  represents 
the  inherent  spatial  resolution  of  the  chamber,  which  is  determined  by  the 
electronic  noise  and  other  system  effects.  The  second  term  is  due  to  electron 
lateral  dispersion,  and  the  third  term  is  due  to  ionization  fluctuations. 

Measurements  are  made  in  the  chamber  with  a  gas  mixture  of  80  percent  argon 
and  20  percent  methane  under  the  following  conditions:  gas  pressure — 

8.5  atm,  drift  electric  field  strength — 75  kV/m,  and  magnetic  field  strength — 
0.4T;  the  results  are  =  (160±2)  ym,  Og  =  (105±6)  ym,  Oq  =  (249±7)  ym. 
Therefore,  when  L  =  0  pnd  a  =  0®,  c^y  =  160  ym,  which  is  very  close  to  the 
predicted  value  of  150  ym.  If  a  superconducting  magnet  (1.5T)  is  used  instead 
of  the  conventional  magnet  (0.4T),  further  improvement  can  be  expected. 

The  measured  resolution  in  the  z  direction  o^,  is  (340±5)  ym,  which  is  better 
than  the  predicted  value  of  350  ym.  The  resolution  of  the  dual  tracks  is 
1  cm. 


2.  Momentum  Resolution 

Under  the  above  conditions,  the  momentum  resolution  based  on  cosmic  ray 
measurements  is: 

dP/P^  =  (6.410.3)  percent  (GeV/c)“^ 

For  a  magnetic  field  strength  of  0i4T,  the  predicted  momentum  resolution  is 
1.6  percent  (GeV/c)“^,  which  differs  from  the  measured  value  by  approximately 
a  factor  of  4.  This  is  due  to  the  imperfect  design  of  the  chamber  which 
causes  the  particle  tracks  to  be  twisted  in  the  vicinity  of  the  inner  and 
outer  walls  because  of  electric  field  distortion.  Consequently,  in  the  data 
analysis  only  the  untwisted  segment  of  the  track  can  be  used,  which  degrades 
the  momentum  resolution.  Later,  improvements  have  been  made  to  the  electric 
field  cage,  and  a  gate  grid  has  been  added  to  the  MWPC  on  the  end  cover,  thus 
eliminating  the  twisting  effect  and  greatly  improving  the  momentum  resolution. 

3.  Resolution  of  dE/dx^®»^^ 

Measurements  are  made  using  cosmic  rays.  A  particle  track  with  155-183  dE/dx 
samples  is  chosen,  but  in  the  data  analysis,  only  65  percent  of  the  samples 
with  the  lowest  pulse  height  are  included  in  the  averaging  process.  Under  the 
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pressures  of  8.5,  4.0,  and  1.0  atm,  the  measured  dE/dx  resolution 
respectively  3.0±0.1  percent,  4.210.3  percent,  and  6.510.5  percent. 

4.  Measurement  of  Relativistic  Rise 

The  relativistic  rise  of  dE/dx  is  measured  using  the  y  mesons  of  the  cosmic 
rays.  The  relativistic  rise  is  represented  by  the  ratio  between  the  3.5  GeV/c 
K-TT  separation  and  the  minimum  ionization.  The  energy  loss  of  the  ir  meson  and 
K  meson  with  a  momentum  value  of  3.5  GeV/c  is  equivalent  to  the  energy  loss  of 
y  mesons  with  momentum  values  of  2.65  GeV/c  and  0.75  GeV/c  in  the  same  gas 
medium: 

K-TT  separation  =  [  (dE/dx)  (2.65GeV/c)-(dE/dx)  (0.75GeV/c)  ]  / 

[ (dE/dx) (minimum  value) ] 

The  measured  results  are  listed  in  the  table  below.  If  we  want  to  achieve  a 
4“Standard  deviation  K-7T  separation,  the  dE/dx  resolution  must  be  better  than 
3.5  percent,  and  the  systematic  error  must  be  smaller  than  1  percent.  There¬ 
fore,  fluctuations  in  the  signal  wire  gain  with  voltage,  gas  pressure,  and 
temperature  must  be  carefully  controlled  and  corrected. 


Table 


Gas  pressure  (atm) 

3.5GeV/C 

k/tt  separation 

Percent 

Standard 

deviation 

8.5 

14.6±0.3 

4.810.25 

4.0 

16 .6±0.5 

3.910.5 

1.0 

18.410.3 

2.810.3 
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[Text]  I.  Introduction 

The  activities  of  nuclear-powered  ships  reach  every  port  around  the  world; 
but  today's  technology  cannot  completely  eliminate  the  possibility  of  an 
accident  which  would  release  radioactive  materials  into  the  environment. 
Consequently,  the  operational  safety  of  a  nuclear  ship  is  an  issue  of  world¬ 
wide  concern.  If  the  nuclear  propulsion  system  is  to  become  a  reliable  power 
source,  not  only  must  its  utilization  efficiency  be  improved,  but  the  safety 
problem  must  be  solved.  The  safety  requirements  of  a  land-based  nuclear  power 
plant  generally  need  only  be  approved  by  the  country  using  the  power  plant, 
but  the  safety  requirements  of  a  nuclear  ship  must  be  approved  by  a  world 
organization.  Therefore,  a  more  stringent  requirement  must  be  imposed  on 
the  safety  of  nuclear  propulsion  systems;  this  requirement  must  be  such  that 
"the  use  of  this  system  will  not  produce  greater  danger  to  the  crew  members, 
the  passengers,  and  the  environment  than  a  ship  using  conventional  propulsion 
system."  Otherwise,  nuclear  ships  may  ultimately  be  forced  out  of  the  ship¬ 
ping  business  due  to  various  limitations  imposed  on  the  design  and  operation 
in  the  name  of  safety. 

This  article  briefly  addresses  several  important  problems  that  must  be  under¬ 
stood  in  carrying  out  the  safety  design  of  nuclear  propulsion  systems,  the 
safety  concept  of  nuclear-powered  ships  and  the  nuclear  propulsion  system,  as 
well  as  the  key  methods  used  in  safety  design. 

II.  Safety  Concept  of  Nuclear-Powered  Ships 

As  a  safety  system,  the  "safety  concept  of  nuclear-powered  ships"  is  developed 
from  a  combination  of  the  traditional  safety  concept  of  sea-going  ships  and 
the  safety  concept  of  nuclear  power  plants. 

The  traditional  safety  concept  of  sea-going  ships  is  as  follows;  most 
disasters  and  accidents  over  the  ocean  produce  primarily  local  consequences; 
they  mainly  endanger  the  ship  Itself,  the  ship  crew  and  passengers,  but  have 
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very  little  effect  on  the  general  public  and  the  environment.  Therefore,  the 
basic  concept  of  treating  conventional  ship  safety  problems  is  "survival  of 
the  ship  is  the  ultimate  guarantee  of  the  crew  members’  lives;  protecting  the 
ship  from  destruction  means  protecting  the  environment,"  i.e.,  "the  ship’s 
safety  comes  first."  The  guidelines  derived  from  this  concept  for  ship 
design  and  operation  are:  1)  minimize  the  probability  of  getting  into  an 
accident;  2)  mitigate  the  consequences  of  a  ship  disaster  and  minimize  the 
probability  of  ship  destruction;  3)  increase  the  rate  of  survival  and  the 
rate  of  rescue  in  case  of  a  ship  destruction  and  abandonment. 

The  basic  concept  of  treating  safety  problems  in  nuclear  systems  is:  "the 
nuclear  reactor’s  safety  comes  first."  It  requires  that  "under  normal  opera¬ 
tion  or  in  case' of  an  accident,  the  amount  of  radioactive  materials  released 
by  the  nuclear  devices  shall  not  exceed  a  level  that  will  endanger  the 
health  of  the  residents  or  cause  pollution  to  the  environment."  Therefore, 
in  the  design  and  operation  of  a  nuclear  facility,  1)  we  must  minimize  the 
probability  of  accident  occurrence;  2)  once  an  accident  has  occurred,  we 
must  first  protect  the  reactor  from  being  damaged;  3)  if  the  reactor  is 
damaged,  we  must  do  everything  possible  to  prevent  radioactive  materials 
from  being  released  into  the  environment. 

From  the  above  discussion,  one  can  see  that  whether  we  adopt  the  nuclear 
safety  concept  or  the  traditional  ship  safety  concept,  the  ultimate  goal  is 
to  protect  human  lives  and  well-being,  as  well  as  properties  and  the  environ¬ 
ment.  But  there  are  two  different  approaches  to  achieve  this  goal:  one 
insists  that  the  safety  of  the  nuclear  reactor  comes  first,  the  other  insists 
that  the  safety  of  the  ship  comes  first. 

The  early  designs  of  nuclear  ships  and  their  propulsion  systems  (e.g.,  the 
Savannah,  the  Otto  Hahn,  and  the  Luau)  basically  adopted  both  safety  concepts 
and  design  guidelines.  Specifically,  the  ship  design  followed  the  guidelines 
and  standards  established  for  ships  powered  by  conventional  fuels,  and  the 
design  of  nuclear  propulsion  system  followed  the  guidelines  and  standards  for 
large  nuclear  power  plants. 

However,  using  two  different  concepts  simultaneously  to  guide  the  design  and 
operation  of  a  nuclear  ship  generally  cannot  produce  the  desired  overall 
safety  standards;  under  certain  conditions  it  may  even  lead  to  handling  errors. 
For  example,  consider  a  scenario  where  the  ship  is  traveling  in  bad  weather  or 
through  a  dangerous  zone  in  the  ocean,  and  at  the  same  time  an  accident  has 
occurred  which  may  cause  damage  to  the  reactor,  using  both  safety  concepts 
in  this  case  will  lead  to  a  chaotic  situation.  If  we  follow  the  rule  that 
reactor  safety  comes  first  and  shut  down  the  reactor,  then  the  ship  will  lose 
power  and  begin  drifting,  thereby  running  the  risk  of  striking  a  reef, 
running  aground  or  even  sinking;  if  we  follow  the  rule  that  ship  safety  comes 
first  and  force  the  reactor  to  continue  operation,  then  we  run  the  risk  of 
damaging  the  reactor  or  causing  leakage  of  radioactive  materials. 

In  recent  years,  a  safety  concept  has  evolved  which  considers  the  ship  and 
its  nuclear  propulsion  system  as  an  integral  unit;  this  concept  is  based  on 
comprehensive  studies  of  the  safety  concept  of  land-based  nuclear  power  plants 
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and  large  ships  and  adopts  the  traditional  approaches  of  handling  the  safety 
problem.  It  treats  the  nuclear  ship  as  a  "floating  nuclear  facility  designed 
to  carry  out  a  special  mission."  This  is  the  "nuclear  ship  safety  concept," 
which  requires  that:  "the  design,  construction  and  operation  of  a  nuclear 
ship  must  meet  the  criterion  that  under  all  possible  conditions,  normal  or 
abnormal,  it  shall  not  release  radioactive  materials  which  may  cause  damages 
to  the  crew  members,  the  general  public,  or  the  environment.  Under  abnormal 
conditions,  the  first  priority  is  to  protect  human  lives. "1  Clearly,  this 
emphasizes  protection  against  radioactivity,  but  it  has  the  implication  that 
ship  safety  comes  first,  and  that  survival  of  the  ship  is  the  most  important 
consideration  in  ensuring  the  ultimate  safety  of  the  reactor  and  preventing 
the  leakage  of  radioactive  materials.  Using  this  concept  to  handle  the 
relationship  between  ship  safety  and  reactor  safety,  the  following  guidelines 
can  be  established: 

1)  The  reactor  is  allowed  to  operate  beyond  its  technical  specifications  or 
to  suffer  limited  damage  if  it  is  judged  to  be  necessary  or  beneficial  for 
the  overall  safety  of  the  nuclear  ship.  Because  it  is  much  safer  for  a  ship 
to  be  powered  by  a  slightly  damaged  (still  controllable)  nuclear  propulsion 
system  than  to  lose  power  and  drift  freely  over  the  ocean. 

2)  If  the  ship  is  already  beyond  the  point  of  salvageability,  then  from  the 
point  of  view  of  preventing  radioactive  leaks,  the  reactor  should  be  shut 
down.  It  is  much  safer  to  have  the  reactor  safely  shut  down  than  to  have  a 
damaged  reactor  if  the  ship  were  to  run  aground  or  sink.  This  implies  that 
the  reactor  is  only  allowed  to  be  shut  down  under  the  following  conditions: 
(a)  the  accident  has  completely  damaged  the  reactor,  thus  power  will  be  lost 
in  any  case;  (b)  even  if  power  is  available  or  the  reactor  can  be  forced  to 
operate  beyond  its  technical  specifications,  the  ship  will  still  sink  or  be 
destroyed  by  other  means. 

III.  Safety  Design  of  Nuclear  Propulsion  System 

1.  Important  Factors  in  the  Safety  Design  of  Nuclear  Propulsion  System 

The  safety  concept  of  nuclear  ship  does  not  imply  that  the  safety  of  a  ship 
reactor  is  not  as  important  as  that  of  a  land-based  nuclear  power  plant  or 
need  not  be  taken  seriously,  neither  does  it  advocate  the  idea  that  the 
reactor  can  be  forced  to  operate  beyond  its  technical  specifications  or  to 
operate  in  a  mode  which  may  be  damaging  to  the  reactor  if  the  ship  is  not 
subject  to  any  danger.  On  the  contrary,  this  concept  requires  that  other 
factors  should  be  considered  in  designing  a  nuclear  propulsion  system; 
specifically,  in  addition  to  a  comprehensive  safety  design  nominally 
required  for  nuclear  power  plants,  one  should  also  consider  the  effects  of 
the  following  factors. 

(1)  Possible  sea  disasters.  They  include  collision,  striking  a  reef,  run¬ 
ning  aground,  water  leakage,  sinking  and  fire.  According  to  the  statistics 
published  by  the  Liverpool  Insurance  Co.  on  sea  disasters  of  ships  above 
20,000  tons  between  1949  and  1959^  (Table  1) ,  the  probability  of  sea 
disaster  that  may  endanger  the  reactor  room  is  as  high  as  7xlO“3. 
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Table  1.  Disasters  That  Occurred  on  Above  20,000-Ton  Ships 


Type  of  disaster 

Number  of 

occurrences 

Number  of  disasters 
that  may  endanger 
the  reactor 

Fire 

9 

8 

Collision 

15 

13 

Running  aground 

19 

11 

Weather 

3 

2 

Mechanical  failure 

5 

0 

Total 

51 

34 

Note:  For  high-speed,  large  ships,  the  probability  of  encounter¬ 
ing  a  disaster  that  may  endanger  the  reactor  is  0.007. 

(a)  The  Effect  of  Collision  and  Striking  a  Reef  and  the  Corresponding 
Protective  Measures 

A  strong  collision  can  damage  the  reactor,  and  also  cause  failure  in  the 
cooling  system  or  other  safety  equipment,  thereby  creating  a  serious  problem 
of  radiation  leakage.  Therefore,  a  nuclear  ship  must  be  equipped  with  a 
special  collision  avoidance  system  and  a  highly  reliable  emergency  cooling 
system;  furthermore,  the  cabin  layout  and  seal  design  must  meet  the  sink- 
proof  requirement . 

There  are  two  types  of  collision  avoidance  measures:  active  collision 
avoidance  and  passive  collision  avoidance.  The  active  measure  involves 
the  installation  of  high-quality  navigation  instruments  and  various  naviga¬ 
tion  systems  to  minimize  the  probability  of  collision  caused  by  human  error, 
the  use  of  highly  reliable  propulsion  system  and  control  system  to  ensure  that 
the  nuclear  ship  will  not  lose  its  mobility  under  any  circumstances,  and  the 
installation  of  emergency  propulsion  system  to  provide  emergency  propulsion 
for  a  sufficiently  long  time  in  case  the  main  propulsion  system  is  damaged. 

The  passive  measure  is  used  to  mitigate  damages  to  the  ship  and  the  reactor 
in  case  of  a  collision.  There  are  two  types  of  commonly  used  passive 
collision  avoidance  structures:  the  deformable  type  structure  and  the 
rigid  frame  (with  shock  reducing  device)  type  structure.  The  former  relies 
on  non-linear  plastic  deformation  to  achieve  shock  reduction,  the  latter 
relies  on  linear  response  to  achieve  shock  reduct ion. ^  In  order  to  minimize 
the  probability  of  failure  of  the  emergency  core  cooling  system  due  to  colli¬ 
sion,  it  has  been  suggested  that  3-4  different  emergency  cooling  systems 
should  be  installed,  each  with  100  percent  capacity;  furtheraore,  they 
should  be  located  in  cabins  which  are  isolated  from  one  another. 

(b)  Effect  of  Fire 

Fire  not  only  threatens  the  ship’s  survivability,  but  may  also  cause  failure 
to  the  reactor  protective  system  and  the  emergency  cooling  system,  and  lead 
to  serious  disaster.  Fire  prevention  on  a  nuclear  ship  should  be  given 
special  attention  because  unlike  a  land-based  power  plant,  external  assistance 
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is  generally  not  available.  Fire  prevention  measures  include  various  tech¬ 
niques  for  preventing  water  loss,  cabin  isolation,  fire  detection  and  predic¬ 
tion  and  a  reliable  fire  extinguishing  system. 

(c)  Effect  of  Water  Entering  in  the  Ship  Cabin 

When  water  enters  in  the  ship  cabin,  the  emergency  power  supply,  the  reactor 
protective  system,  and  various  safety  equipment  may  all  be  adversely 
affected. 

(2)  Loss  of  Mobility 

We  must  fully  realize  the  serious  consequences  that  once  the  nuclear  propul¬ 
sion  system  stops  working,  the  ship  will  lose  its  mobility  which  may  lead  to 
destruction  of  the  ship  as  well  as  the  nuclear  propulsion  system.  For  this 
reason,  the  nuclear  propulsion  system  must  be  designed  with  much  higher 
reliability  than  a  land-based  power  plant.  In  addition,  to  ensure  the  over¬ 
all  safety  of  the  ship  when  it  is  traveling  under  dangerous  conditions,  the 
propulsion  system  may  be  required  to  operate  beyond  its  technical  specifica¬ 
tions  or  to  operate  at  the  risk  of  being  slightly  damaged.  Therefore,  the 
design  safety  margin  of  the  propulsion  system  should  be  sufficiently  high. 

(3)  Effect  of  Other  Factors 

(a)  We  must  realize  that  because  of  the  limited  space  on  a  ship,  it  is  not 
possible  to  install  a  safety  system  as  sophisticated  and  as  large  as  that  of 
a  land-based  power  plant.  Consequently,  we  should  not  and  cannot  use  the 
safety  measures  of  nuclear  power  plants  to  achieve  the  safety  standards  of 
nuclear  ships. 

(b)  The  ship  has  limited  resources  for  handling  an  accident;  in  general  it 
cannot  rely  on  external  assistance. 

(c)  The  living  conditions  on  a  ship  are  far  inferior  than  those  on  land; 
hence  the  mental  conditions  of  ship  reactor  operators  are  generally  not  as 
sound  as  their  counterparts  working  in  a  land-based  power  plant,  and  the 
possibility  of  human-induced  accident  is  greater. 

(d)  Because  the  possible  combinations  of  ship  disasters  and  nuclear  acci¬ 
dents,  the  probability  of  occurrence,  the  resulting  consequences,  the  method 
of  control,  and  the  ultimate  social  risks  are  all  difficult  to  predict,  in 
practice  it  is  impossible  to  implement  a  perfect  emergency  plan.  Therefore, 
the  ultimate  authority  to  allocate  safety  measures  between  the  ship  and  the 
reactor  should  be  reserved  for  the  captain,  who  has  the  best  knowledge  about 
the  nature  and  seriousness  of  the  accident,  the  development  process,  and 
the  possible  social  risks. 

2.  Ways  To  Improve  the  Safety  of  Nuclear  Propulsion  System 

The  "In-Depth  Defense"  method  used  to  prevent  accidents  from  occurring  in 
land-based  reactors  can  also  be  used  in  the  safety  design  of  ship  reactors. 
Generally,  such  a  defense  system  contains  three  layers  of  defense. 
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The  first  line  of  defense  is  to  prevent  accidents  from  occurring  by  mini¬ 
mizing  the  probability  of  malfunctions  and  errors  that  may  cause  accidents. 
The  important  technical  guidelines  include  the  following;  1)  the  design, 
construction  and  management  of  each  system  must  strictly  comply  with  very 
high  standards;  2)  the  design  malfunction  rate  and  failure  rate  of  each 
system  must  be  extremely  low;  3)  all  systems  and  equipment  must  have  high 
degree  of  tolerance  for  operator  errors;  4)  all  materials  and  components 
must  be  thoroughly  tested;  and  5)  all  monitoring  and  control  systems  should 
have  sufficient  margin. 

The  second  line  of  defense  is  to  prevent  a  minor  malfunction  or  error  from 
developing  into  a  major  accident.  The  key  approach  is  to  maintain  safe 
operation  of  the  reactor  or  to  execute  safe  reactor  shutdown  in  case  of  an 
operator  error  or  equipment  malfunction.  The  practical  measures  to  carry 
out  this  mission  Include  the  following;  1)  provide  multi-channel  electric 
power  supply  for  the  nuclear  facility;  2)  install  a  sensitive  detection 
system  that  can  identify  a  damaged  housing  of  the  fuel  elements  or  broken 
pressure  boundary  of  the  cooling  system;  3)  Install  an  automatic  or  manual 
system  that  can  shut  down  the  reactor  based  on  information  obtained  from 
various  detection  systems. 

The  third  line  of  defense  is  to  prevent  radioactive  materials  from  being 
scattered  into  the  environment  once  a  major  accident  has  already  occurred. 

The  key  measures  for  accomplishing  this  task  Include  the  following;  1)  a 
sealed  "safety  container"  for  isolating  the  nuclear  Island,  on  a  ship  there 
is  also  a  safety  isolation  cabin;  2)  an  emergency  cooling  system;  3)  a  heat 
removal  system  in  the  safety  container;  and  4)  a  system  inside  the  safety 
container  for  disposing  radioactive  materials. 

Because  of  the  various  considerations  mentioned  above,  the  relative  emphasis 
for  the  different  lines  of  defense  and  the  actual  implementation  of  the 
safety  concepts  and  measures  for  a  ship  nuclear  propulsion  system  differ 
from  those  of  a  land-based  nuclear  power  plant.  For  example,  because  of  the 
possibility  of  collision  and  water  leakage,  the  probability  of  loss  of  power 
is  much  higher  for  a  nuclear  propulsion  system  than  for  a  land-based  power 
plant;  because  of  the  limited  space  on  a  ship,  it  is  not  possible  to  install 
a  sophisticated  safety  system,  hence  the  probability  of  failure  of  the  system 
is  higher.  Consequently,  we  cannot  rely  on  the  second  and  third  lines  of 
defense  to  ensure  the  propulsion  system’s  safety;  rather,  we  should  strive 
to  strengthen  the  first  line  of  defense,  i.e.,  minimizing  the  probability 
of  accident  occurrence  and  increasing  the  reactor's  own  ability  to  guard 
against  accidents,  in  other  words,  increasing  the  reactor's  inherent  safety. 

The  inherent  safety  of  a  reactor  can  be  improved  by  using  the  following 
approaches;  1)  Using  reliable  fuel  elements  and  increasing  the  thermodynamic 
and  hydrod3mamic  design  margin,  e.g.,  using  low-temperrature  and  low  heat 
capacity  fuel  elements  (increasing  the  surface  area  to  fuel  volume  ratio) , 
lowering  the  surface  thermal  loading  of  the  elements,  and  lowering  the 
irunning  parameters  of  the  heat -carrying  agents.  2)  Using  a  more  reliable, 
reserve-reactivity  control  technique,  so  that  reserve  reactivity  is  intro¬ 
duced  gradually  as  the  fuel  is  being  burned  up.  This  technique  reduces  the 
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use  of  control  rods,  and  therefore  reduces  the  probability  of  accidents  due 
to  control  rod  failure.  The  best  technique  is  to  use  a  solid  diffusion 
type  burnable  poison  for  reactivity  control.  This  technique  not  only  pro¬ 
vides  high  reserve  reactivity,  but  occupies  no  heat-release  surface  area; 
with  the  proper  design,  even  cold  sea  water  can  enter  the  reactor  core 
without  causing  critical  accidents.  3)  Increasing  the  natural  circulation 
of  the  primary  circuit,  thus  reducing  or  eliminating  the  requirement  for 
external  emergency  cooling  system  in  case  of  an  accident.  A  good  design 
will  provide  a  certain  level  of  output  power  for  emergency  propulsion. 

IV.  Safety  Analysis  of  Nuclear  Propulsion  System 

Because  the  mission  and  operating  environment  of  a  nuclear  propulsion 
system  are  different  from  those  of  nuclear  power  plant,  the  operational 
conditions  and  accidents  encountered  by  the  two  systems  are  also  different. 
Therefore,  different  procedures  are  used  for  safety  analysis  and  evaluation. 
Ref.  4  presents  a  method  for  classifying  the  types  of  operational  conditions 
and  gives  a  description  of  each  of  the  conditions.  The  bases  for  classify¬ 
ing  operational  conditions  are  the  frequency  of  accident  occurrence,  the 
possible  consequences,  and  the  degree  of  difficulty  for  recovery.  In  con¬ 
sideration  of  the  unique  features  of  ship  accidents,  the  operational  condi¬ 
tion  IV  is  further  divided  into  two  sub-conditions  IV  A  and  IV  B  (Table  2) . 


Table  2.  Bases  for  Classifying  the  Operational  Conditions  of  Nuclear 
Propulsion  Systems 


Operational 

condition 

I 

II 

III 

IV 

Illustration 

Normal  opera¬ 

Minor  accident 

Major  accident 

Serious  accident 

tion 

Frequency  of 
occurrence 

Continuous  or 
frequent 
(many  times/ 
year) 

Infrequent 

(once/year) 

Rare 

Extremely  rare 
(lO^^/year) 

Consequence 

classifi¬ 

1 

2 

3 

4 

cation 

Consequence 

treatment 

Possible 
reduction 
in  power 

Reduction  in 
power  or 
reactor 
shutdovm 

Fast  shutdown 
plus  residual 
heat  removal 

Fast  shutdown 
plus  emergency 
cooling  or 
residual  heat 
removal 

Restarting 

Immediate 

After  repair 

After  dock 
repair 

Not  possible 

Release  of 

Less  than  normal  limiting  value 

Less  than 

radio-  accident  limit- 

activity  value 

Note:  Operational  condition  IV  can  be  divided  into  two  categories:  IV  A  is 
the  case  with  design  power  supply;  IV  B  is  the  case  without  design  power  supply. 
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The  definition  of  consequence  classification  in  Table  2  is  as  follows: 

Consequence  level  1:  the  level  of  radioactivity  of  the  ship  and  its  steam 
supply  system  is  within  the  specified  limits  for  normal  operation. 

Consequence  level  2:  it  is  caused  by  an  accident,  but  operation  of  the 
ship  is  not  seriously  affected,  and  the  level  of  radioactivity  is  lower 
than  the  allowable  radiation  level  for  reactor  technicians . 

Consequence  level  3:  it  is  caused  by  single  or  multiple  accidents;  the 
ship  structure  or  propulsion  system  is  damaged,  and  there  are  human  injuries; 
the  ship  is  partially  put  out  of  service  and  requires  external  assistance, 
but  the  level  of  radioactivity  does  not  exceed  the  allowable  radiation 
dosage  for  reactor  technicians. 

Consequence  level  4:  it  is  caused  by  single  or  multiple  accidents;  it  may 
lead  to  death  of  personnel  and  sinking  of  ship,  but  the  level  of  radio¬ 
activity  in  the  hull  and  on  the  deck  above  the  reactor  does  not  exceed  twice 
the  allowable  radiation  dosage  for  reactor  technicians. 

Table  3  presents  detailed  descriptions  of  the  four  classes  of  operational 
conditions.  The  purpose  of  safety  analysis  is  to  provide  quantitative 
answers  to  the  consequences  of  these  operational  conditions  and  accidents. 

Table  3.  Examples  of  Classification  of  Operational  Conditions  of  a  Nuclear 
Propulsion  System 

Classification  of 
operational 

conditions  _ Name  of  operational  conditions _ 

I  Starting,  stopping,  low-power  operation,  full-power  opera- 

(normal)  tlon,  variable  load,  regular  maintenance,  bad  weather, 

loading  and  unloading  cargos,  refueling,  and  overload 
operation,  which  has  already  been  considered  in  the  design 
and  operation  procedures. 

II  Component  failures  or  errors  which  lead  to  loss  of  power 

(infrequent  or  loss  of  mobility;  e.g.,  sudden  shutdown  of  generator, 

minor  accidents)  sudden  shutdown  of  a  main  engine,  sudden  shutdown  of  a 

condenser,  cut-off  of  sea  water  cooling  unit,  cut-off  of 
water  supply,  shut-off  of  main  steam  valve,  malfunction 
in  the  main  electric  system,  malfunction  in  the  steering 
mechanism  or  anchoring  mechanism;  accidental  starting  of 
main  coolant  pump  or  water  pump,  primary  operator  error, 
malfunction  in  the  reactor  protective  system  or  safety 
system,  error  in  operating  the  pressure  relief  valve  or 
safety  valve  of  the  primary  system,  accidental  withdrawal 
of  control  rod,  temporary  power  outage,  and  minor  leak 
in  the  primary  system. 


[Continued  on  following  page] 


Classification  of 
operational 

conditions  _ Name  of  operational  conditions _ 

III  Leaks  in  the  primary  circuit  pressure  boundary  which  does 

(rare  occurrence)  not  cause  pressure  drop  in  the  primary  circuit  but 

requires  activation  of  the  emergency  cooling  system  and 
safety  shield  for  Isolation;  leaks  in  the  radioactive 
waste  treatment  system  and  storage  system;  leakage  of 
primary  water  into  the  secondary  system,  which  requires 
reactor  shutdown;  malfunction  in  the  control  rod  drive 
mechanism;  jammed  control  rod;  malfunction  in  the  forced 
circulation  system  of  the  primary  circuit;  malfunction 
in  the  water  supply  system  of  the  secondary  circuit; 
exhaust  from  the  primary  system;  malfunction  in  the 
coolant  supply  system  in  the  docks;  striking  a  reef  with 
no  damage  to  the  ship  integrity  and  no  loss  of  coolant; 
collision  which  causes  water  leakage  between  two  neigh¬ 
boring  cabins;  fire  or  explosion  with  no  serious  damage 
to  the  reactor  room;  fire  in  the  machine  room  or  opera¬ 
tions  room;  unusually  bad  weather  along  the  designated 
course. 

Loss  of  water  accident;  destruction  of  the  waste  storage 
tank;  striking  a  reef,  which  may  cause  tilting  of  the 
ship  and  intermittent  loss  of  coolant,  but  still  maintain 
Integrity;  extremely  bad  weather;  ejection  of  control  rod; 
running  aground,  which  can  cause  partial  damage  to  ships 
with  large,  double-layered  bottoms  or  ships  with  long 
hulls;  collision  which  causes  fire  or  explosion. 

IV  B  Capsizing;  ship  sinking  in  deep  or  shallow  waters,  caus- 

(rare  occurrence  ing  failure  in  emergency  power  supply  system, 
of  extreme  acci¬ 
dents) 

V.  Concluding  Remarks 

We  see  from  the  previous  discussions  that  nuclear  ship  is  a  unique  system 
which  has  the  technical  and  safety  problems  of  both  the  shipping  industry  and 
the  nuclear  power  industry.  In  view  of  the  public’s  concern  about  nuclear 
safety,  this  system  should  be  considered  as  a  “floating  nuclear  facility  with 
a  special  mission."  The  ultimate  goal  of  safety  measures  should  be  to  pre¬ 
vent  radiation  damages,  but  we  must  assign  first  priority  to  the  survival  of 
the  ship  and  establish  an  overall  nuclear  ship  safety  concept. 

In  case  there  is  a  conflict  between  reactor  safety  and  ship  safety,  the 
overall  system  safety  should  be  the  ultimate  guideline.  Unless  the  reactor 
is  completely  damaged  or  it  will  destruct  even  if  the  ship  is  still  under 
power,  the  reactor  should  continue  operation. 


IV  A 

(rare  occurrence 
of  extreme  acci¬ 
dents) 
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Because  a  ship  disaster  and  a  reactor  failure  affect  one  another,  an  accident 
in  either  will  lead  to  very  grave  consequences.  Therefore,  the  safety 
requirement  for  a  nuclear  ship  and  its  propulsion  system  should  be  higher 
than  that  for  a  ship  with  conventional  power  or  a  nuclear  power  plant.  The 
key  approach  to  improving  the  safety  of  a  nuclear  propulsion  system  is  to 
improve  the  safety  of  the  reactor  itself.  The  safety  procedures  of  a  nuclear 
propulsion  system  cannot  simply  be  copied  from  those  used  for  nuclear  power 
plants;  they  should  be  determined  by  careful  study  of  the  particular  ship 
application. 

Since  one  cannot  predict  every  type  of  accidents  and  establish  protective 
measures  to  guard  against  all  possible  combinations  of  potential  problems, 
the  ultimate  authority  to  allocate  safety  measures  between  the  ship  and  the 
reactor  should  be  reserved  for  the  captain,  who  has  the  best  understanding 
of  the  seriousness  of  the  accident,  the  development  process  and  the  possible 
consequences. 

It  should  be  pointed  out  that  the  analyses  presented  in  this  article  are 
primarily  aimed  at  nuclear-powered  commercial  ships.  For  a  military  vessel, 
the  effect  of  mission  safety  (which  is  often  more  important  than  the  safety 
of  the  propulsion  system)  should  also  be  analyzed;  but  the  final  result  is 
similar. 
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PULSE  CALIBRATION  TECHNIQUE  OF  X-RAY  DETECTOR 

Beijing  WULI  XUEBAO  [ACTA  PHYSICA  SINICAl  in  Chinese  Vol  35  No  7,  Jul  86 
pp  864-873 

[English  abstract  of  article  by  Sun  Jingwen  [1327  2529  2429]  of  the  Southwest 
Institute  of  Fluid  Physics,  Chengdu] 

[Text]  A  pulsed  CuK  X-ray  source  with  flux  intensity  10^®-10^^  K  X-rays  per  sr 
per  s  has  been  produced  using  the  high  current  electron  beam  technique.  It 
has  been  employed  to  calibrate  the  response  of  the  silicon  PIN  diode  detector 
to  pulsed  X-ray  photons.  An  absolute  X-ray  monitor,  a  PIO  gas  pulse  ionization 
chamber,  has  been  used  as  a  standard  for  the  pulsed  X— ray  flux  density.  The 
automatic  pulsed  charge  measuring  instrument  is  controlled  and  corrected  in 
real  time  by  a  microcomputer.  The  measured  uncertainty  of  the  flux  density  is 
±5  percent.  This  chamber  is  suitable  for  the  energy  flux  rate  range  of 
4  X  10-“^  -  2  X  10^  W/cm^  and  the  photon  energy  range  of  1.5  -  10  keV.  The 
calibrated  accuracy  of  the  detector  is  ±7.0  percent.  It  is  found  that  the 
pulsed  sensitivity  of  the  silicon  PIN  diode  is  about  30  percent  higher  than  the 
sensitivity  calibrated  by  steady  X-ray  beams. 


94 


STABLE  REGION  OF  TOKAMAK  DEVICE  FEEDBACK  CONTROL  (I).  THEORETICAL 

Beijing  WULI  XUEBAO  [ACTA  PHYSICA  SINICA]  in  Chinese  Vol  35  No  7,  Jul  86 
pp  874-881 

[English  abstract  of  article  by  Zheng  Shaobai  [6774  1421  4101]  and  Shen 
Zhongqing  [3088  0112  0615]  of  the  Institute  of  Physics,  Chinese  Academy  of 
Sciences] 

[Text]  In  this  paper,  the  evolution  equations  of  feedback  control  of  plasma 
current  and  horizontal  displacement  simultaneously  in  a  tokamak  device  are 
presented.  Under  the  assumption  of  slight  disturbance  and  linear  approxima¬ 
tion,  the  stability  criterion  and  the  stable  region  of  the  CT-6B  tokamak 
feedback  control  are  given. 
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STE^Y-STATE  ELECTRIC  FIELD  STRUCTURE  AND  DENSITY  PROFILE  IN  A  PLASMA 
IRRADIATED  BY  DOUBLE-FREQUENCY  LASER 

Beijing  WULI  XUEBAO  [ACTA  PHYSICA  SINICA]  in  Chinese  Vol  35  No  7,  Jul  86 

pp  882-888 

[English  abstract  of  article  by  Zhu  Shitong  [2612  5535  6639]  of  Shanghai 
Institute  of  Optics  and  Fine  Mechanics,  Chinese  Academy  of  Sciences;  and 
Shen  Wenda  [3088  2429  66711  of  the  Department  of  Physics,  Shanghai  University 
of  Science  and  Technology] 

[Text]  The  analytic  expressions  for  steady-state  electric  field  intensity 
and  local  scale  length  in  a  plasma  irradiated  by  a  double-frequency  laser  are 
derived.  It  is  shown  that  the  ponderomotive  force  due  to  the  double-frequency 
laser  field  would  lead  to  the  local  density  dip  near  the  critical  surface  and 
the  anomalous  increase  and  fluctuation  of  the  critical  scale  length  in  some 
region  of  field  intensity  value. 
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COMPARISON  OF  TWO  COLLISION  MODELS  IN  PLASMA 

Beijing  WULI  XUEBAO  [ACTA  PHYSICA  SINICA]  in  Chinese  Vol  35  No  7,  Jul  86 
pp  947-952 

[English  abstract  of  article  by  Zhang  Chengfu  [1728  2110  4395]  of  the 
Department  of  Physics,  Beijing  University] 

[Text]  In  this  paper,  the  effects  of  collision  on  the  dispersion  equation  of 
plasma  are  discussed  for  two  kinds  of  collision  models  (the  BGK  model  and  the 
LB  model).  Analysis  shows  that  the  main  effect. of  collision  is  to  change  the 
usual  dispersion  function  Z(5)  into  a  generalized  one  W(5,U).  The  properties 
of  the  W  function  have  been  analyzed  and  compared  for  the  two  models.  The 
effects  of  collision  on  the  low  frequency  drift  wave  (local  mode)  and  ion- 
acoustic  wave  are  discussed.  It  is  shown  that  the  LB  model  is  not  only  more 
reasonable  physically  than  the  BGK  model,  but  also  simpler  in  calculation  and 
expression. 
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ELECTRONMICROSCOPIC  IMMUNOFERRITIN  TECHNIQUE  AND  ITS  APPLICATION 

Shijiazhuang  ZHONGHUA  WULI  YIXUE  ZAZHI  [CHINESE  JOURNAL  OF  PHYSICAL 
MEDICINE]  in  Chinese  Vol  8  No  2,  25  Jun  86  pp  88-90 

[English  abstract  of  article  by  Yu  Chengzhi  [0151  3397  0037]  of  the  Institute 
of  Microbiology  and  Epidemiology,  Academy  of  Military  Medical  Sciences, 
Beijing] 

[Text]  Immunoferritin  indirect  staining  of  ultrathin  sections  in  electron 
microscopy  was  established  on  BHK-13  cells  infected  with  Sindbis  virus.  The 
author  observed  HFRS  virus  particles  isolated  in  Vero  E-6  cells  from  a 
patient  with  HFRS  using  this  method.  In  addition,  this  method  has  been 
successfully  applied  to  the  observation  of  the  surface  structure  of  the 
Jinghe  strain  of  spotted  fever  rickettsiae  in  the  yolk  sac  section  which 
could  not  be  detected  with  ordinary  electron  microscopy.  (Paper  received 
April  1985.) 
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RATE-CONTROLLED  LIQUID  NITROGEN  COOLING  EQUIPMENT 


Shijiazhuang  ZHONGHUA  WULI  YIXUE  ZAZHI  [CHINESE  JOURNAL  OF  PHYSICAL 
MEDICINE]  in  Chinese  Vol  8  No  2,  25  Jun  86  pp  93-96 

[English  abstract  of  article  by  Wan  Guoying  [8001  0948  2019],  et  al.,  of  the 
Academy  of  Military  Medical  Sciences] 


[Text]  A  rate-controlled  cooling  apparatus  has  been  reported.  It  can  be 
used  for  freezing  bone  marrow  cells,  blood  cells,  hematopoietic  cells  of 
embryonic  liver,  cornea  and  skin,  etc.  After  cryopreservation  of  cells, 
the  cell  viability  is  more  than  80  percent.  The  equipment  consists  of  three 
parts:  (1)  a  program-controller;  (2)  a  system  for  storage  and  transfer  of 

LNa;  (3)  a  freezing  chamber.  There  are  seven  cooling  rates  in  the  equipment. 
Start-freezing  temperature  and  end-freezing  temperature  can  be  adjusted 
continually  and  chosen  freely.  The  cooling  rate  is  controlled  automatically 
during  the  freezing  process.  The  temperature  of  the  freezing  sample  and  in 
the  freezing  chamber  is  shown  digitally.  The  equipment  is  applicable  to 
cryobiological  research  and  cryopreservation  of  biological  materials. 

(Paper  received  1  February  1985.) 


1.  Mazur  P  and  Schmidt  J  J.  Interactions  of  cool- 
>ng  Telocit7,  temperature, and  warming  velocity 
on  the  survival  of  frozen  and  thawed  yeast. 
Cryobiology  1868,  5(l)il. 

2.  •  Mazur  P,  et  al.  A  two-factor  hypothesis  of 

freezing  injury.  Exp  Cel]  Res  1972,  7(2), S45. 

3.  Leibo  S  P,  et  al.  Factors  affecting  survival  of 
mouse  embryos  during  freezing  and  thawing.Exp 
Cell  Res  1974,  89(I)i79. 

4.  Cavins  J  A,  et  al.  Recovery  of  lethally  irradi¬ 
ated  dogs  following  infusion  of  autologous  mar¬ 
row  stored  at  low  temperature  in  dimethyl-sul- 
phoxide.  Blood  1962,  20(6), 730. 


REFERENCES 

5.  Leibo  S  P  and  Mazur  P.  The  role  of  cooling 
rates  in  low- temperature  preservation.  Cryobio¬ 
logy  1971,  8(5), 447. 

6,  Kelley  F  D,  et  al.  Controlled-rate  liquid  N** 
fflicrowave  biological  freeze-thaw  device.  Cryo¬ 
biology  1982,  19(4),372. 

)?tt:^I984,  12,21. 


9717 

CSO:  4009/1005 


99 


Nuclear  Medicine 


POSITRON  EMISSION  TOMOGRAPHY — ^A  NEW  NUCLEAR  MEDICINE  TECHNIQUE 

Shanghai  ZIRAN  ZAZHI  [NATURE  JOURNAL]  in  Chinese  Vol  9  No  4,  Apr  86  pp  269-271 

[Article  by  Yu  Renkang  [0205  0088  1660],  Jiaotong  University,  Shanghai] 

[Abstract]  PET  stands  for  positron  emission  tomography,  which  merges  the 
biochemical  evaluation  of  pathology  in  early  nuclear  medicine  methods  with 
the  reconstructed  precise  positioning  using  computer  graphics.  There  are 
many  unique  features  in  medical  diagnosis  and  physiological  research j  there¬ 
fore,  research  on  its  clinical  applications  is  vigorously  under  way  in  West 
Europe  and  the  United  States.  Research  on  PET  will  widely  develop  in  China 
with  the  radioactive-isotope  technique  and  computer  tomography.  This  paper 
presents  the  fundamental  principle  of  PET  and  its  tracers.  The  PET  system 
requires  a  medical  cyclotron  and  the  corresponding  chemical  equipment,  a 
detector  array  for  receiving  gamma  photons,  and  a  computer.  Cooperation  is 
required  among  personnel  specialized  in  physics,  chemistry,  computer, 
physiology,  and  medicine.  More  than  40  PET  research  centers  have  been 
established  throughout  the  world.  The  author  saw  a  PET  model  IV  system 
at  the  Medical  College  of  the  University  of  Washington  in  the  United  States. 
The  system  was  jointly  developed  by  the  university  and  a  private  company. 

The  PET  system  can  be  used  to  examine  the  human  head  and  trunk.  In  the  newer- 
generatlon  PET  system,  an  important  research  direction  is  the  addition  of 
photon  transition-time  data  in  the  graphic  reconstruction  process.  This 
requires  an  accurate  measurement  of  time  difference  of  two  trapped  photons 
between  their  arrivals  at  the  receiver  in  order  to  enhance  the  signal-to- 
noise  ratio  and  the  graph's  resolving  power.  PET's  development  prospects  are 
promising  in  biology  and  medicine.  Five  figures  show  gamma  photons  picked  up 
by  a  ring-shaped  detector,  a  positron  transmission  source  at  the  intersection 
of  a  group  of  straight  lines,  a  projection-recording  radioactivity  distribu¬ 
tion,  wave  filtration,  and  the  reconstruction  of  radioactive  distribution 
with  an  inverse  projection. 
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STUDY  OF  THE  CHARACTERISTICS  OF  PATHS  IN  AtiF*  FORMATION  PROCESSES* 

Shanghai  GUANGXUE  XUEBAO  [ACTA  OPTICA  SINICA]  in  Chinese  Vol  6  No  8,  Aug  86 
pp  679-687 

[English  abstract  of  article  by  Gu  Zhiyu  [7357  0037  3768],  et  al.,  of  Anhui 
Institute  of  Optics  and  Fine  Mechanics,  Chinese  Academy  of  Sciences,  Hefei] 

[Text]  The  kinetic  processes  in  Ar/He/Fa  mixtures  excited  by  fast  discharge 
have  been  analyzed.  A  mathematical  description,  which  could  reveal  the 
characteristics  of  paths  in  the  kinetic  processes,  is  given  based  on  a  sim¬ 
plified  kinetic  model.  Characteristics  of  paths  in  ArjF*  formation  processes 
are  studied  by  comparing  the  experimental  data  with  the  results  of  numerical 
calculations.  (Received  28  February  1985;  revised  31  January  1986.) 

*  Supported  by  the  Science  Fund  of  the  Chinese  Academy  of  Sciences . 
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OFF-RESONANT -PUMPING  AC  STARK  EFFECT  IN  OPTICALLY  PUMPED  FIR  LASER 

Shanghai  GUANGXUE  XUEBAO  [ACTA  OPTICA  SINICA]  in  Chinese  Vol  6  No  8,  Aug  86 
pp  688-694 

[English  abstract  of  article  by  Lin  Yikun  [2651  6318  1024],  et  al.,  of  the 
Department  of  Electronics,  Zhongshan  University,  Guangzhou] 

[Text]  The  interactions  among  R35,  T31,  and  S54  processes  in  a  six-level  system 
have  been  studied.  Near  the  cross  point  of  the  resonant  curves  of  these 
processes  in  the  multiphoton  resonance  diagram  of  the  six-level  system,  such 
interactions  should  be  very  strong.  We  found  theoretically  that  in  the  case 
of  off-resonant-pumping,  the  interactions  among  these  multiphoton  processes 
should  also  lead  to  splitting  of  the  laser  spectral  line,  which  is  called 
the  of f“re sonant —pumping  AC  Stark  effect.  Because  the  influences  of  the 
saturation  effect  and  electric  dipole  transition  enhancement  effect  were  not 
the  same  for  different  multiphoton  processes,  contributions  of  different 
multiphoton  processes  to  the  off-resonant  AC  Stark  effect  could  be  distin¬ 
guished  by  the  injection  of  the  FIR  signal  into  the  system.  (Received 
7  October  1985;  revised  4  February  1986.) 
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NEW  METHOD  FOR  ON-RESONMT  TWO-PHOTON  ABSORPTION  SPECTROSCOPY 

Shanghai  GUANGXUE  XUEBAO  [ACTA  OPTICA  SINICA]  in  Chinese  Vol  6  No  8,  Aug  86 
pp  708-713 

[English  abstract  of  article  by  Zhu  Zhenhe  [2612  2182  0735]  of  the  Department 
of  Mathematics  and  Physics,  Central  Nationalities  Institute,  Beijing] 

[Text]  A  new  experimental  method  for  on-resonant  two-photon  absorption 
spectroscopy  is  presented  in  this  paper.  The  experimental  set-up  is  the  same 
as  that  of  the  two-photon  fluorescence  experiment  used  in  measuring  durations 
of  ultrashort  laser  pulses.  A  picosecond  laser  pulse  is  divided  into  two 
beams  which  then  travel  in  opposite  directions  in  the  sample.  The  signal 
intensity  measured  along  the  laser  trace  will  be  a  wide  plateau  overlapped  by 
a  narrow  peak.  The  ratio  of  the  on-resonant  two-photon  absorption  section 
to  the  two-step  single-photon  absorption  section  can  be  deduced  frcxn  the 
contrast  ratio  of  the  peak  intensity  to  the  shoulder  intensity.  The 
theoretical  calculation  shows  that  this  method  is  valid  for  various  cases. 
(Received  2  August  1985;  revised  13  February  1986.) 
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GROWTH  AND  CHARACTERISTICS  OF  ALEXANDRITE  CRYSTALS 

Shanghai  GUANGXUE  XUEBAO  [ACTA  OPTICA  SINICA]  in  Chinese  Vol  6  No  8,  Aug  86 
pp  734-743 

[English  abstract  of  article  by  Ma  Xiaoshan  [7456  4562  1472],  et  al.,  of 
Shanghai  Institute  of  Optics  and  Fine  Mechanics,  Chinese  Academy  of  Sciences] 


[Text]  The  opaque  sample  and  transparent  alexandrite  crystals  have  been 
analyzed  by  different  techniques,  such  as  the  fluorospectrometer ,  X-ray  fluoro- 
spectrometer  and  ESCA.  The  measured  data  support  the  authors*  previous  sug¬ 
gestion  concerning  the  formation  mechanism  of  Lienard's  *'melt  aging."  It 
follows  that  the  transparent  crystal  can  be  grown  with  good  reproducibility. 


Inclusions  and  dislocations  in  the  crystals  have  also  been  investigated, 
leading  to  the  growing  of  high  quality  alexandrite  crystals.  The  optical 
homogenicity  of  laser  rods  is  0.25X/in  to  0.40X/in  and  is  free  |rom  scattering 
centers.  The  laser  threshold  is  34. IJ.  The  line  width  is  O.OIA  after  Q-switch 
with  a  color  center  crystal . 

The  annealing  effect  pressure  effect  of  Rj  and  Rg  fluoroline  and  striations 
in  alexandrite  crystals  have  also  been  studied.  (Paper  received  19  December 
1985;  revised  28  February  1986.) 
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A  PICOSECOND  GaAs  PHOTOCONDUCT IVE  SWITCH 

Shanghai  GUANGXUE  XUEBAO  [ACTA  OPTICA  SINICAJ  in  Chinese  Vol  6  No  8,  Aug  86 
pp  744-747 

[English  abstract  of  article  by  Zhu  Shidong  [2612  0013  2767] ,  et  al.,  of 
Changchun  Institute  of  Optics  and  Fine  Mechanics,  Chinese  Academy  of  Sciences] 

[Text]  The  development  of  a  picosecond  GaAs  photoconductive  switch  is 
reported.  The  switch  consists  of  a  microstrip  transmission  line  formed  on 
Cr-doped  semi-insulating  GaAs  with  a  small  gap  (~20  pm)  in  the  stripline.  As 
the  source  of  optical  pulses  used  for  preliminary  evaluation  of  the  switch, 
a  synchronously  pumped  Rh6G  dye  laser  was  used.  The  rise  time  of  the  signal 
observed  on  a  sampling  oscilloscope  was  120  ps  and  the  FWHM  was  230  ps.  We. 
measured  the  dependence  of  peak  voltage  of  the  switch  electric  pulse  on  the 
energy  of  the  incident  laser  pulse  for  a  bias  voltage  of  30  V  and  incident 
energy  range  of  38  pJ  ~  320  pj.  (Received  9  January  1986;  revised 
13  February  1986.) 
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APPLICATION  OF  MONTE  CARLO  METHOD  IN  RADIATIVE  HEAT  TRANSFER 

Shanghai  ZIRAN  ZAZHI  [NATURE  JOURNAL]  in  Chinese  Vol  9  No  4,  Apr  86  pp  250-254 

[Article  by  Bian  Bohui  [0593  0130  4940]  and  Li  Hao  [2621  3185]  of  Qinghua 
University] 

[Abstract]  The  Monte  Carlo  method  relies  on  random  variable  statistical 
testing  to  solve  problems  in  mathematics,  physics,  and  engineering  tech¬ 
nology.  For  nearly  100  years  this  random  simulation  method  has  a  history  of 
use  in  approximate  numerical  calculations.  In  its  early  history,  however, 
there  were  limitations  in  simulation-testing  tools;  relatively  few  practical 
problems  could  be  solved.  In  the  past  2  decades,  what  with  the  rapid  progress 
in  computer  technology,  propitious  conditions  were  provided  for  the  Monte  Carlo 
method.  Gradually,  the  method  is  being  applied  to  the  problem  of  calculating 
radiative  heat  transfer,  as  illustrated  in  five  figures.  To  a  great  extent, 
the  method's  accuracy  depends  on  the  number  of  sampling  values  of  the  random 
variables  in  the  simulation.  The  more  the  sampling  values  in  the  random 
variables,  the  higher  the  accuracy.  Therefore,  the  method  demands  a  large 
storage  capacity  in  the  computer  used;  this  is  a  drawback  of  the  method. 
However,  the  result  accuracy  in  random  simulation  is  not  related  to  the 
number  of  dimensions;  thus,  the  method  is  suitable  to  calculating  multi¬ 
dimensional  and  complex  problems.  The  basic  features  of  the  method  are  as 
follows:  First,  the  surface  radiative  heat  energy  is  considered  as  composed 
of  a  large  number  of  energy  beams  (possessing  the  same  energy  for  each  beam) 
emitted  from  the  surface.  The  random  emission  process  of  each  energy  beam  is 
traced  until  its  absorption  by  the  surface  or  its  emission  into  the  ambient 
space.  Finally,  the  number  of  energy  beams  absorbed  by  each  surface  is 
statistically  determined;  this  is  the  radiative  heat  transfer  between  surfaces. 
According  to  References  2  through  4,  the  Monte  Carlo  method  can  also  be  used 
in  calculating  gas  radiation.  One  table  shows  a  comparison  between  radiation 
transfer  coefficients  and  radiation  angle  coefficients;  in  the  calculation, 
the  total  number  of  energy  beams  is  10,000. 
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